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FOREWORD 


This  report,  OSURF  Report  Number  2183-2,  was  prepared  by 
The  ElectroScience  Laboratory,  Department  of  Electrical  Engineering, 
The  Ohio  State  University  at  Columbus,  Ohio.  Research  was  conducted 
under  Contract  AF  19(628)-5929.  Dr.  John  K.  Schindler,  CRDG  of  the 
Air  Force  Cambridge  Research  Laboratories  at  Bedford,  Massachusetts 
was  the  Program  Monitor  for  this  research. 


ABSTRACT 


In  this  report j  jthe  exact  solution  for  the  electromagnetic  field 
diffracted  by  a  perfectly  conducting  plane  angular  sector  is  determined. 
The  problem,  is  a  three-dimensional  vector  problem  and  the  solution  is 
presented  In  the  form  of  a  dyadic  Green's  function,  which  is  the  roost 
general  form  of  solution  possible.  Thus  the  vector  fields  as  well  as 
the  current  on  the  sector  may  be  determined  for  any  given  source  exci¬ 
tation. 

The  comer  angle  of  the  plane  angular  sector  is  arbitraiy,  varying 
between  zero  and  tt.  Special  cases  of  the  plane  angular  sector  are  the 
half  plane  and  the  quarter  plane.  To  find  the  fields  for  larger  angles, 
such  as  at  the  corner  of  an  aperture,  Babinet's  principle  can  be  used. 

The  dyadic  Green's  function  is  composed  of  vector  wave  functions, 
which  in  turn  are  composed  of  scalar  wave  functions,  The  problem  is 
solved  in  a  sphero-conal  coordinate  system.  In  this  system,  the  plane 
angular  sector  is  one  of  the  coordinate  surfaces,  so  that  the  separation 
of  variables  technique  is  used  to  find  the  scalar  wave  functions.  They 
consist  of  spherical  Bessel's  functions  and  Lame  functions.  The  Lame 
functions  are  solutions  of  two  coupled  differential  equations.  These 
equations  are  solved  for  the  special  case  of  a  quarter  plane  scatterer. 
The  first  192  eigenvalues  and  eigenfunctions  are  computed  and  tabulated. 

The  fields  and  currents  close  to  the  tip  of  the  quarter  plane  are 
presented.  These  fields  and  currents  have  been  the  subject  of  much 
conjecture  by  several  authors.  It  is  shown  that  the  dominant  field 
behaves  as  where  the  lowest  value  of  the  eigenvalue  v  is  0.296. 

The  far  fields  for  infinitesimal  dipole  sources  very  close  to  the  tip  are 
also  determined  and  several  patterns  are  presented. 

As  with  any  exact  solution  of  a  complex  problem,  the  results  are 
not  simple.  It  is  felt,  however,  that  the  exact  solution  obtained  here 
lays  the  foundation  for  subsequent  work.  For  instance,  it  should  be 
possible  to  use  this  work  to  determine  an  asymptotic  approximation  and 
thvis  derive  a  "diffraction  coefficient"  for  the  tip.  This  would  be  very 
useful  in  Keller's  "Geometrical  Theory  of  Diffraction."  Without  any 
additional  work,  there  are  sufficient  numerical  results  presented  in  this 
report  to  determine  the  fields  within  approximately  one  wavelength  of 
the  tip  for  any  source,  or  the  fields  everywhere  for  a  source  within  one 
wavelength  of  the  tip. 


TABLE  OF  CONTENTS 


■#- 


Chapter  Page 


I 

INTRODUCTION 

1 

II 

COORDINATE  SYSTEM 

7 

III 

SOLUTION  OF  THE  SCALAR  WAVE  EQUATION 

20 

IV 

SOLUTION  OF  THE  VECTOR  WAVE  EQUATION 

74 

V 

DETERMINATION  OF  THE  DYADIC  GREEN'S 

FUNCTION 

79 

VI 

DISCUSSION  OF  THE  FIELDS  AND  CURRENTS 

FOR  SOME  SPECIAL  CASES 

100 

VII 

CONCLUDING  REMARKS 

144 

APPENDIX 

A 

THE  SELF-ADJOINT  AND  POSITIVE-DEFINITE 
PROPERTY  OF  THE  TWO-DIMENSIONAL 
STURM-LIOUVILLE  OPERATOR 

146 

B 

DETERMINATION  OF  THE  EIGENVALUES 

151 

C 

DETERMINATION  OF  THE  EIGENFUNCTIONS 

173 

D 

VARIATIONAL  METHOD  FOR  DETERMINING 
EIGENVALUES  AND  EIGENFUNCTIONS 

184 

REFERENCES 

189 

iv 


CHAPTER  I 


INTRODUCTION 

The  primary  purpose  of  this  work  is  to  determine  the  exact  solution 
for  the  electromagnetic  field  scattered  by  a  perfectly  conducting  plane 
angular  sector.  A  plane  angular  sector  is  the  section  of  an  infinite 
plane  bounded  by  two  interesting  straight  edges  which  terminate  at  a 
corner.  This  is  a  three-dimensional  vector  problem,  so  the  final  solu¬ 
tion  is  in  the  form  of  a  dyadic  Green’s  function.  The  physical  config¬ 
uration  is  shown  in  Figure  1.  The  corner  angle  is  arbitrary  in  the 
range  zero  to  n.  For  larger  corner  angles,  such  as  the  corner  of  an 
aperture,  Babinet’s  principle  can  be  used  to  obtain  the  fields  dif¬ 
fracted  by  the  coit^jlementary  structure. 

Diffraction  by  objects  with  edges,  corners,  tips,  etc.,  has  occu¬ 
pied  the  attention  of  several  authors.'*^'  The  classical  case  is  the  half¬ 
plane,  which  was  solved  by  Sommerfeld,  and  has  since  been  studied  by 
many  authors.*  The  problem  considered  here  is  more  general,  the  half- 
plane  being  a  special  case  of  the  plane  angular  sector.  The  solution 
near  the  corner,  or  tip,  is  of  primary  interest  in  the  angular  sector 
problem.  A  half -plane  edge  is  a  discontinuity  in  the  scattering  sur¬ 
face  since  the  normal  to  the  surface  cannot  be  defined.  A  tangent  can 

♦ 

be  defined,  however.  The  corner  of  a  plane  angilLar  sector  is  a  "double" 
discontinuity  sitxce  even  the  tangent  cannot  be  defined  there. 


*For  extensive  bibliographies,  see  Reference  17,  Ch.  12,  and  Reference  19. 
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Several  authors  [1>  2,  3,  4,  5]  have  considered  the  field  varia¬ 
tion  at  this  "double"  discontinuity  and  have  made  con^jectures  based  on 
approximations  and  physical  reasoning.  Two  of  these  cases  are  compared 
with  the  exact  solution  in  Chapter  6.  The  effect  of  the  corner  on  the 
far  field  is  ailso  of  interest.  By  using  the  well  kno\'m  asymptotic 
solution  of  the  edge  diffraction  problem,  and  the  reflection  from  a 
plane,  it  should  be  possible  to  identify  a  "diffraction  coefficient" 
for  the  tip.  This  has  not  been  done  in  this  work;  however,  it  is  felt 
th8.t  the  exact  solution  presented  here  is  a  good  starting  point  for 
such  an  endeavor. 

As  far  as  the  authors  could  determine,  there  is  no  published  work 
on  the  electromagnetic  diffraction  by  a  plane  angular  sector.  The 
scalar  problem  has  been  solved  by  Kraus  [1],  Radlow  [2],  and  Kraus  and 
Levine  [6].  Radlow  has  considered  the  problem  of  the  diffraction  of  a 
scalar  plane  wave  by  a  quarter  plane.  He  determined  a  two  variable 
integral  representation  of  the  scattered  field.  Using  a  generalization, 
or  extension,  of  the  classical  one  variable  Wiener-Hopf  method,  he  then 
found  the  transform  solution  that  forces  the  total  field  to  zero  on  the 
quarter  plane.  It  should  be  emphasized  that  his  scalar  solution  is  not 
applicable  to  the  electromagnetic  case. 

Kraus  has  approached  the  problem  from  another  point  of  view.  In 
his  dissertation  he  developed  a  "uniformized"  sphero-conal  coordinate 
system,  and  determined  a  scalar  Green’s  function  in  this  coordinate 
system.  His  scattering  body  is  a  plane  angular  sector,  which  is  a 
degenerate  case  of  one  of  the  elliptic  cone  coordinate  surfaces.  An 


u 

eigenfunction  ejcpansion  was  used  to  determine  the  Green’s  function  for 
both  the  hard  and  soft  boundary  conditions  on  the  sector.  The  angular 
functions  in  the  eigenfunction  expansion  are  called  Lame  functions  and 
are  solutions  of  two  coupled  Lame^  equations.  The  first  eigenvalue  of 
these  equations  for  both  boundary  conditions  was  determined  approxi¬ 
mately  for  several  sector  angles. 

The  article  by  Kraus  and  Levine  v;as  published  several  years  later. 
It  contains  some  refinement  of  the  problem  and  the  general  solution  for 
an  elliptic  cone  scattering  body;  no  numerical  results  were  presented. 

The  "uniformized"  sphero-conal  coordinate  system  is  used  in  this 
paper  to  solve  the  vector  problem.  This  is  one  of  the  six  coordinate 
systems  in  which  the  vector  wave  equation  is  separable.  The  coupled 
Lame^  equations  which  occur  are  solved  in  this  paper  for  the  special 
case  of  the  quarter  plane.  The  L-ame^  functions  are  a  very  general  class 
of  functions.  They  depend  on  an  ellipticity  parameter  and  two  eigen¬ 
values.  As  the  ellipticity  goes  to  one  and  the  elliptic  cones  become 
circular  cones,  the  Lame^  functions  reduce  to  Legendre  functions  and 
trigonometric  functions.  For  other  limiting  values  of  the  ellipticity 
and  the  eigenvalues,  they  become  Tchebycheff  polynomials  and  Mathieu 
functions . 

This  report  has  been  organized  in  the  same  order  as  the  problem 
was  solved.  In  order  to  determine  the  dyadic  Green’s  function,  it  is 
necessary  to  have  a  complete  set  of  vector  wave  functions,  and  in 
order  to  have  a  corrplete  set  of  vector  wave  functions  for  this  problem, 
two  conplete  sets  of  scalar  wave  functions  are  needed.  Before  the 


scalar  wave  functions  car.  be  determined,  it  is  necessary  to  understand 
the  coordinate  system. 

Consequently,  this  report  begins  with  a  discussion  of  the  sphero- 
conal  coordinate  system  in  Chapter  2.  The  peculiarities  of  the  coordi¬ 
nate  system  are  discussed  and  the  effect  of  these  peculiarities  on  the 
scalar  wave  functions  is  discussed. 

In  Chapter  3  the  scalar  wave  equation  for  this  coordinate  system 
is  presented  and  then  separated  into  the  spherical  Bessel’s  equation 
and  the  coupled  Lame^ equations .  The  boundary  conditions  are  examined 
and  the  problem  is  separated  into  four  different  boundary  value  prob¬ 
lems.  The  method  for  solving  these  problems  is  discussed  and  then  the 
first  192  eigenvalues  and  eigenfunctions  are  determined  for  the  quarter 
plane . 

The  vector  wave  functions  are  determined  in  Chapter  4,  and  in 
Chapter  5  the  dyadic  Green’s  function  is  derived.  Since  the  vector 
wave  functions  a.re  a  new  set  of  functions,  it  is  necessary  to  investi¬ 
gate  them  in  some  detalJ.  in  order  to  determine  the  dyadic  Green’s  func¬ 
tion.  This  investigation  is  chiefly  concerned  with  orthogonality  pro¬ 
perties  and  normalization.  A  suggestion  concerning  normalization  of 
/ 

the  Lame  functions  is  also  included  in  Chapter  5* 

Some  numerical  results  for  the  fields  diffracted  by  and  the  cur¬ 
rents  on  a  quarter  plane  are  presented  in  Chapter  6.  The  dominant 
behavior  of  the  fields  and  current  near  the  tip  and  along  the  edges 
is  examined.  This  is  done  for  both  an  infinitesimal  dipole  and  a 
plane  wave  source.  The  far  field  due  to  an  infinitesimal  dipole  very 


close  to  the  tip  is  also  discussed  and  some  patterns  are  presented. 

The  appendices  are  all  connected  with  Chapter  3.  Appendix  A 
includes  a  discussion  of  some  of  the  self-adjoint  properties  of  the 
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scalar  wave  equation.  Appendices  B  and  C  describe  the  calculation  of 
the  eigenvalues  and  eigenfunctions  in  detail,  and  Appendix  D  describes 
another  approach  for  calculating  the  eigenvalues  and  eigenfunctions. 
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CH/tPTER  II 


COORDINATE  SYSTEM 

This  chapter  contains  a  description  of  the  coordinate  system  and 

its  peculiarities,  and  the  effect  of  these  peculiarities  on  solutions 

of  the  wave  equation. 

Description  of  the  Coordinate  System 

The  uniforraized  sphero-conal  coordinate  system  (r, 

introduced  by  the  following  coordinate  transformation: 

[l>6]  is 

X  =  r  cos  0<>/l-k*^  cos^  cp 

(2.1a) 

y  =  r  sin  0  sin  cp 

(2.1b) 

z  =  r  cos  <p^/l-k^  cos^  0 

(2.1c) 

where 


k’2  =  i_k2,  0  <  <  1 

0  ^  0  <  n 
0  ^  9  <  2n 

r  >  0 

and  X,  y,  and  z  are  the  usual  Cartesian  coordinates. 

It  is  derived  from  the  standard  sphero-conal  system  in  the  follovf- 
ing  manner. 

Consider 

r^  =  +  y^  +  z^  (2.2a) 


T 
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(2.2b) 


=  0 


if  -  -Zf _  -  _ =  0 

b^-v^  c^-v^ 


and  the  inverse  transformation 


x  =  £e: 
be 


(2.2c) 


(2.3a) 


y  =  Vb^-v^ 

z  =  n/c^-v^ 

c  n/  c^-b^ 

where 

0  <  <  b^ 


(2.3b) 

(2.3c) 


b^  ^  ^ 

The  surfaces  r,  ii,  and  v  are  mentioned  by  Morse  and  Feshbach  [7]> 
Byerly  [8],  Moon  and  Spencer  [9],  and  several  other  authors  under  the 
heading  of  conical  coordinates. 

In  general,  the  surfaces  r,  ji,  and  v  intersect  at  eignt  points, 
which  introduces  an  eight  fold  ambiguity  in  fixing  a  point  in  space.  In 
order  to  achieve  a  one  to  one  coi’respondence  between  (x,  y,  z)  and 
(r,  p,  v),  the  variables  0  eind  cp,  and  the  parameter  k  are  used.  They 


cos  0  =  v/b 

0  <  0  <  IT 

(2.4a) 

.  u^-b^ 

sin  CD  =  -i- - 

^/c2-b2 

0  <  cp  <  2n 

(2.4b) 

are  defined  so  that 
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k  =  b/c  (2.4c) 

In  terms  of  the  new  variables,  equations  (2.3)  become  equations  (2.1). 
For  a  more  detailed  discussion  of  this  "uniformization",  see  Kraus  and 
Levine  [1,6]. 

The  surface  0  =  is  an  elliptic  cone  with  its  axis  the  x  axis 
and  its  tip  at  the  origin.  (See  Figure  2).  The  angle  between  the  x 
axis  and  the  surface  of  the  cone  in  the  z  =  0  plane  is  0^.  The  angle 
beviween  the  x  axis  and  the  line  passing  through  the  origin  and  one  of 
the  foci  of  the  ellipse  in  the  y  =  0  plane  at  the  plane  x  »  constant  is 
given  i.y  e  and  is  related  to  k  by 


1  +  tan^  e  cos^  0i 


(2.5) 


Note  that  for  0  =  0,  the  surface  is  a  plane  angular  sector  in  the 
y  =  0  plane,  centered  around  the  positive  x  axis.  For  0  =  Il/2,  it  is 
the  entire  x  =  0  plane.  For  0  =  H,  it  is  again  a  plane  angular  sector 
in  the  y  =  0  plane,  but  it  is  centered  around  the  negative  x  axis. 

When  0  =  n,  €  is  the  semi-angle  of  the  angular  sector,  and 


1  +  tan^  e 


cos^  e 


(2.6) 


Thus,  0  =  n,  k^  =  1/2  corresponds  to  the  quarter  plane  shown  in  Figure  3. 
The  quarter  plane  lies  in  the  y  =  0  plane  and  is  symmetric  around  the 
negative  x  axis.  The  corner  angle  is  II/2.  When  k^  =  0,  e  =  T[/2,  the 
corner  angle  becomes  11,  and  the  plane  angular  sector  is  a  half  plane. 
Whan  k^  =  1,  e  =  0,  and  the  plane  angular  sector  reduces  to  the  nega¬ 


tive  X  axis. 
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Note  that  for  =  1,  the  surface  0  =  0^^  is  a  right  circular  cone. 
Actually  for  this  value  of  the  coordinate  system  becomes  the  familiar 
spherical  coordinate  system. 

Now  consider  the  elliptic  half  cone  sketched  in  solid  lines  in 
Figure  4;  it  is  the  surface  (p  =  q)i.  The  surfaces  cpi  and  2Jl  -  cpi  com¬ 
pose  a  complete  elliptic  cone  with  its  axis  the  z  axis.  Its  character¬ 
istics  sire  described  in  the  same  manner  as  the  0  =  0^  surface  using  the 
parameter  k’^.  For  <p  =  0,  the  surface  is  a  plane  angular  sector  in  the 
y  =  0  plane,  centered  around  the  z  axis.  For  cp  =  ]l/2,  it  is  the  half 
plane  z  =  0,  y  >  0.  For  cp  =  H,  it  is  again  a  plane  angular  sector  in 
the  y  =  0  plane,  centered  around  the  negative  z  axis.  For  cp  =  3Il/2, 
it  is  the  half  plane  z  =  0,  y  <  0,  and  for  cp  =  2n,  it  coincides  with 
the  cp  =  0  surface. 

The  intersection  of  0  =  0i  and  cp  =  cpi  is  a  line  in  the  r  direction. 
The  intersection  of  this  line  with  the  sphere  r  =  r^  yields  the  point 
(rxj  01,  cpi)  as  shown  in  Figure  5« 

From  equations  (2.1)  the  unit  vectors  and  metric  coefficients  are 
determined  using  standard  techniques. 

r=|x+fy  +  |  z  (2.7a) 


<p  =  cos  9  cos  <p  sin  (p 
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sin  9  cos  cpyl-k*^  cos*' 


(k'=  sin*^  @  +  k*  sin*^  cp 


(2.7c) 


sin  (p^/l-k^  cos  0  -s/l-k*^  cos^ 


^/k‘^  sin*^  e  +  k*^  sin^  CD 


The  vectors  r,  9,  0  form  a  right-handed  orthogonal  system  with  each 
vector  pointing  in  the  direction  of  the  increasing  coordinate.  (See 
Figure  5.) 

The  metric  coefficients  are 


hj^  —  1 


(2.8a) 


r\/k^  sin^  0  +  k*^  sin^ 
\'^l-k^  cos^  0 


(2.8b) 


r  Vk^  sin^ 


r»2  Rin^ 


fl-k’^  cos^  cp 


(2.8c) 


Using  equations  (2.8),  the  gradient  operator  is 


fl-k^  cos^  0 


r  -Jk^  sin^  0  +  k'^  sin^  cp  30 


_ yi-k*^  cos^  cp  B 

r  -Jk^  sin^  9  +  k’^  sin^  cp  Bcp 


(2.9) 


and  the  Laplacian  is 
V®  =  i  ^  1 

r^  3r  '  3r 


r^(k^  sin^  0  +  k’^  sin^  9) 


\fl-k^  cos^  9  —  ( - 
30  \ 


+  \l  1-k’^  cos^  9  ^  (\/l-k’^  cos^  9  ~ 

39  \  39, 


'l-k^  cos^  0  A 
30 


(2.10) 


The  Effect  of  Peculiarities  of  the  Coordinate 
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System  on  Solutions  of  the  Wave  Equation 

From  the  description  of  the  cp  coordinate  surface  in  the  previous 
section,  it  is  seen  that  the  variable  cp  is  periodic  with  period  2II. 

Thus  in  the  absence  of  any  physical  boundaries  in  the  cp  direction  it  is 
necessary  that  solutions,  \Jc(r,  0,  q>),  of  the  wave  equation  satisfy  the 
following  periodicity  condition. 

\!c(r,0,cp)  =  +  2n)  (2.U) 

This  is  the  same  condition  that  occurs  in  cylindrical  and  spherical 
coordinate  systems;  it  is  a  necessary  condition  if  the  solution  is  to 
be  a  single-valued  function  of  cp. 

Now  consider  the  y  =  0  plane.  (See  Figure  6).  It  is  divided  into 
four  sectors  ^ich  axe  the  surfaces  0  =  0,  0  =  H,  cp  =  0,2!I,  and  cp  =  IT. 
The  cp  =  n  surface  is  regular  and  presents  no  difficulties.  The  cp  =  0,211 
surface  is  described  in  two  different  ways,  but  this  is  taken  care  of 
by  the  periodicity  requirement  just  discussed.  The  0  =  II  surface  is  the 
scattering  body,  and  will  have  boundary  conditions  prescribed  by  the 
nature  of  the  physical  problem.  The  0  =  0  surface  is  a  singular  coordi¬ 
nate  surface.  A  point  on  it  is  described  in  two  different  ways  depend¬ 
ing  on  whether  the  surface  is  approached  from  above  or  below.  For 
y  =  O'*",  a  point  is  described  by  (r,0,cp).  For  y  =  0“,  the  same  point  is 
given  by  (r,0,2II-cp) .  At  (r,0,cp),  0  =  y,  and  $  has  a  negative  z  compo- 

A  A  A 

nent.  At  (r,0,2n-cp),  0  =-y,  and  cp  has  a  con^ionent  in  the  positive  z 
-irection.  This  is  easily  seen  by  observing  that  the  unit  vectors 
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point  in  the  direction  cf  the  increasing  coordinate.  It  is  also  evident 
from  equations  (2.7). 

In  order  to  have  continuous  scalar  fields,  it  is  necessary  that 
\lf(r,e,cp)  and  V\}r(r,9,cp)  be  continuous.  For  continuity  of  \jf(r,0,(p)  on 
the  0  =  0  surface. 


>j;(r,0,q))  =  \|f(r, 0,211  -  cp)  (2.12) 

Using  the  periodicity  in  cp,  this  can  be  written 

\!f(r,0,(p)  =  \if(r,0,  -  cp)  (2.13) 

Inspection  of  equation  (2.9)  indicates  that  the  following  two  equations 
must  be  satisfied  for  the  gradient  to  be  continuous  at  0  =  0. 


II 

o 

-  (>-,0, 

-  cp) 

(2.14) 

80 

80 

^  (r,0,9)  = 

-  ^  (r,0, 

-  9) 

(2.15) 

8cp 

8cp 

Equation  (2.15)  is  a  necessary  consequence  of  equation  (2.13),  so  it  is 
not  an  independent  equation. 

Next  consider  the  four  lines  described  by  (0  =  n,  cp  =  0,  2Il), 

(0  =  n,  cp  =  n),  (0  =  0,  cp  =  0,  2n),  and  (0  =  0,  cp  =  n).  The  two  lines 
bordering  0  =  H  are  part  of  the  scattering  surface  and  will  have  bound¬ 
ary  conditions  prescribed  by  the  physical  problem.  Consider  the  line 
(0  =  0,  cp  =  O).  On  the  cp  =  0  surface,  cp  =  y  and  0  is  in  the  y  =  0 
plane.  On  the  0  =  0  surface  0  =  y  and  cp  is  in  the  y  =  0  plane.  In 
order  to  have  a  continuous  gradient  at  the  line  (0  =  0,  cp  =  0), 


it  is 
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necessary  that  various  coraponents  of  the  gradient  behave  properly  as 
the  line  is  approached  from  different  directions.  For  example,  con¬ 
sider  the  component  of  the  gradient  which  is  perpendicular  to  the 
y  =  0  plane.  On  the  0  =  0  surface,  it  is  given  by  Vijr  *  0.  On  the 
cp  =  0  surface  it  is  given  by  *  (p.  As  the  line  (0  =  0,  cp  =  0)  is 
approached;  it  is  necessary  that 

lim  V\|f  •  q)  =  lim  V\j/  *  0 

0  -^O'*'  cp  ^0+  (2.16) 

cp  =  0  0  =  0 

Kraus  and  Levine  have  investigated  all  of  the  necessary  conditions  on 
the  components  of  the  gradient  and  have  determined  that  they  are  auto¬ 
matically  satisfied  for  solutions  of  the  wave  equation.  For  details  on 
this  rather  ingenious  derivation,  see  their  paper  [6]. 

To  summarize,  the  "boundary  conditions"  imposed  by  the  coordinate 
system  are  given  by  equations  (2.11),  (2.13),  sind  (2.14). 

\!;(r,0,cp)  =  >!r(r,0,cp  +  2n)  (2.11) 

\lf(r,0,(p)  =  \lf(r,0,  -  cp)  (2.13) 

^  (r,0,(p)  =  -  ^  (r,0,  -  cp) 

Be  he 


(2.14) 


CHAPTER  III 


SOLUTION  OF  THE  SCALAR  HAVE  EQUATION 

This  chapter  contains  a  description  of  the  scalar  problem  and  the 
separation  of  the  scalar  wave,  equation  into  a  sphericsil  Bessel's 
equation  and  two  coupled  Lame  equations.  The  angular  boundary  values 
are  discussed  and  it  is  seen  that  the  Lame  problem  can  be  decomposed 
into  four  boundary  value  problems,  the  solutions  of  which  comprise 
two  complete  sets  of  functions.  Each  of  these  four  problems  is  dis¬ 
cussed  in  detail  and  then  the  solutions  are  tabulated. 

Description  of  the  Problem 

As  mentioned  earlier,  the  primary  purpose  of  this  report  is  to 
determine  a  dyadic  Green's  function  for  the  plane  angular  sector.  In 
order  to  do  this,  solutions  of  the  vector  wave  equation  are  needed, 
and  in  order  to  solve  the  vector  wave  equation,  solutions  of  the 
scalar  wave  equation  are  needed.  The  solution  of  the  scalar  wave 
equation  is  described  in  this  chapter. 

The  scalar  wave  equation  is 

(V2  +  k2)  i|«(r,8,(^)  =  0  (3.1) 

where  ic  is  the  usual  wave  number.  Using  the  expression  for 
(equation  2.10)  in  the  sphero-conal  coordinate  system,  this  becomes 
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Before  taking  a  closer  look  at  equation  (3.5),  the  restrictions 
imposed  by  the  coordinate  system  "boundary  conditions"  will  be 
Investigated.  From  equation  (2.11) 

+  2ti)  (3.7) 

From  equation  (2.13) 

0(0)  <!>((}.)  =  0(0)  <I>(-4>)  (3.8) 

From  equation  (2.ll») 

0'(O)  =  -0'(O)  4>(-(^)  (3.9) 


To  see  what  the  last  two  conditions  mean,  separate  the  functions 
0(0)  8uid  <{'(((1)  into  even  and  odd  parts.  Since  0  is  restricted  in 
the  range  0  _<  0  £  ir,  the  terms  even  and  odd  here  simply  mean  that 
0^(0)  =  0  and  0q(O)  =  0. 


0e(9)  +  0^(0) 

(3.10a) 

4'e((j))  +  ♦0(4') 

(3.10b) 

Using  equations  (3.10)  and  the  properties  of  even  and  odd  functions, 
equation  (3.8)  can  be  written 


0e(O)  +  0e(o)  '’’o('fr)  =  0e(O)  'J>e(<J>)  -  0e(O)  <I'o('i>)  • 


(3.11) 


This  can  be  satisfied  if 


or  if 


4’o((J>)  =  0,  i.e.  $((!))  is  even 


(3.12) 


0e(O)  =  0 


(3.13) 
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This  last  condition,  along  with  the  fact  that  0^  (O)  =  0,  implies 
that  Qq  (6)  =  0.  Thus  equation  (3.13)  implies  that 

0(e)  is  odd  ,  (3.14) 

Using  equations  (3.10)  and  the  properties  of  even  and  odd  functions, 
equation  (3.9)  can  he  written 

0q(O)  ♦g((J))  +  0q(O)  4>o((|))  =  -  0q(O)  «te((b)  +  0o(O)  «>o(4>)  .  (3.15) 


This  can  be  satisfied  if 

<I'e(<J>)  =  0,  i.e.  <!>(({>)  is  odd 


(3.16) 


or  if 

0o(O)  =  0  .  (3.17) 

This  last  condition,  along  with  the  fact  that  0o(O)  =  0,  implies 
that  0o(0)  =  0.  Thus  equation  (3.17)  implies  that 


0(0)  is  even 


(3.18) 


In  summary,  equation  (3.3)  can  be  satisfied  if  is  even  or  if 
0(9)  is  odd.  Equation  (3.9)  can  be  satisfied  if  '5>((j>)  is  odd  or  if 
0(e)  is  even.  The  only  non-trivial  combination  of  these  conditions 
is  both  0(0)  and  $((p)  even  or  both  0(0)  and  <!>((}>)  odd.  Thus  the 
coordinate  imposed  "boundary  conditions"  require  the  solution  of 
equation  (3-5)  to  be  written  in  the  form 


2h 

The  physical  boundary  conditions  are  imposed  at  the  plane  angular 
sector  and  will  be  conditions  on  0(8)  since  the  plane  angular  sector 
is  described  by  0  =  :r.  It  will  be  seen  in  the  next  chapter  that  the 
solution  to  the  vector  problem  requires  both  Dirichlet  and  Neumann 
scalar  solutions;  that  is,  two  types  of  0(6)  functions  are  needed, 
those  that  satisfy  0(ii)  =  0,  and  those  that  satisfy  0'(Tt)  =  0. 

To  completely  specify  the  problem,  it  is  necessary  to  normalise 
the  functions.  A  convenient  normalization  is  accomplished  by  setting 

0e(o)  =  «^e(0)  =  1  (3.20a) 

and 

0«(O)  =  4>'(0)  =  1  .  (3.20b) 

o  o 

Equation  (3.5)  can  be  separated  into  the  following  two  equations 

|l  -  k2  cos2  e  ^  (  j 1  -  cos2  0  ^  ) 

'  l0  '  de 

+  (v(v+l)  k^  sin^  0  +  y)  0  =  0  (3.21) 

jl  -  k'2  cos^  ij)  —  ( 1 1  -  k^  cos^  ij>  — ) 

d(Ji  d({) 

+  (v(v+l)  k'^  sin^  <{>  -  u)  4>  =  0  .  (3.22) 


Considering  the  boundary  conditions,  the  solution  of  these  equations 
is  actualiy  four  separate  problems. 


I.  Dirichlet 
A.  Even 

Solve  equations  (3. 21)  and  (3.22)  subject  to  the  boundary 


conditions 

(1)  0ei(O)  =  1 

(2)  $el(0)  =  1 


0e]U)  =  0 
^eiU)  =  ‘5’elU  + 
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B.  Odd 

Solve  equations  (3«2l)  and  (3.22)  subject  to  the  boundary 

conditions 

(1)  0ol'(O)  =  1  eol(’f)  =  0 

(2)  4>ol'(0)  =  1  *ol('}>)  =  +  27r) 

n.  Ne’jmann 

A.  Even 

Solve  equations  (3. 21)  and  (3.22)  subject  to  the  boundary 

conditions 

(1)  0e2  (0)  =  1  0e2’(’f)  =  0 

(2)  4162  (0)  =  1  'J>e2('i’)  “  ‘•’e2('J>  2tt) 

B.  Odd 

Solve  equations  (3. 21)  and  (3.22)  subject  to  the  boundary 

conditions 

(1)  0o2'(O)  =  1  0o2'(it)  =  0 

(2)  4'o2'(0)  =  1  *^02^)  =  'I’o2('^  27r) 

Before  investigating  each  of  these  problems  separately,  some 
general  results  are  discussed.  First,  note  that  each  problem  is  two 
two-parameter  Sturm-Liouville  problems.  Thus,  for  each  value  of  v  in 
equations  (3-21/  and  (3.22)  there  are  an  infinite  nuaiber  of  p's  which 
can  satisfy  each  equation.  Kraus  and  Levine  [6]  have  shown  that  only 
a  finite  nuniber  of  p '  s  are  needed  for  each  v  in  order  to  have  a 
complete  set  of  solutions.  Also  note  that  equation  (3.5)  is  a  two- 
dimensional  Sturm-Liouville  type  equation.  It  can  be  shown  that  the 
two-dimensional  Sturm-Liouville  type  operator  is  self  adjoint  and 
positive  definite  and,  therefore,  that  the  eigenvalues  v  are  all 
positive.  These  results  will  be  used  to  find  the  eigenvalues  v  and  p. 
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It  follows  from  the  self  adjoint  properties  of  the  differential 
equations  that  the  solutions  to  the  Dirichlet  and  Neumann  problems 
are  orthogonal.  The  orthogonality  relationship  takes  the  form 

Is  0n(6)  '*>nU)  0p(e)  dS  =  0  n  ^  p  (3.23) 

where  the  subscripts  n  and  p  indicate  that  the  eigenfunctions 
correspond  to  the  eigenvalue  pairs  (vn,  Pn)  ^^d  (vp,  pp) .  The  surface 
S  is  spherical.  It  is  understood  that  all  of  the  eigenfunctions  in 
equation  (3.23)  belong  to  either  the  Dirichlet  set  or  the  Neumann  set. 
Dirichlet  eigenfunctions  are  not  necessarily  orthogonal  to  Neumann 
eigenfunctions . 

The  self-adjoint  property  and  the  positive-definite  property  of 
the  two-dimensional  Sturm-Liouville  type  operator  are  proved  in 
Appendix  A.  For  proof  of  completeness  and  orthogonality,  see  Kraus 
and  Levine [1,6]. 

Method  of  Solution 

lA.  Even  Dirichlet  Problem 

Ihe  usual  approach  to  an  unfamiliar  differential  equation  is  to 
try  a  power  series  solution.  By  assuming  a  solution  to  equation  (3.2l) 
of  the  form 

00 

0el(8)  =  I  cosm0/2 

ra=0 

it  is  seen  that  there  are  two  independent  solutions,  one  with  m  even 
and  the  other  with  m  odd.  By  imposing  the  Dirichlet  boundary  condition, 
the  solution  with  m  even  can  be  eliminated.  Instead  of  writing  the 


solution  in  the  form 
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0el(6)  =  I  4i  cos“0/2 
ni—l )  3  ^  •  w  • 

it  is  found  to  be  more  convenient  to  use 

0ei(e)  =  I  Am  cos  {2m  -  1/2)  0  (3.2U) 

m 


This  set  of  equations  can  be  vritten  in  matrix  form,  and  a  determinant 
can  be  identified  which  must  be  zero  in  order  to  have  a  non-trivial 
solution.  Ince[10]  has  considered  a  similar  problem  which  must  have 
the  same  solution  as  this  problem.  He  sets  up  the  problem  in  the  same 
way  and  identifies  an  infinite  determinant  which  must  be  zero.  With 
some  rather  straightforward  manipulations,  an  infinite  determinant  of 
the  type  considei’ed  here  can  be  written  in  the  form  of  an  infinite 
continued  fraction.  The  fraction  for  this  problem  is 

h  ~  ^  *  3)  k^/9  16  (2v  -  3)(2v  +  3)  k^/22^ 

i*  1  +  k2  -  lih/9  +  1  +  k2  -  lih/25 


36  (2v  -  ;3)(2v  +  7)  k^/1229 
+  1  +  k2  -  !4h/i»9  + 


where 


p  =  h  -  v(v+l)  k^  . 


(3.26) 

(3.27) 
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equations  can  be  written 
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n  = 


-(v^-l)(v)(v4-2)  k'V8  (v2-9)(v-g)(v+l<)  k'V(l6)2 
2  -  k’2  -  nA  -  2  -  k'2  -  n/l6 

(v2-25)(vA)(v+6) 


2  -  k’2  -  n/36 
for  equations  (3.30)  and 

n  =  v(v+l)  k'2/2  +  2  -  k’2  - 


(v^-4)(v-1)(v+3)  k’^/(6)' 


2  -  k'2  -  n/9 
(v^-16)(v-3)(v+5)  k'^/(30)2 


2  -  k’2  -  n/25 
for  equation  (3.3l)  where 


(3.32) 


(3.33) 


V  =  1/2  (-n  +  v(v+l)  k»  )  . 


(3.34) 


In  order  to  find  the  eigenvalues  for  problem  lA,  it  is 
necessary  to  solve  equations  (3.26)  and  (3.32),  and  equations  (3.26) 
and  (3.33)  simultaneously.  The  solutions  are  the  eigenveilues  vv,u) 
of  the  even  Dirichlet  problem.  The  manner  in  which  this  is  done  is 
shown  in  Appendix  E.  The  eigenvalues  are  tabulated  for  the  quarter 
plane  problem  (k2=l/2)  in  Table  1  (pages  4l  to  44  ). 

Onv,e  the  eigenvalues  are  determined,  the  eigenfunctions  can  be 
found  by  solving  the  recurrence  relations.  This  again  involves  the  use 
of  continued  fractions  and  is  explained  in  Appendix  C.  The  form  of  the 
eigenfunctions  is  shown  in  Table  1,  and  the  values  of  the  coefficients 
are  given  in  Table  2  (pages  4')  to  56  ). 

IB.  Odd  Dirichlet  Problem 

Solutions  to  this  problem  may  be  expressed  in  terns  of  the 
stEuidard  Lame  polynomials.  The  eigenvalues  v  are  integers,  and  there 


are  exactly  v  p's  for  each  value  of  v.  The  solutions  can  be  found 
in  the  same  way  as  in  problem  lA.  Series  solutions  are  assumed,  and 
it  is  found  that  the  series  are  finite  for  integral  values  of  v  if 
they  6ire  written  in  the  following  forms.  For  equation  (3.21), 

N 

001^®^  =  1  ^2m+l  sin  (2m-H)  9  (3- 

m=0 
N 

0Oi(0)  =  I  A2m  sin  2m  6  (3.3( 

m=l 

I -  N 

0Ql(0)  =  jl  -  coB^  0  ^  A2111+1  sin  (an+l)  0  (3.31 

m=0 

I -  N 

0Ol(0)  =  (1  -  cos^  0  I  A21J1  sin  2m  0  (3.3( 

m=l 

where  N  is  an  integer  that  depends  on  the  eigenvalue  v.  Equations 
(3.35)  and  (3.37)  describe  one  independent  solution;  equations  (3.36) 
and  (3.38)  describe  the  other  independent  solution. 

Each  of  these  solutions  gives  rise  to  a  recurrence  relation  and 
a  continued  fraction  equation  for  the  eigenvalues.  The  recurrence 
relations  for  equation  (3.35)  are 

Ai  ((ii  -1)  .  3v(v.l)  ^  (6-v(v+l))  =  0  (3.39^ 

2  4  4 

^2m-l  ^  ((2ni-l)  (2ra)  -  v(v+l))  (3.39' 

4 

, ,  , ,2  v(v+l)  k^ 

^2m+l  ((2m+l)  (-^  -i)  +  ^  +  p) 

k2 

+  Apjn+3  -r  ((2m+3)(2m+2)  -  v(v+l))  =  0  m  >_  1 
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Each  of  these  recurrence  relations  can  be  written  in  matrix  form  and  the 
corresponding  determinantal  equations  found.  These  equations  written 
in  the  form  of  continued  fractions  are,  for  equation  (3.35)  and 
equation  (3.37) j 
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n 


-v(v+l)  k2/2  +  2  -  k2  - 


(/-U)(v-1)(v4-5)  k'7(6)^ 
2  -  k2  -  n/9 


(v^-16)(u-3)(v+5)  kV(30)^ 
2  -  k^  -  n/25 


(3.U3) 


where 

u  =  1/2  (n  -  v(v+l)  k^)  (3.Ulf) 

and  for  equation  (3.36)  and  (3.33), 


„  =  1.(2.1:2)  -  (v^-9)(M-2)(vtli)  kVcS)^ 

2  -  k^  -  n/l6 

(v^-23)(v-U)(vt-6)  kV(lt8)^ 

2  -  k^  ~  n/36  -  . 

where  p  is  again  given  by  equation  (3.^*t). 

The  same  forms  of  solution  are  employed  in  equation  (3.22) 
with  the  coefficients  B  instead  of  A. 

N 

'J’Ol^'j*)  =  1  B2m+1  (2m+l)  (j) 

m=0 

N 

•fOlU)  =  1  B2ni  sin  2m  <|) 

m=l 


'J’OlU)  =  jl  -  cos2  (j)  ^  52^1+3^  sin  (2m+l)  (|) 

m=0 


'i’OlU)  =  jl  -  k'^  cos2  .jl  ^  B2n  sin  2m  <J>  . 

m=0 


(3.i<5) 


(3.J+6) 

(3.>t7) 

(3.48) 

(3.49) 


The  recurrence  relations  are  the  same  with  k^  replaced  by  k'^  and 
p  replaced  by  -p .  The  eigenvalue  equations  are  the  same  except  that 
k^  is  replaced  by  k’^.  The  definition  of  p  is  just  the  negative  of 
equation  (3.44). 

p  =  1/2  (-n  +  v(v+l)  k'^) 


(3.50) 
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The  eigenvalues  are  found  hy  simultaneously  solving  equations  (3.UU) 
and  (3.50).  As  previously  mentioned,  the  eigenvalues  v  are  integers, 
and  there  are  v  y's  for  each  value  of  v.  For  the  actual  method  of 
computation,  see  Appendix  B.  The  eigenvalues  are  given  in  Table  1. 

With  the  eigenvalues  determined,  the  eigenfunctions  are  found 
by  using  the  recurrence  relations.  This  is  done  in  Appendix  C.  The 
form  of  the  eigenfunctions  is  shown  in  Table  1,  and  the  values  of 
the  coefficients  are  given  in  Table  2. 

IIA.  Even  Neumann  Problem 

Solutions  to  this  problem  are  also  standard  Lame  polynomiaD-S. 

The  eigenvalues  v  are  integers,  and  there  are  v+1  n’s  for  each  VEilue 
of  V.  The  even  Nevimann  solutions  can  be  found  in  the  same  way  as  the 
previous  solutions.  The  solutions  are  expressed  as  suitable  series, 
and  it  is  foxmd  that  the  series  are  finite  for  integral  values  of  v 
if  they  are  written  in  the  following  forms 

N 

0e2(e)  =  I  A2n,+i  cos  (2m+l)  0  (3.5I) 

m=0 

N 

®e2(0)  =  I  ^2m  cos  2m  0  (3.52) 

m=0 

0e2(e)  =  1  -  cos2  e  I  A2j„+i  cos  (2m+l)  0  (3.53) 

m=0 

j - N 

0e2(e)  =  ll  -  k^  cos^  0  I  k2m  cos  2m  0  (3.5>‘*) 

m=0 


where  again  M  is  an  integer  that  depends  on  the  value  of  v. 


With  a  few  exceptions,  the  recurrence  relations  are  the  same  as 
for  the  odd  Dirichlet  problem.  For  equation  {3.1?l),  the  recurrence 
relation  is 

k2 

A2m-1  T  ((2m-l)(2»")  -  «(v-:l))  (3. 

4 

+  A2in+1  ((2m+l)2  (—  _i)  +  --  +  p) 

2  2 

k2 

+  A2ni+3  —  ((2m+3)(2m+2)  -  v(v+l))  =  0 
4 

A_i  =  0,  m  ^  0 


For  equation  (3.52),  the  recurrence  relations  are 

,  /v(v+l)  k^  ,  ^  ^  k2  , 

Aq  ( - - - V  p)  +  A2  -jp  (2  -  v(v+l))  =  0 

Ao  {i,(!£ -d  .  +  „) 

2  2  2 


(3.5 


For  equation  (3.53),  the  recurrence  relation  is 


Aan-l  T"  ((2m+l)(2m)  -  v(v+l)) 

4 


(3.57) 


+  A2ni+3  ~  ((2rn+l)(2m+2)  -  v(v+l))  =  0 


A_i  =  0,  ID  ^  0 


For  equation  (3.5^),  the  recurrence  relations  are 


Ao  +  u)  +  Ap  (lyirriL^-!  )  =  0 


(3.58a) 


Ao  —  (2  -  v(v+l))  +  A2  (U(^  -1)  +  +  p) 


(3.58b) 


+  Ai,  —  (6  -  v(v+l))  =  0 


^2111-2  ~  ((2m)(2in-l)  -  v(v+l)) 

4 

.A2„((a.)2  (if 


(3.f8c) 


+  Apm+p  ~  {(2m)(2in+l)  -  v(v+l))  =  0  m  ^  2 


As  before,  each  of  these  recurrence  relations  can  be  written  in  matrix 


form  and  determinantal  equations  found.  These  equations,  written  in 
the  form  of  continued  fractions  are,  for  equations  (3*51)  and  (3.53), 


n  =  v(v+l)  k"^/2  +  2  -  - 


2  (v^  -4)(v-l)(v^3)  kV(6)2 


2  -  k?  -  n/9 


(v^  -l6)(v-3)(v+5)  k^/(3Q)^ 
2  -  k?  -  ri/25 


(3.5y) 
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and  for  equations  (3.52)  and  (S.5^)t 

^  -(v^-l)(v)(v+2)  kVs  (v^-9){v-2)(v+k)  kV(l6)^ 
^  2  -  k2  -  n/i*  -  2  -  -  n/l6 

(v2..25)(u-k)(v+6)  k^/{k8f 

2  -  -  n/36  -  ..... 


where 

V  =  1/2  (n  ~  v(vH)  k^)  . 


(3.61) 


Note  that  these  are  the  same  as  equations  (3*33)  and  (3.32)  except 
for  k^  and  k'2. 

The  same  form  of  solutions  is  assxuned  for  equation  (3.22)  with 
the  coefficients  B  instead  of  A. 

N 

‘*’e2U)  *  I  B2m+1  cos  (2m'M)  fj>  (3.62) 

m=0 


N 

‘I’e2(^)  =  I  B2in  cos  2m  (j>  (3.63) 

m=0 


^ - N 

4'e2U)  =  {  1  -  k'2  cos2  tf,  ^ 

m=0 


~  N 


<i>e2(<l>)  =  1 1  -  k'^  cos^  (j)  I 

m-O 


B2m+1  cos  (2m+l)  (t> 


B2m  cos  2m  i{) 


(3.6U) 

(3.65) 


The  recurrence  relations  are  the  same  with  k^  replaced  by  k'^  and  u 
replaced  by  -u.  The  eigenvalue  equations  are  the  same  except  that 
k^  is  replaced  by  k'^.  The  definition  of  w  is 

p  =  1/2  (-n  +  v(v+l)  k*^)  (3.66) 
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The  eigenvalues  v  and  w  are  found  in  the  same  way  and  are  tabulated 
in  Table  3  (pages  57  to  60  ) .  The  form  of  the  eigenfunctions  is 
also  shown  in  Table  3,  and  the  values  of  the  coefficients  are  given 
in  Table  k  (pages  6l  to  72). 

Together  with  the  solutions  of  the  odd  Dirichlet  problem,  the 
eigenfunctions  compose  a  complete  set  of  functions  that  correspond  to 
the  Legendre  polynomials  in  the  spherical  coordinate  system.  In  fact, 
as  the  parameter  approaches  1,  the  sphero-conal  coordinate  system 
degenerates  into  the  spherical  coordinate  system,  and  the  Lame  poly¬ 
nomials  degenerate  into  Legendre  polynomials.  For  each  value  v  =  n 
there  are  2n+l  standard  Lsme  polynomials.  In  principle,  the  solution 
to  this  problem  could  be  found  entirely  in  terms  of  these  polynomials, 
since  they  are  complete  and  can  be  used  to  represent  any  piecewise 
continuous  function.  In  actual  practice,  the  difficulties  encountered 
because  of  the  boundary  conditions  would  have  made  this  approach 
intractable.  The  Lame  polynomials  have  been  studied  by  several 
authors.  They  are  discussed  by  Ince(lO,ll],  Erdelyi(l2],  Prasad[l3], 
and  several  others,  and  are  tabulated  as  power  series  by  Arscott[l^t] . 

IIB.  Odd  Neumann  Problem 

This  problem  is  very  similar  to  the  even  Dirichlet  problem.  A 
series  solution  to  equation  (3.21)  ir  assumed,  and  it  is  found  that  the 
one  independent  solution  can  be  written  in  the  form 

002(9)  =  I  Am  sin  (2m  -  1/2)  0  (3.67) 

m 

The  recurrence  relation  is  given  by  equation  (3.25)  and  the  eigenvalue 


equation  by  equation  (3.26).  Tv/o  independent  solutions  of  equation 

(3.22)  are  found  in  the  forr. 

00 

=  1  ®2m  sin  2m  iji  (3.68) 

ni=l 

00 

^02^^)  =  I  B2ni+1  sin  (2ra+l)  (3.69) 

m=0 

The  recurrence  relations  and  eigenvalue  equations  are  the  same  as 
for  equations  (3.36)  and  (3*35)  with  replaced  by  k*^  and  u 
replaced  by  -y.  The  simultaneous  solution  of  the  eigenvalue  equations 
yields  the  eigenvalues.  The  computational  method  is  the  same  and  is 
shown  in  Appendix  B.  The  eigenvalues  are  tabulated  in  Table  3.  The 
eigenfunction  coefficients  are  determined  in  Appendix  C  and  are 
tabulated  in  Table  U.  The  odd  Neumann  problem  was  also  solved  using 
a  variational  method;  this  is  described  in  Appendix  D.  The  first 
few  eigenvalues  and  eigenfunctions  were  determined  and  compared  with 
the  solutions  obtained  by  the  more  exact  method.  The  comparison  is 
reasonably  good  and  is  shown  in  Appendix  D. 

Discussion  of  Tables 

Tables  1  through  U  show  all  of  the  eigenvalues  and  eigenfunctions 
for  V  less  than  9j  and  k-  =  1/2.  In  Table  1,  both  the  even  and  odd 
Dirichlet  eigenvalues  are  tabulated,  as  well  as  the  form  of  the 
corresponding  eigenfunctions.  Note  the  grouping  of  the  eigenvalues. 

If  n  is  an  integer  defined  such  that  n  -  1/2  v  <_  n  +  1/2,  there  are 
2n  +1  y’s  for  each  v.  Also  note  the  pattern  that  is  forming  for  the 
vedues  of  v.  As  v  increases,  the  value  of  v  corresponding  to  large 
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positive  values  of  y  is  approaching  n  +  1/2.  As  p  decreases  for  any 
one  group  of  eigenvalues,  the  value  of  v  decreases  until  for  large 
negative  y,  v  is  approaching  n.  The  same  type  of  pattern  is  noticed 
in  Table  3  where  both  the  even  and  odd  Neumann  eigenvalues  are 
tabulated.  For  the  Neumann  eigenvalues,  there  are  again  2n  +1  y’s 
for  each  v;  v  approaches  n  -  i/2  for  large  positive  y,  and  n  for 
large  negative  y.  Also  note  the  pattern  forming  for  the  eigenvalues 
of  the  odd  Dirichlet  and  even  Neumann  problems  (Lame  polynomials). 

As  y  increases,  the  eigenvalues  for  these  two  problems  are  beginning 
to  coincide. 

In  Table  U,  one  should  be  careful  when  interpreting  the 
coefficients  of  Am  sin  (2m  -  1/2)  0.  Coefficients  with  zero  and 
negative  subscripts  correspond  to  sine  functions  with  negative  angles. 
Thus  if  the  eigenfunctions  are  expanded  in  series  of  sines  with  positive 
angles,  it  is  necessary  to  reverse  the  signs  of  all  coefficients  with 
zero  and  negative  subscripts.  This  problem  does  not  arise  in  the 
^  cos  (2m  -  1/2)  0  eigenfxmctions  since  the  cosine  function  is  even. 

The  eigenfunctions  in  Tables  2  and  1*  are  normalized  such  that 
0e(O)  =  <$>5(0)  =  1  and  0^(0)  =  <I'q(0)  =  1.  Coefficients  with  magnitudes 
less  than  5  x  10“^  are  not  given.  Eigenvalues  should  be  accurate  to 
within  +  5  X  10~3  even  in  the  cases  of  least  accuracy.  This  accuracy 
is  not  limited  by  the  method  and  can  be  improved  if  desired.  Since 
some  of  the  eigenfunction  coefficients  are  very  sensitive  to 
inaccuracies  in  the  eigenvalues,  it  is  expected  that  their  accuracy 
is  not  as  good.  The  continued  fraction  method  used  to  find  the 
eigenfunction  coefficients,  explained  in  Appendix  C,  tends  to  minimize 


»40 


this  problem.  In  the  case  of  the  Lar.e  polynomials,  the  continued 
fraction  method  is  not  used,  and  it  is  expected  that  some  of  the 
coefficients  with  large  subscripts  are  not  precise.  This  occurs 
because  the  computation  of  the  large  subscripted  coefficients  usually 
involves  the  subtraction  of  two  large  numbers  in  order  to  determine 
a  small  number.  These  and  other  computational  problems  are  discussed 
in  the  appendices. 

Two  other  eigenvalues  were  computed  for  different  values  of  k^. 
The  lowest  eigenvalues  were  computed  for  =  0.1,  corresponding  to  a 
plane  angular  sector  with  a  corner  angle  of  1^3.11*°,  and  for  k^  =  0.9, 
corresponding  to  a  plane  angular  sector  with  a  corner  angle  of  36.86°. 
These  eigenvalues  are  given  in  Table  5  (page  73  ).  The  dominant 
behavior  of  the  vector  fields  near  the  tip  of  the  plane  angular  sector 
is  governed  i>y  these  eigenvalues;  this  is  discussed  in  Chapter  VI. 

All  of  the  computations  for  the  eigenvalues  and  eigenfunctions 
were  done  on  the  IBM  360/75  computer. 


Table  1  -  Eigenvalues  of  the  Dirichlet  Problem  and  the  Form  of  the  Angular  Eigenfunctions 


Ul 


H 

+ 


H 

+ 

&■  ^ 
•H 

_  G  W 
©•  ‘H 

W  H 

’2  ^  ^ 

+  m  pq 

©•^  9-  WW 

<§  W  <§  fl 
O  -H 

m  o  c  (a 

O  -ri 
O  H  to  H 

m  «  m  « 

WWWC4:^I 

I!  II  II  II  II  II 

MX!  -H  -o  H  Pi 


ni 

CO 

o 

o 

w 

M 

I 

r-> 

■> 


<D 


CD 

H 

& 

C9 

o 

o 

I 

w 


Q> 

H 

+ 


G 

•H 

Vi 


<D 

H 

+ 

i 


G 

•H 

CO 


2 

•H 

P  to 
to  H 

I  i 

«a!  < 

II  II  II  II 
H  Hi  Ph 


CD 

W 

CO 

O 

O 


CD 

|W 

CO 

o 

c; 


3 


CD 

o’ 


-P 

Pc 


bO  tsD'OXJ  hDPiXS-HbO  tjOTSXSHbO'r^X! 


Pc  PtlHPcPlPC  Plh)Pct-qPlh)p£' 


M  MOW  MOMOM  WOMOWOM 


O 

CT\ 

O 


V\  O  lA 
H  O  lA 
ON  O  J- 

d  d  d 

I 


O  O  lA  o  lA 

t-  O  I-I  o  o 

VO  LA  CVJ  lA  o 

•  •  •  •  • 

C\J  H  O  H  H 

I  I 


o  ro  lA  o  lA  m  o 
J-  i>-  roo  cvj  l>-xi- 
CO  lA  O  CO  CO  ON 
«»•»*«• 
lA  O')  H  O  O  ro  CO 

I  I 


vs 

lA  O 

O 

o  o  o  o  o 

lA  o 

ON 

(M  O 

CO 

COCfJNOxf 

O'.  Q 

CM 

J-  O 

r*l 

xr  O  (M  O  O 

xj-  o 

H  O  J-  O  H 
xl-OHOO 

Hr-iH  cvicvjcvjajcvj  cooooororororo 


tiD  Pi 


W  H 


M  O 


lA  O 
(M  <3N 
■M  H 

•  • 

CJN  l>- 


ON  O 
CJN  O 

^  o 

•  • 

-G*  -G* 


4.470  3.790 


42 


XI-HbjftXi'HM  M'n>J2HWox:rHbD‘OXJ  bOPiX:'HbOP4X!-HbOP<X!-Ht5D 


HOWOWOW  HOWOWOpqOWOW  WOWOHOWOMOWOpq 


OlTNOvnOlTN  COONOVOOOOVOir^.  ONOO  lACO  ir\  J-  OCVjlTNCVJOJ-OCOCO 

CTv  on  ON  t-- o\  o  HOD  H  ON  H  Q  t- ON  coco  ON  ocoi^NHinirn.;j-mcOr-iNococ.- 

H  OOH  ifNH  OJ  O.^  H  H  f-i  O  H  H  onj- J-  CO  C~- J-  H  00  00  00  on  H  H  t- oo 

(^IoC\joJ^-^-  J■r-(t~l^N<^]ot^lrNU^r-^H  ON  VO  H  ON.:t  CVJ  O  C\J  m  ON  ON  VO  V.’) 

rHi—i  IH  r^^-l 

I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


O  O  O  O  O  irv 
O  00  O  UN  o  o 
O  CVl  o  o  o  o 


-d-  J-  J- 


OOOOONOOOITNQO 

OOONOONOUNOHOO 

ifNO-^O-d-OHOOOO 


l^\l^N^^^lrN^rNu^u^^l^l^\lrN^rN 


OOONOirNOCNJOirNOUNOC 

OOONOVpOOOOUNOOOO 

UNOJhOJSfOCVJOOOOOO 


VOVOVOVOVOVOVDVOVDVOVOVOVO 


Table  1  -  (continued) 
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Table  2  -  Birchlet  Eigenfunction  Coefficients 


E  V: 

=0.296 

U:!0.090 

Ao  = 

1.048 

Bo 

=  1.036 

A-,  =■ 

-0.057 

B2 

=-0.034 

A_r 

0.011 

Bli 

=-0.002 

h  =■ 

-0.003 

A-2= 

0.001 

E  v=  1.425  y=0.915 

0  v=1.000  w=0.000 

E  v=l.l30  y=-0.455 

Ao  =  0.173  Bo  =  1.417 
Ai  =  0.845  B2  =-0.405 
iLi=-0.021  Bij  =-0,011 
A2  =  0,004  Bg  =-0.001 
A-2=-0.001 

Al  =  1.000  B^  =  1.000 
Am  =  0  ,m>l  Bm  =  0  ,m>l 

A  =  X.372  B^  =  1.010 

aV  =-0.299  B3  =-0.009 
A_i=-0.064  B5  =-0.001 
A2  =-0.006 

A_2=-0»C03 

E  V =2.480  y=2.670 

0  v=:2*000  p=1.500 

E  v=2.290  y=0.215 

Ao  =  0.121  Bo  =  2.470 
Al  =  0.043  Bg  =-0.015 
A_i=  0.842  B^  =  0.040 
A2  =-0.006  Bg  =  0,002 
A_2=  0.001 

Ag  =  0.500  Bi  =  1.4l4 

Aju  =  0  ,in>2  Bjjj  =  0  ,in>l 

Aq  =  0.343  B,  =  1.166 
Al  =  0.831  B3  =-0.160 
A^l=-0.113  B5  =-0.005 
A2  =-0.056 

A_2--0.001 

A3  =-0.C03 

Table  2  -  (continued) 


0  V=2.000  u=-1.500 

E  v=:2.0i40  p=-1.705 

=  I.U1I4  32  =  0.500 

Ajn  =  0,nf»l  =  0,i!i>2 

Aq  =  2.297  Bq  =  1.1*09 
Ai  =-1.000  Bg  =  0.862 
A_i=-0.325  Bj;  =-0.003 
A2  =  0.033 

A_2= -0.007 

A3  =  0.002 

A_3=-0.001 

E  v=  2.k93  U=  5.1*40 

0  v=3.000  M=3.873 

E  v=3.410  y«1.535 

Aq  =  0.023  Bq  =  4.941* 
Ai  =  0.104  B2  =-4.343 
A>i=  0.009  Bi^  =  0.391 
A2  =  0.865  Bg  =  0.007 

A.2=-0.001  B3  =  0.001 

A-l  =  0.032  B,  =  1.968 
A3  =  0.323  B3  =-0.323 
Ajj.  =  0,ib>3  Bm  =  0,itf>3 

Aq  =  0.248  Bi  =  1.539 

Ai  =  0.122  B3  =-0.552 

Kr  0.737  B5  =  0.013 

Ag  =-0.044  By  =  0.001 

A_2=-0.059 

A_3=-0.003 

0  v=3.000  y=0.000 

E  v=3.i45  y=-0.825 

0  v=3.C00  p=-3.873 

A2  =  0.707  Bg  =  0.707 

A^  =  0,m>2  Bn^  =  0,in>2 

Aq  =  0.503  Bq  =  0.276 
Aj  =  0,997  B2  =  0.820 
A_i=-0.34l  Bi,  =-0.092 
A2  =-0.172  Be  =-0,004 
A_2=  0.014 

A3  =-0.003 

A_3=  0.001 

At  =  1.968  Bi  =  0.032 
A3  =-0.323  B3  =  0.323 
=  0,in>3  Bjjj  =■  0,ni>3 

Tali3.e  2  -  (continued) 


E  v= 

=  3.010 

W= 

-3.9UO 

^  “ 

U.575 

Bl 

=  o.iui 

Ai  ■=- 

P' 

=  0.860 

/Li=- 

-1.152 

B5 

=-0.001 

A2  = 

0.283 

A_2= 

0.030 

^3  = 

0.006 

A- 3= 

0.001 

E  v=lt,li99  P=9.22.5 

0  v=4.000  ii=7.190 

E  v=4.470  p=3.790 

Ao  *  0.020  Bo=  10.54^ 
A,  =0.005  B2=-11.260 
iLi=  0.098  Blj=  1.752 
Ag  =  0.002  B,^=-  0.035 
A.2=  0.876  B3=-  0.001 

Ao  =  0.026  =  3.421 

AlJ  =  0.237  B3  =-0.669 

A^  =  0,np4  Bj„  =  0,nii>3 

Aq  =  0.077  Bl  =  2.377 

Ai  =  0.208  B3  =-1.525 

A.i=  0.04l  B5  =  0.145 

A2  =  0.757  B7  =  0.003 

A>2=-0.015 

A3  =  0.065 

Ai,  =-0.003 

0  v=4.000  u=2.190 

E  v=4.280  y=0.335 

0  v=4.000  m=-2.190 

Ai  =  0.109  Bg  =  0.806 

A3  =  0.435  =-0.153 

Am  =  0,iii>3  Bjj  =  0,m>4 

Aq  =  0.374  Bq  =  0.410 

=  0.192  B2  -  0.849 

A_i=  0.701  Bl,  =-0.263 

A2  =-0.l40  B5  =  0.004 

A_2=-0.132 

A3  =  0.007 

A_3=-0.001 

A;5  =  0.806  Bl  =  0.109 

A4  =-0.153  B3  =  0.435 
A„  =  0,in>4  Bjjj  =  0,ni>3 

E  v=l*.050  p=-2.575 

0  v=l*.000  p=-7.190 

E  v=li.005  p=-7.205 

Ao  =  0.751  Bl  =  0.313 
Ai  =  1.1*82  B3  =  0.761* 

A_,  =-0.889  Be  =-0.071* 
Ag  =-0.1*61*  By  =-0.003 
/L2=  0.103 

A3  =  O.Oll* 

A_^=  0.002 

All  =  0.001 

A,  =  3.1*21  B2  =  0.026 

A3  =-0.669  Bi|  =  0.237 

Aj^  =  0,nS>3  B^  =  0,n^l* 

Aq  =  0.001  Bo  =  0.023 

Ai  =-7.029  B2  =  0.107 
A_i=-3.5l*6  B^  =  0.871 
A2  =  1.257 

A_2=  0.289 

A3  =-0.027 

A_3=  0.005 

Ai,  =-0.001 

E  v=5.500  y=ll*.013 

0  v=5.000  p=ll.l*89 

1 

1  ..  .  .  .  . 

Aq  =  0.006  Bq=  23.021 
Ai  =  0.016  B2=-27.689 

A„i=  0.002  Bl|=  6.01*6 

Ag  =  0.093  Bg  =-0,372 
A_2=  0.001  B3  =-0.006 

A3  =  0.883 

Ajjj  =  0,in>3 

Ai  =  0.003  Bl  =  5.809 

A3  =  0.020  B3  =-1.915 
A5  =  0.188  B5  =  0.187 
A^  =  0,ta>5  Bj„  =  0,in>5 

0  v=5.000  p=5.196 

E  v=5.1*99  p=2,110 

Ag  =  0.095  B2  =  1.319 
All  =  0.306  Bl|  =-0.306 

Aji  =  0,m>l*  B„  =  0,m>li 

i 

Aq  =  0.061  Bq  =  0.658 
A,  =  0.379  Bg  =  0.992 

A  1=  0.032  Bl,  =-0.722 
Ap  =  O.70I*  Bg  =  0.071 

Al2=-0.026  Bq  =  0.001 

A3  =-0.151 

A_3=  0.002 

E  v=5.‘*90  y=7.110 


A  ,=  0.195  B5  =  0.651 
A2  =  0.011  BJ  =-0.013 
&_2*  0.T92  Bg  =-0.001 
A3  =~0.00U 
A_3=-0.069 
A_^=-0.003 


O  v=5«000  ii=0.000 


T  =  0.125  B-.  =  0.125 

3  =  O.J437  B3  =  0.1*37 

5  =-0.087  B5  =-0.087 

=  0,m>5  Bm  =  0,ra>5 
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Table  2  -  (continued) 


E  v=5.i50  u=-i.ri’o 

0  v=5.000  y=-5.196 

— 

E  v=5.015  y=-5.335 

Ao  =  0.551  Bi  =  0.U9U 

=  0.2lt5  B3  =  0.676 

A_l=  0.778  B5  =-0.172 

A2  =-0.366  By  =  0.001 

A_2=-0.264 

A3  =  0.01*6 

0*009 

All  =  0.001 

A2  =  1.319  Bg  =  9.1<72 
Ajj  =-0.306  =  0.506 

Aj„  =  0,m>U  Bj^  =  0,m>i* 

Aq  =  1.220  Bq  =  0.062 

A,  =  2.615  Bg  =  0.228 

A  ,=-2.127  Bit  =  0.785 
A2  =-1.21*7  B5  =-0.071 
A_2=  0.1*38  Bg  =-0.003 
A3  =  0.110 

A  g=-0.010 

Alt  =  0.002 

0  v=5.000  u=-11‘^89 

E  v=5.000  u=-ll.l*93 

Ai  =  5.809  Bi  =  0.003 

A3  =-1.915  B3  =  0.020 

A5  =  0.187  B5  =  0.188 

Am  =  0,m>5  =  0,m>5 

Ao=  21.  Oil*  Bi  =  0.020 
A3=-15.82l*  B3  =  0.096 
A_i=-9.110  B5  =  0.881* 
A2  =  3.952 

A_2=  1.261 

A3  =-0.261* 

A_3=~0.023 

Alt  “O.OOl* 

A_it=-0.001 

E  v=6.50':' 

ii=19.805 

0  v=6.000 

y=l6.783 

E  v=6.1*99 

y=11.455 

Aq  =  0.005 
A,  =  0.001 
A^l=  0.015 
A.2=  0.091 
A_3=  0,893 

Ajj^  =  0,!n<-3 

1*1*.  205 

B2=-57.J+27 
Bit=  15.81*0 
B6=-  1.61*1* 

Bq=  0.025 

Bio=  0.001 

Ag  =  0.002 

kl^  =  0.016 
A6  =  0.155 

Am  =  0,Tn>6 

Bi=  11.11*1 

B,=-  3.973 
B^=  0.1+38 

Bn=  0,in>5 

Aq  =  0.136 
A,  =  0.01*0 
A_i=  0.050 
A5  =  0.153 

a:2=  0.028 
A3  =  0.665 

A_3=-0.010 
Alt  =-0.059 
A5  =-0.003 

Bi  =  7.730 

Bg  =-8.81*0 

Be  =  2.259 
Bt  =-0.1’*7 

B9  =-0-002 

Table  2  -  (continued) 


0  v=6.000  w=»9«ll** 

E  y  =6.1*65  m=1«.850 

0  v=6.000  H=2.832 

Al  =  0.015  Bg  =  1.813 
A3  =  O.O7U  Bit  =-0.780 
A5  =  0,235  Bg  =  0.083 
Ajn  =  0,in>5  Bm  =  0,ib>6 

Ao  =-0.137  Bq  =  0.995 
Al  =-0.080  Bg  =  1.360 

A_i=  0.1*19  Bij  =-1.61*1 
Ag  =-0.0l*5  Bg  =  0.293 
a_2=  1-016  Be  =-0.006 
A3  =  0.032 

A_3=— 0.202 

Aij  =-0.002 

Ao  =  0.103  Bi  =  0.268 
aJ  =  0.261*  B3  =  0.652 
=-0.057  B5  =-0.162 
Ajjj  =  0,m>6  Bm  =  0,m>5 

E  v=6.282  y=0.1tU5 

0  v=6.000  »i=-2.832 

E  v=6.055  y=-3.380 

Aq  =-0.U11  Bi  =  0.688 

A,  =  0.688  B3  =  0,629 

A_i=-0.189  B5  =-0.3^3 
A2  «  1.002  Bi  -  0.026 

A_2=  2.790  B9  =  0.001 

A3  =-0.3l»7 

A_3=-0.037 

Alt  =  O.OlU 

A5  =  0.001 

Al  =  0.268  Bg  =  0.103 

A3  =  0.652  Bit  =  0.281* 

A5  =-0.162  B5  =-0.057 

Aj„  =  0,in>5  Bm  =  O.nP’^ 

Aq  =  1.056  Bq  =  0.028 

Al  =  0.396  Bg  =  0.391 

A.i=  1.256  Bit  = 

Ag  =-1.381*  Bg  =-0.163 

A-2=-0.683 

A3  =  0.303 

0.063 

Alt  =-0.008 

A_it=  0.001 

0  v=6.000  y=-9.1lU 

E  v=6.005  u=-9.160 

0  v=6.000  y=-l6.783 

A2  =  1.813  B,  =  0.015 
Alt  =-0.780  Bi  =  O.O7U 
Ag  =  O.C82  B5  =  0.235 
Aju  =  0,m>6  Bin  =  0,m>5 

Aq  =  1.781  Bi  =  0.069 
Al  =  1*.1*62  B3  =  0.199 
A„i=-3.923  Be;  =  0.808 
A2  =-2.850  =-0.072 

A„2=  1.186  bA  =-0.003 

A3  =  0.1*39 

A_3=-0.085 

Alt  =-0.009 

iLit=-0.001 

Ai=  ll.lUl  Bg  =  0.002 
A3=  -3.973  Bit  =  0.016 
A-=  0.1*38  Bg  =  0.155 
a;=  0,in>5  Bm  =  0,m>6 

Table  2  -  (continued) 


E  V 

=6.000 

p=. 

-16.783 

^0  ~ 

1+0.690 

®0 

=  0.003 

32.01+5 

B2 

=  0.015 

A_3= 

-20.297 

=  0.092 

A2  = 

10.063 

Eg 

=  0.890 

A_2= 

3.937 

h  = 

-1.139 

-0.221+ 

Al+  = 

0.018 

A-l+= 

-0.003 

A5  = 

0.001 

E  v=7.500  v-26.391  0  V=7.000  p=23.0T5  E  v=7.500  u=i6.820 


Aq  =  0.007  Bq=.  55.3**6 
Aj!  =  O.OOU  B2=-75.696 
A_1=  0.002  Bij=  24.78U 
Ag  =  O.Ollt  Bg=-  3.603 
Ao  =  0.088  B3=  0.167 
A4  -  0.886  0.002 


U.115  Aq  =  0.: 

6.380  Ai  =  0.( 

1.261+  A_i=  0.( 

Ag  =  0.( 

A_2~  0 . . 

A3  =  0.( 

A  3=  0.' 

A.  =-0.( 
A_4=-0.062 
A_^=-0.i 
A^  =  0,in:>lt 
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Table  2  -  (continued) 


0  v=7.000 

U =14.000 

E  v=7.492 

y=8.695 

0  v=7.000 

y'-6.444 

Ag  =  0.015 
Alj  =  0.059 
A^  =  0.191 

Ajjj  =  0,m>6 

Bg  =  3.314 
Bit  =-1*591 
Bg  =  0.191 
Bm  =  0,bi>6 

Aq  =-2.035 

A^  =-0.759 
A_, =-1.136 
A2  =  1.665 

A_o=-0 . 660 
A3  =  4.368 
A_3=  0.453 

Alt  =-0.865 
A_it=-0.029 

A_5=-0.001 

Bq  =  2.647 
Bg  =  3.500 

Bi,  =-6.936 
Bg  =  1.917 
Bg  =-0.127 

B3^q=-0.002 

A,  =  0.018 
A3  =  0.077 
Ae  =  0.208 
Ay  =-0.042 
A„  =  0,in>7 

B,  =  0.338 

B3  =  0.747 

B^  =-0.374 

B^  =  0.04l 
b'  =  0,n..7 

E  v=7.393 

y=2.660 

0  v=7.000 

y =0.000 

E  v=7.157 

y =-1.500 

Aq  =  0o3lt3 
Ai  =  1.527 
A.1=-0.358 
Ag  *  0.758 
A_2=-0.578 
Ao  --1.02k 
A_3=  0.199 
A]^  =  O.ll+O 
A  i.=-0.009 

A5  =  0.001 

Bi  =  0.970 
B3  =  0.644 
Be  =-0.736 
By  =  0.123 

Bp  =-0.002 

Ap  «■  0.221 

Aj.  =  0.389 
Ag  =-0.097 
Aj^  =  0,m>6 

Bg  =  0.221 

Bi^  =  0.389 
Bg  =-0.097 

%  =  0,in>6 

Aq  =-1.751 
A,  =  0.852 

A  T =-0.615 
A^  =  0.936 
A  2=  2.628 
A,  =-0.527 
A  3=-0.594 

aI  =  0.051 

A  1^=  0.017 
Ac  =  0.001 
A_5=  0.001 

B«  =-0.082 

Bp  =  0.644 

BT,  =  0.745 
Bg  =-0.329 

Bg  =  0.021 

B^q=  0.001 

0  v=7.000 

y =-6.441* 

E  v=7.020 

y =-6.660 

0  v=7.000 

y=-l4.000 

Ai  =  0.338 
Ao  =  0.7>t7 
Ac  =-0.374 

Ay  =  O.OUl 
Ajn  =  0,m>7 

B^  =  0.018 
B3  =  0.077 
B^  =  0.208 

Brr  — 

%  =  0,m>7 

Aq  =-0.501i 
A,  =-0.194 
A  =-0.547 
A'^-  2.708 
A:2=  0.442 
A3  =-0.903 

A_3=-0.072 

Al^  =  0.067 

A_it=  0.001 

A5  =  0.001 

Bi  =-0.086 
B,  =  0.378 
Be  =  0.890 
B^  =-0.183 

Ao  =  3.314 

a:  =-1.591 

Ag  =  0.191 

Ajjj  =  0,m>6 

Bg  =  0.015 

B^  =  0.059 

Bg  =  C.I9I 

=  0,m>6 
m 
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Table  2  -  (continued) 


E  v=7.000  u  =-1)1.010 

0  v=7.000  P=~23.075 

E  v=7.000  P=-23.075 

Aq  =  1.580  Bq  =  0.075 

A,  =lU.115  B.  =  O.OOU 

A.  =51.657  B,  =  0.003 

=  5.199  Bp  =  O.OU)* 

A3  =-6.380  B3  =  0.023 

aJ'=-141.656  b::  =  0.01I4 

A_i=-3.923  Bi*  =  0.175 

Ac  =  I.26U  B^  =  0.12)4 

A  ,=-28.1427  BiJ  =  0.090 

Ag  =-U.075  Eg  =  0.778 
A^2-'  1*556  Bg  =-0.069 

A3  =  0.905  Bj^q=-0.003 
A_3=-0.192 

Al^  =-0.05h 

A_ij=  0.003 

A5  =-0.001 

Ay  =-0.0)42  B^  =  0,ni>7 
Ajj  =  0,m>7 

Ag  =15.)H9  =  0.89)4 

A_2=  8.990 

A3  =-2.'kk6 

A_3=-0.670 

Ai*  =  0.121 

A_l;=  0.010 

A^  =  0.002 

E  v=8.500  y=3l4.389 

0  v=8.000  u=30.367 

E  v=8.500  y=23.191 

A_i=  0.002  Bq  =)4)4.092 

Ajj  =  0.002  Bi  =33.l)t2 

Aq  =  0.391  Bi  =15.87)4 

A_2=  0.013  B2=-62.390 

Ag  =  0.011  B3=-15.59)* 

Ai  =  0.101  B3=-22.791 

A_^=  0.087  Bi,  =23.282 

Ag  =  0.116  B5  =  3.291 

A_i=  0.109  Be  =9-335 

A_l*=  0.898  Bg  =-)4.316 
Ajjj  =0jni<— )4  Bq  =  0.336 
Bio=”0. 00)4 

Ajn  =  0,in>8  Bj  =-0.200 
Bjn  =  0,m>7 

A2  =  0,030  B7  =-1.1492 

A„2=  0.039  Bn  0.072 
A3  =  0.13)4  B,,=  0.001 

a:3=  0.022 

Al,  =  0.632 

A_i;=-0.008 

A5  =-0.056 

Ag  =—0.002 

Ajn  =0,in<-)4 

Tab3.e  ?  -  (continued) 


0  v=8.000  u=19.876 

E  v=8.U97  u=13.590 

0  v=8.000  ..'■•=10.9l»9 

Ai  =  0.002  Bg  =  li.6l4 

Aq  =-5.59‘‘  Bq  =-2.223 

Aj  =  0.017  3,  =  0.705 

A3  =  0.012  =-2.999 

A^  =-1.931  Bp  =-3.223 

aJ  =  0.059  B  =  1.233 
Ag  =  0.165  Be  =-0.726 

A5  =:  O.OI19  Bg  =  0.678 

A_3^=-2.01*1*  B),  =  9.701 

Ay  '  '2  Bg  =-0.037 

Ap  =-1.0l*2  Bg  =-3.670 

Ag  =-0.032  By  =  0.092 

A„  =  0,in>7  Bj„  =  0,m>8 

A_j=  3.631*  Bg  =  0.1*28 
A-j  =-0.616  B,q=-0.007 
a:.3=10.201* 

Ai,  =  0.1*11 

A^!,=-1.995 

A5  =-0.026 

Ag  =— 0 . OOl 

=  0,m>8  Ej„  =  0,ra>7 

E  v=8.1;63  u=5.885 

0  v=8.000  y=3.1|l*l 

E  v=8.275  y=0.535 

Aq  =  0.378  B,  =  1.957 

A,  =  0.0l*8  Bp  =  0.251 

A^  =  2.827  B.=-26.1*90 
=  0.573  Bp  =23.^183 

A^  =  0.170  B3  =  1.051 

A3  =  0.163  bJ  =  0.386 

A_i=  O.76U  Be  =-2.712 
Ag  =-0.127  By  =  0.753 

A^  =  0.267  B.  =-0.206 

A  ,=-1.267  B.  =17.981* 

Ay  =-0.065  Bg  =  0.023 

Ap  =  0.221  3g=-l6.021 

A_2=  0.1*00  Bn  =-0.0U9 
A3  =-0,22l*  B^,  =-0.001 

A_3=— 0 . 503 

Al^  =  0.076 

A_i,=  0.070 

A^  =-0.00U 

Ajjj  =  0,in>7  Bjjj  =  0,m>8 

A  p=-l.l*79  Bg  =  2.370 
A3  =-0.776  B  =-0.025 

A:3=  0.815  Bjp=-0.001 
A^  =  0.177 

AJ4=-0.083 

Ae  =-0.006 
a;.5=-o.ooi 

0  v=8.000  y=-3.1*l*l 

E  v=8.073  y=-l*.170 

0  v=8.000  y=-10.9l*9 

A2  =  0.251  =  0.0>*8 

Aq  =  0.2l*5  B;^  =-0.660 

Aj^  =  0.705  Bg  =  0.017 

A\^  =  0.386  B3  =  0.163 

A  =-0.005  B  =  O.85I* 

A  =  1.233  B,  =  0.059 

Ag  =-0.206  Be  =  0.267 

a\=  0.071  B?  =  1.25l» 

a;^  =-0.726  Bg  =  0.165 

AI  =  0.092  Bg  =-0.032 

Ag  =  0.023  By  =-0.065 

Ap  =-0.00.’*  By  =-0.1*76 

Ajn  =  0,m>8  B„  =  0,rii>7 

A_2=  l.Oll*  Bp  =  0.027 
Ao  =  O.OOl*  Bf.=  0.001 
A_3=-0.352 

Aj,  =-0.001 

A_i,=  0.027 

a;  =  0,m>7  Bnj  =  0,m>8 
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Table  2  -  (continued) 


E  v=8.005  u=-11.01 


0.136 

0.063 

0.lll6 

1.306 

-0.163 

-0.601 

0.039 

0.078 

-0.002 

-0.001 


Bq  =-0.312 


Bg  =  1.312 


0  v=8.000  u =-19. 876 

1 - - -  1 

E  v=8.000  *'=-19.879 

Ag  =  l*.6ll*  B,  =  0,002 
Ai,  =-2.999  B3  =  0.012 
Ag  =  0.678  Be  =  0.0!*9 
Aq  ---0.037  Bj  =  0.162 
Aju  =  0,m>8  B„,  =  0,nf7 

Af.  =  0.81*7  B,  =  0.127 
a!:  =  1*.1*22  B^  =  0.037 

A^^ =-2.183  =  0.161 

a"  =-l*.015  B^  =  0.71*1* 

A  9=  1.061  Bg  =-0.066 
A3'=  1.163  bJ, =-0.003 
A_3=-0.187 

Aj,  =-0.120 

A_i,=  0.010 

A^  =  0.002 

A,  =33.11*2 
A:>=-15.59l* 

Ac  =  3.291 
Aj  =-0.200 

Ajj  =  0,m>7 


0.002 

0.011 

0.116 

0  .in>  8 


A  =1*2.222 
a;=-31*.393 
A-^=-2l4.936 

a;^=11*.399 
A  2=  7.316 
A3  =-2.908 

A_3=-0.982 
Ai,  =  0.21*1 
A  1,=  0.0l*3 
Ac  =-0.003 
A_5=  0.001 


0.002 

0.013 

0.088 

0.896 


0  v=9.000  p=38.660 

0  v=9.000  u=26.751 

0  v=9.000  y =16. 1*13 

Ae  =  O.OOl  B  =-11.517 

Aj  =  0.010  B3  =  6.208 
Aq  =  0.103  B^  =-1.776 
a1  -  0,m>9  b2  =  0.2l*7 
Bq  =  0.116 
B;Ji  =  0,in>9 

A  =  0.002  B  =10.015 
A?  =  0.010  B?  =-6.916 
Ag  =  O.OHl  Bg  =  1.685 
Ag  =  0.11*0  Bg  =-0.133 
Ajn  =  0,m>8  Bj^_  =  0,m>8 

A  =  0.003  B  =  1.021 
A^  =  O.Oll*  B3  =  1.612 
A^  =  0.0l*7  B^  =-1.1*05 
Aj  =  0.137  B'  =  0.31*6 
Ag  =-0.027  Bg  =-0.028 
a;  =  0.m>9  B^  =  0,in>9 
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Table  3  -  Eigenvalues  of  the  Neumann  Problem  and  the  Form  :>f  the  Angular-  Eigenfunctions 


H  -o  x:  to  -o  Oi  -H  ta  H  -o  jc  -H  H  -o  x: 


MOW  WOMOM  MOMOMOW 


0 

0x1-0 

CM  lA  0  (M  <M 

0 

0  ONO 

CO  ON  0  lA  CO 

0 

lAH  lA 

t-  C-  0  lA  t- 

0* 

0*  0’ d 

1  1 

H  d  d  H  H 

1  1 

O  lA  O  O  O  O  O 
lA  c\j  lA  lA  lA  0\  lA 
a\\D  CO  ON  CO  ON 


CO  CVJ  O  O  O  CO  CO 

.  I  I  I 


0 

0  X/  0 

0  lA  0 

lA  0 

0  0  0  CO  0  0  0 

0 

0 

OHO 

0  ONO 

lA  0 

0  CM  0  0  0  CJN  0 

0 

0 

oco  0 

0  lA  0 

C3NO 

0  lAOCO  0  CTNO 

0 

d 

r— i  0 

CM  H  CM 

H  CM 

CO  CM  CO  CM  CO  CM  CO 

X}- 

3.505 


Table  3  -  (continued) 


e 


CD 


>> 

+3 

•rl 
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ft  -H  bO  ’O  P)  -H  bO  H  T-J  Xi  -H  r-l  -o  ^  -rl  ft  TO  XJ  bO  ‘O  ft  -rl  bO  TO  ft  ft  bO  TO  ft  ft  bO 


HOWOMOW  HOHOHOWOWOP-q  WOHOWOWOMOWOM 


VO  O  IfVVO  vCN  ft 

xj-  VO  ft  xt-  ON  ft 

VO  CVI  COVO  ft  CVi 

•  •  •  •  •  • 

CVJ  H  O  CVI  cy  t>-  t- 

I  I  I  I 


xtoo  coo  CJMTVOVO  COftxJ- 
(JV  ft  VO  00  VO  C7VVO  V^OVO  OV  o\ 
x^  cy  covo  coxf  CO  CVI  coxi-  xt 

cH  On  ir\  CO  o  rH  ir\  ir\  r— l  r-l 

rH  r— 1  rH 

I  I  I  I  I  I 


-d■CV]U^O^-CVJOlf^t^-OlA  CO-d- 
COHvOVpOONOHOCVJVOCOOO 
O  H  o  iTvVO  O  O  ir\  H  H  t>- 
**•••••»*•«*• 
VD-^  CJVD-COft  o  cococrvovvovo 

ft  >-<  ft  ft 

I  I  I  I  I  I 


o  c~o  O  OOO  o 
O  ft  Oxf  O  CJVO 
OVO  O  CDO  CJVO 


OOOOOCrvOxJ-OOO 

OOOCOOChOCOOOO 

OlCNOlCNOt-OOVOOO 


oooc~-ot>ot>-ot>ooo 

OOOOOCV!OCV)OCJVOOO 

OicvOirvOVOOOVOCTvOOO 


xf  coxt-  cox^  coxf  lAxt  lc^x^  irvxj-  irvx}-  umcmcv  vo  mwo  icvvo  iavo  iavo  irvvo  vo  vo 


Parity _ 0(0) 


H  -O  X3  -H  H  'i-a  -H  H  ^  -H  i-«  43  bO  "O  ft  -H  bO  'O  ft  -H  bO  -ra  ft  .H  bO  Tj  ft  -H  bO 


WOWOWOWOHOWOWOM  WOWOWOWOWOHOWOMOH 


vi)-3--3-oooir\oooooocviJ-  irwo 
t-oHJ'O  t-t>-fni>-oooHt—t>- 
O  CO  O  j-  t>-vO  C^VO  t>-  iTv  o  o  o  o 


CO  O  VO  'X>  lA-d-  lA  O  ir\  J-  lAVO  O  CO  CO 

VO  c7\co  H  m  CM  CO  CT\  o  CO  rovo  ro  t-co  VO  VO 
ro  lAco  coovocvTOOVorocJvOcocoroco 


CO  OV  J-  HVOJ-H  OHVOVO^J-COrO 

cvjhhh  HHcvjcvj 

I  I  I  I  I  I  I  I 


O  VO  (TWO  H  CO  ^  CM  O  CO  O 
CO  CM  r-i  H  r-'  I— I 

I  I  I 


H  CJV  (7\  O  O 
ft  rM  rH  CO  CO 
I  I  I  I  I 


OOOftOOOCMOOOOOOO 

OOOOOft-OCjNOCOQOOOO 

OlAOlAOlAOC^OCTvOOOOO 


OQOOOOOCMOCOOirvOOOOO 

OOOOOftOcOO^OCTvOOOOO 

OlAOlAOlAOVOOONOCTvOOOOO 


C^vo  i>-VO  t>-VO  t>VO  c^vo  C-  l>-  1>-  CO  C^CO  l>-CO  C^CO  t-CO  [>-03  L-^CO  CO  CO  CO  CO 


Table 


Table  h  -  Neumann  Eigenfunction  Coefficients 


E  v=0.000  ji=0.000 


Aq  =  1.000  B  =  1.000 

Ajj,  =  0,m>0  =  0,m>0 


E  v=1.000  y=0.500 

0  v=0.8lU  u=-0,19i* 

E  v=l.C00  u=-0.500 

Ai  =  1.000  Bq  =  l.klk 
Ajj  =  0,m>l  Bjjj  =  0,m>0 

r— . —,1—  , 

Aq  =-1.1j73  B,  =  0.961r 
Ai  =  0.21U  =  0.010 

fCi=  0.030  Be  =  0.001 
A2  =  0.007 

A.2=  0.002 

A3  »  0.001 

Aq  =  l.klk  B^  =  1.000 
Ajjj  =  0,m>0  Bjj^  =  0,m>l 

E  v=2.000  y=1.732 

0  v=i.595  y =0.795 

E  v=2.000  y =0.000 

Aq  =  0.131*  Bq  =  1.866 
A2  =  0.866  B2  =-0.866 
Aj„  =  0,m>2  Bm  =  0,m>2 

Aq  =  O.17U  Bi  =  1.160 
A^  =  0.689  B3  =-0.0148 
A_3^=-0.026  Be  =-0.003 
A2  =-0.005 

A_2=-0.001 

A^  =  1.41I4  Bi  -  1.414 
Ajj  =  0,m>l  B^  =  0,m>l 

0  v=1.955  y=-l.552 

E  v=2.000  y =-1.732 

Aq  =-1.360  Bo  =  0.496 
Ai  =  0.568  bJ  =  0.002 
A_i=  0.179 

A2  =-0.016 

A_2=  0.004 

A3  =-0.001 

Aq  =  1.866  Bq  =  0.134 

Ap  =-0.866  Bp  =  0.866 

Ajj;  =  0,m>2  Bpj  =  0,in>2 

t 


Table  -  (continued) 


E  v=3.000 

y=3.950 

Ai  =  0.138 

Bg  =  3.856 

A3  =  0.862 

Bp  =-2.442 

A,2  =  0,in>2 

=  0.m>2 

0  v=2.803 

j 

y =-0.350 

Ag  =  0.346 

Bg  =  0.583 

Ai  =  0.718 

B^  =-0.038 

A_i=-0.180 

Br  =-0.002 

Ag  =-0.091 

A.2=  0.003 

A3  =—0.003 

0  v=2.520  U=2-625  E  v=3.000  u=0.950 


3.856  Aq  =-0.060  B,  =  1.601  Ag  =  0.389  B,  =  1.364 

2.442  A,  =-0.022  Bt  =-0.193  Ap  =  1.025  B.,  =-0.364 

0,ai>2  A_i  =-0.397  Be  =-0.004  A“  =  0,in>2  B'^  =  0,in>3 

Ag  =  0.033  ^  “ 

A  2=  0*001 


E  v=3.000  y =-0.950  0  v=2.990  y =-3.890 


=  1.364 
=-0.364 
=  0,ih>3 


=  0.389 
=  1.025 

=  0,m>2 


E  v=3.000  y=-3.950 


Ag  =  3.856  B,  =  0.138 
Ag  =-2.442  B3  =  0.862 
Ajjj  =  0,ni>2  B^  =  0»ni>3 


E  v=4.000  y=7.211 

0  v=3.505  y=5.430 

Ag  =  0.019  Bg  -  7.246 
Ap  =  0.106  Bp  =-7.120 
Aii  =  0.875  Bi,  =  0.874 
Am  =  0,m>4  B^  =  0,m>4 

Ag  =  0.013  Bi  =  2.505 

■-  0.035  B3  =-0.524 
k^2r  0-0^5  Be  =  0.012 
Ap  =  0.275  Bj  =  0.001 
Al2=-0.001 

E  v=4.000  y=2.646 


.479 

B;l  =  2.350 

.935 

B  =-0.936 

,in>3 

BjJ  =  0,in>3 
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Table  U  -  (continued) 


0  v=3.617  p=1.260 


v=l».000  u=0.000 


[O  v=3.9l*0  u=-2.315 


-0.106 

-0.002 


Aq  =-0.173 
=-0.088 
A_i=-0.1»55 
k2=  0.038 
A_2=  O.OU7 
A_3=  0.002 


v=l4.000  w=-2.6U6 


0.375 

0.833 

-0.208 

0,Bl>l4 


0.375 

0.833 

-0.208 

0,m>U 


Aq  =  0.365 
k  -  O.71J+ 
A^,=-0.U09 
k  =-0.211 
Af2=  O.OUU 

A3  =  0.005 

A_3=  0.001 


0.071 

0.363 

-0.030 

-0.001 


v=3.998  ii=-7.195  E  v=U.000  ^=-7.211 


A,  =  2.350 

A3  =-0.936 
Ajj  =  0,in>3 


O.U79 

0.935 

0,in>3 


A-  =-2.639 
a!:  =  1.836 

A  ^=  0.927 

Ag  =-0.327 

A  2=-C.075 
A3  =  O.GO7 

A_3=-0.001 


0.026 

0.237 


A  =  7.2l<6 

a;  =-7.120 


=  0.019 
=  0.106 


Aj^  =  O.87U  Bf  =  0.875 

Ajjj  =  0,Tn>U  =  0,in>l| 


E  v=5.000  g=li.U9l' 

0  v=5*.500  u=9.2l8 

E  v=5.000  u=5.363 

k^  =  0.020  Bq  =16.709 

A3  =  0.096  Bp=-17.796 
A5  =  0.885*  Bij  =  2.501 
A^  =  0,in>5  B^  =  0,m>l) 

Aq  =-0.003  =  U.28I 

Aj  =-0.001  B3  =-1.250 
A_2^=-0.023  B5  =  0.092 
A_2=-0.210  B^  =  0.002 

Ajji  =0  ,in^  -2 

Aq  =  0.106  B,  =  3.133 
A2  =  0.363  B3  =-2.1»68 
kl^  =  O.9U6  B3  =  0.335 
Ajjj  =  0,m>l»  B^  =  0,in>5 

Table  4  -  (continued) 


0  v=4.530  y=3.680 

E  v=5.000  y=1.369 

0  v=4.795  y=-0.495 

Aq  =  0.028  Bg  =  0.969 
A]_  -  0.086  =-0.245 

A_i=  0.015  Bg  =  0.006 
Ag  =  0.303 

A_g=-0 . 006 

Ao  =-0.028 

Aij  =-0.001 

A,  =  0.471  Bq  =  C.704 

A3  =  0.685  Bg  =  1.151 

Ac  =-0.156  B.  =-0.442 

A^  =  0,in>5  bJ  =  0,m>4 

Aq  =-0.295  B.  =  0.112 

A,  =-0.l4l  =  0.420 

A  ^ =-0.459  B^  =-0.073 

a"  =  0.157  B^  =-0.001 

A  g=  0.125 

A3  =-0.015 

A_3=-0.002 

E  v=5.000  y =-1.369 

0  u=4.984  y=-5.230 

E  v=5.000  y=-5.363 

Aq  =  0.704  B,  =  0.471 

Ag  =  1.151  Bo  =  0,685 
=-0.442  Be  =-0.156 
Ajj  =  0,in>4  B^  =  0,m>5 

Ao  =  0.406  B«  =  0.059 
aI:  =  0,867  Bf  =  0.256 

A  ,=-0.702  B^  =-0.022 
aZ  =-0.409  Bo  =-0.001 
Ala*  0.143  ^ 

A3  =  0.035 

A_3=-0.003 

A4  =  O.OOl 

A,  =  3.133  B  =  0.106 
a::  =-2.468  Bp  =  0.363 
A^  =  0.335  bJ  =  0.946 
A^  =  0,m>5  Bj^  =  0,in>4 

0  v=5.000  y=-ll,491 

E  v=5.000  y=-11.494 

Aq  =-4.535  B,  =  0.003 
Al  =  3.417  B3  =  0.020 

A  ,=  1.967  Be  =  0.188 

Ag  =-0.854  ^ 

A_g=-0.272 

A3  =  0.057 

A_3=  0.005 

Aj^  =  0.001 

Aq  =16.709  B,  =  0.020 
Ap=-17.796  B-  =  0.096 
A^  =  2.501  B^  =  0.884 
Aj^  =  0,m>4  bJ  =  0,in>5 

Table  U  -  (continued) 


E  v=6.000 

11=16.784 

Aq  =  0.003 
Ag  =  0.015 
Aj^  =  0.092 
A5  =  0.890 

\  =  0,m>6 

b.  =35.528 
B^=-44.652 

Bj^  =11.067 
Bg  =-0.943 

B„  =  0,m>6 

El 

0  v=5.507 

y=7.060 

Aq  =— 0 . 0l4 
A,  *-0.005 
A  ,=-0.055 
Ap  =-0.003 
A:2*-0.223 
A3  =  0.001 
A_3=  0.020 

A_J=  0.001 

Bg  =  1.439 

Bj.  *-0.534 
Bg  *  0.042 
Bq  s  0.001 

E  v=6.000 

11=0.000 

A^  =  0.884 
A3  =  0.795 

Ac  =-0.265 

a;  =  0,ni>5 

B,  =  0.884 
B3  *  0.795 
Be  =-0.265 

B„  =  0,m>5 
m  ’ 

0  v=5.500  u=lU,012 


0.003 

B 

0.018 

B 

0.170 

B 

0,m>3 

B 

E  v=6.000  y=3.507 


Aq  =  0.124 
Ap  =  0.347 
Af  =  0.672 
Ag  =-0.142 
A,,  a  0,in>6 


Bp  =  1.104 

at  =-1.034 
Bg  =  0.142 
B„  =  C,in>6 


0  v=5.927  W 


A  =-0.301 

A?  =-0.115 
A  ,=-0.365 
Ap  =  0.352 

a:2=  0.187 

A3  =-0.072 
A_3=-0.016 
A^  =  0,001 


Bp  =  0.099 
bJ  =  0.282 
Bg  =-0.054 


E  v=6.000 

u=9.i65 

A  =  0.118 
A^  =  0.324 
A^  =  0.972 
A^  =  0,m>5 

3  =  6.600 
B^  =-6.158 
B^  *  0.972 
bJ  =  0,in>5 

0  v=5.627 

y =1.692 

Aq  *-0.032 

A,  *  0.172 

A  ,*-0.017 

Ap  =  0.330 

Af2=  0.015 
A3  =-0.076 

A_3=-0.001 

B,  =  0.171 
B.  =  0.503 
B^  *-0.143 
B'  =  0.005 

E  v=6.000 

y=-3.507 

Aq  =  O.7G9 
Ag  =  1.104 
Ar  =-1.034 
Ag  =  0.142 

An,  =  0,m>6 

Bq  =  0.124 

Bp  =  0.347 

B.  =  0.672 
Bg  =-0.142 

Bja  =  0,m>6 

t 


Table  Ij  -  (continued) 


0  v=5.997  u=-9.120  E  v=6.000  }j=-9-165  0  v=6.000  ^=-16.783 


Aq  =  0.552 
Ai  =  1.376 
A_i=-1.217 
Ag  =-0.878 
A„2=  0.367 
A3  =  0.135 
A_3=-0.026 
=-0.003 


0.008 

0.0li5 

0.197 

-0.017 

-0.001 


A,  =  6.600 
a::  =-6.158 
A^  =  0.972 
^  =  0,in>5 


Table  -  (continued) 


0  v=6.501 

E  v=7.000  g=6.708 

0  v=6.5l*0  g=l*.675 

Aq  =  O.sltl  B2  =  2.299 

Ai  =  0.010  Br  =-1.125 
A_2=  0.013  Bg  =  0.151* 
Ag  =  O.OU3  3«  =-0.003 
A_2=  0.007 

A3  =  0.187 

A_3=-0.003 

Aj,  =-0.017 

A5  =-0.001 

Ajl  =  0.ll*3  Bq  =  1.693 

A-  =  0.296  B  =  1.763 

Ac  =  O.70I1  Bl  =-2. 1*1*6 

tZ  =~0.ll*3  B.  =  0.1*03 

aJ  =  0,in>7  B°  =  0,in>6 

Aq  =  0.037  B,  =  0.266 
aI:  =  0.023  b;:  =  o.6i*6 
A-*-  =-0.100  B^  =-0.271* 

A“-‘-=  0.013  B'  =  0.023 

A  2=-0.23l*  ' 

A3  =-0.009 

A_3=  0.01*9 

Aij  =  0.001 

E  v=7.000  u=1.770 

0  v=6.792  y=-0.630 

E  v=7.000  g =-1.770 

Aq  =  0.283  B,  =  0.852 

Ap  =  0.638  B-  =  0.6li 

Ai,  =  0.692  Be  «-0.53lt 

Ag  =-0.200  Br,  =  0.071 

Ajq  =  0,in>6  Bj^  =  0,in>7 

Aq  =-1.1*92  Bp  =  0.ll*8 

A  =  1.21*4  Bl  =  0.315 

A-^,  =-0.589  Bg  =-0.099 

a“  =  1.520  Bo  =  O.OOl* 
A:2^  1.1*86  ^ 

A3  =-0.703 

A_3=-0.27l* 

Ai,  =  0.051 

A_^=  0.003 

A~  =  0.001 

A  =  0.852  B  =  0.283 

A^  =  0.611  B°  =  0.638 

A^  --O.53I*  Bf  =  0.692 

A^  =  0.071  B?  =-0.200 

Ajjj  =  0,m>7  =  0,iii>6 

0  v=6.980  g=-6. 500 

T’.  v=7.000  g =-6.708 

0  v=T.000  g=-ll*.00i! 

Aq  =-0.113  B  =-0.056 
A;  =-o.oU3  b::  =  0.082 
a;;;j^=-o.i22  Bc  =  0.223 
A2  =  0.556  B^  =-0.0l*l< 
A_2=  0.097 

A3  =-0.183 

A  =-0.016 

Ai[  =  0.013 

A  =  1.693  B  =  0.11*3 

Ap  =  1.763  B^  =  0.296 

Af  =-2. 1*1*6  B^  =  O.70I* 

Ag  =  0.1*03  B^  =-0.l!f3 

A„  =  0,m>6  Bjg  =  0,n’.>7 

A  =  1.176  B  =  0.008 
a!:  =  3.852  Br  =  0.036 
A-^  =-2.930  Bg  =  0.160 

Ap  =-3.022  Bg  =-0.0ll* 

A  =  1.163  B°  =-0.001 
A'3  =  0.672 

A  3=-0.11*3 

Ar  =-0.0l*O 

A_j^=  0.002 

A5  =-0.001 

Jk  ' 


Table  -  (continued) 


E  v=7.000  p=-lU.0llt 


v=7.000  u =-23. 075 


^1 

=10.906 

0.022 

A  = 

-17.244 

0.004 

=84.943 

B 

A3= 

-13.548 

^2  = 

0.089 

=13.906 

4  = 

0.023 

=-110.135 

=  4.005 

B^  = 

0.310 

=  9.490 

b5  = 

0.124 

=29.321 

4 

=-0.363 

Bg  = 

0.993 

Ap^ 

=-5. 147 

7 

4 

=-2.715 

4 

=  0,ni>7 

m 

0,m>6 

& 

=-2.333 
=  0.816 

A^ 

m 

=  0,m>6 

B^ 

m 

v=7.000  y =-23. 076 


0.003 

0.014 

0.090 

0.894 

0,in>7 


A_3=  0.224 
=-0.040 

A  .=-0.003 

A^  =-0.001 


v=8.000  p=30.368  0  v=7.500  y=26.597  E  v=8.000  w=19.880 


Ap  =  0.002 

=  0.013 

Ag  =  0.088 

Ao  =  0.896 

Ajjj  =  0,m>8 


Bq  =89.332 

B2=-126.26 

2^  =44.337 


=  0.001 


-126. 26S  A"  =  0.001 

=44.337  Ap  =  0.002 

Bg  =-7.370  A^  =  0.012 
B«  =  0.970  aJ  =  0.124 
"  =  0,ra>8  -  0,m>4 


B  =28.305 

B^=-13..473 
=  2.988 
=-0.266 
B^  =  0.004 


0.025 

0.080 

0.303 

1.006 

0,ni>7 


B  =24.157 
B-‘-=-32.109 
b3  =10.373 
b5  =-1.008 

b”^  =  0,in>7 


0  v=7.500  'J=l6.8l6  E  v=8.000  y=11.036 


v=7.510  u=8.625 


K  =-o- 

A°  =-0. 

a\=-o. 

Ag  =-0. 

A_p=-0. 
Ao  =-0. 
Ai,  =  0. 
A_i^=  0, 
A  5=  0. 


3.889 

■2.329 

0.458 

-0.025 


A  =  0.030  B 
A°  =  O.iOl  B° 
Af  =  0.276  Bf 
Ag  =  0.739  bJ 
Ao  =-0.l46  Bo 
Am  =  B 


2.464 

1.961 

-4.877 

1.597 

-0.145 

0,m>8 


A  =-0,084 
aJ  =-0.031 
A  =-0.047 

Ap  =  0.063 

A  p=-0.028 
A“  =  0.164 

a: 3=  0.019 

Ajj  =-0.033 
A_i*=-0.001 


0.426 

0.910 

-0.537 

0.077 

-0.001 
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Table  li  -  (continued) 


Table  4  -  (continued) 


0  v=8.000  y =-30.368 

E  v=8.000  y =-30. 368 

Aq=-38.665  Bi^  =  0.002 
A^  =31.492  Bg  =  0.011 

A  ,=22.833  Efl  =  0.116 
a;=-13.184  ° 

A_2=-6.699 

A  =  2.663 

A  3=  0.899 

Aj"  =-0.221 

A  ^=-0.039 

A^  =  0.003 

A,  =89.332  B  =  0.002 
A^=-126.269  Br  =  0.013 
Af  =44.337  B;  =  0.088 
Ag  =-7.370  B«  =  0.896 
Aft  =  0.970  B^  =  0,nt>8 
aJj,  =  0,m>8 

E  v«9.000  y=38.660  0  v=8.500  y=3U.389  E  vs9.000  y=26.752 


Ao  =  0.002 

Ac  =  0.013 
AX  =  0.087 
Ai  =  0.898 
A^  =  0,in>9 


Bo=833.672 

B2=-li97.523 

B,=U38.»i50 

b2=-76.676 

Bo  =  3.1*90 

B„  =  0,m>8 

m 


A  2=-0.002 
A  >-0.011 
k~?=-0.109 
Ajj^=0,m<-1* 


B  -k&.klQ 
B>-25.087 

bI  ^  6.661* 

B^  =-0.831 
Bp  =  0.03U 


A  =  O.OOl* 
A°  =  0.019 
Af  =  0.075 
Ag  =  0.299 
Ao  =  1.016 
\  =  0,ni>8 


Bj^  =50.366 

B, =-77. 957 
b:  =3l*.70l* 

B^  =-6.362 

Bg  =  O.2U9 

B„  =  0,in>9 
m  ’ 


0  v=8.500  y=23.190  E  v=9.000  y=l6.l*38  0  v=8.503  y  =13. 570 


Aq  =  0.087 

An  =  0.022 

a12=  0.02l» 
Ag  =  0.007 
A_2=  0.009 

Aj  “  0 . 0.3O 

A^3=  0.005 

Ai,  =  0.ll*2 
A_4=-0.002 
Ac  =-0.013 
Ag  =-0.001 
A„  =C,n<--l* 


B.  =  6.709 

B^  =-i*.67i+ 
Bg  =  1.195 
Bo  =-0.il3 

B3^q=  0.002 


A^  =  0.033 
A3  =  0.0&l» 
Ac  =  0.268 

A7  =  0.765 

Ap  =-0.ll*9 
4ju  =  0,in>9 


B  =  5.802 
Bg  =  3.126 
Bi,=-11.2U5 
Bg  =  l*.llt7 

bJ  =-0.hl7 

B„  =  0,m>8 


A  =  0.087 

=  0.030 
A_l=  0.032 
A  =  0.016 

A  2“0.056 
A“  =  0.010 
A:3=-0.157 

ftl;' =-0.031 
A-i,=  0.031 


0.727 

1.1*71* 

-l.ll*2 

0.2l*l 

-O.Oll* 


i 


Table  k  -  (continued) 


E  v=9.000  y=8.00U 


v=8.5*t!>  y=5.630  E  v=9.000  y=2.158 


Aq  =  0,098 
Ap  =  0.265 
A4  =  0.1*90 
Ag  =  0.71^1 

A3  =-0.180 

Ajn  =  0,in>8 


1.936 

O.U6O 

-2.027 

0.689 

-0.059 

0,ni>9 


=-0.11*6  =  O.31I* 


Air  =-0.066 
A  ,=-0.522 
Ap  =  0.01*5 
A  p=-0.270 

a”  =  0.078 
A^-=  0.358 
A^  =-0.028 
A  i,=-0.052 

A5  =  0.001 


B?  =  0. 1*1*6 

=-0.292 

Bg  =  0.0l*3 
B£q=-0 ,001 


A  =  0.365 

A^  =  0. 1*1*0 
A^  =  0. 1*1*2 

A'  =-0.280 
Aq  =  0.033 
\  =  0,in>9 


Bq  =  0.521 
Bg  =  0.91*2 
B^  =  0.461* 
Bg  =-0.607 
Bo  =  0.093 

B  =  0,in>8 


0  v=8.787  p =-0.770  E  v=9.000  y =-2.158  0  v =8.980  y=-7.710 


Aq  =— 0 . < 

A,  =-0.( 
A_i=  0.( 
A2  =— 0 .  ( 

A_2~  0 .  ( 

A3  =  0.: 
A_3=-0.( 

Aii  =  0.( 

A.4=  0.( 


B,  =10.313 
B-  =-1,546 
B5  =-2.326 
B„  =  1.138 
Bg  =-0.111 

Bi, =-0.001 


Aq  =  0.521 
A-  =  0.942 
A^  =  0.464 
Ag  =-0.607 

Ao  =  0.093 
A  =  0,in>8 


0.365 

0.440 

0.442 

-0.280 

0.033 

0,in>9 


Aq  =— 0 . 291 

A,  =-0.060 

A  ,=-0.090 
C  =-0.047 
A  p=-0.346 
a'  =  0.064 

Af3=  C.I69 
A^  =-0.016 

A_i4=-0.022 

A^  =  0.001 


Bp  =-0.269 
Bi^  =  0.143 
Bg  =  0.275 
Bg  =-0.092 
B  Q=  0.005 


E  v=9.000  y=-8.004  0  v=9.000  y=-l6.420  E  v=9.000  y=~l6.438 


Ai  =  1.936 
A3  =  0.460 

Ac  =-2.027 

Ay  =  0.689 
A9  =-0.059 
Am  =  0,m>9 


Bq  =  0.098 
Bg  =  0.265 
Bi^  =  0.490 
Bg  =  0.741 
Bg  =-0.180 

Bjn  =  0,m>8 


Aq  =  0.136 
A,  =  0.080 
A  ,=  0.155 
Ap  =  0.718 
A  p=-0.225 
A3  =-0.426 
A_3=  0.073 
Ai,  =  0.080 

A_i4=-0.078 

A5  =-0.004 


=-0.099 
=-0.036 
=  0.060 
=  0.173 

=-0.034 


Aq  =  5.802 
Ap  =  3.126 
Af=-11.245 
Ag  =  4.147 
Ag  =-0.417 
=  0,in>8 


0.033 

0.084 

0.268 

0.765 

-0.149 

0,in>9 


CHAPTER  IV 


SOLUTION  OF  THE  VECTOR  WAVE  EQUATION 

In  this  chapter  the  electric  field  in  a  scurce-free  region 
containing  a  perfectly-conducting  plane  angular  sector  is  determined. 
The  medium  surrounding  the  plane  angular  sector  is  linear,  isotropic, 
and  homogeneous.  Thus  E  satisfies  the  equation 

V  X  V  X  E  -  ic2  E  =  0  (U.l) 

with  the  boundary  condition 

n  X  E  =  0  (1*.2) 

on  the  surface  of  the  plane  angular  sector  and  a  radiation  condition 
as  r  ->  *.  The  \mit  vector  n  is  normal  to  the  plane  angular  sector. 
Since  n  is  not  defined  at  the  edges  and  tip  of  the  plane  angular 
sector,  there  are  aJ.so  conditions  which  must  be  satisfied  there. 

Tliese  will  be  discussed  in  Chapter  VI. 

Following  the  method  described  by  Morse  and  Feshbach[l5]  define 
the  following  three  independent  vector  wave  fvuic.tions. 

L  =  Vi|»  {h.3) 

Mn  =  V  X  '{'2  R 
N]_  =  -  V  X  V  X  ’J'l  R 

K 


(1^.5) 
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where  <l'i,  and  '{'2  are  scalar  wave  functions  and  R  =  r  r  is  a  radius 

vector.  Since  V  <■  E  =  0,  it  is  clear  that  only  the  latter  two  can  he 
solutions  of  equation  (4.1).  Checking  M2  and  in  equation  (4.2), 
it  is  seen  that  <42  must  satisfy  the  Neumann  boundary  condition  and 
the  Dirichlet  ho?J,ndary  condition.  Thus  'I'l  and  4>2  are  the  complete  sets 
of  scalar  wave  fixnctions  detemined  in  Chapter  III.  Since  the  sets 
'f'l  and  ^>2  are  complete,  the  sets  M2  and  Ni  are  complete,  and  thus  any 
divergeneeless  vector  can  be  expanded  as  a  sum  of  these  functions , 
Performing  the  indicated  curl  operations  in  equations  (4.4)  and 
(4.5),  the  vector  wave  functions  can  be  written  in  component  forms  as 


Zv£2 

k^  sin^  e-fk'S  sin^  ^ 


1  -  cos^  6 


(‘•.6) 


(•J-)  <i> 


1  -  k*^  cos^ 


_  zvjii  (icr) 

V  “  ‘'■''I''*’  — 


(h.l) 


(r  zv£.i  (icr))' 


ter  [k^  sin^  0  +  k'2  sin2  ({,  ^ 


1  -  k’2  cos2  (|,  ©e^fe) 


Sill' 


4-6, •,(*?>)  <l>  +  |l  -  k2  cos2  6  ©1(6)  (<},)  0 

''  '■  0*1  0*1 


The  subscript  Z  is  an  ordering  index.  It  is  an  integer  which 
identifies  each  (v^,  y^)  eigenvalue  pair,  and  is  useful  not  only  for 
identification  but  as  a  summation  index. 
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(icr)  is  a  spherical  Bessel  function,  its  type  depending  on  the 
bound8u:y  condition  on  r.  All  primes  indicate  differentation  with 
respect  to  the  independent  variables  r,  6,  and  (f. 

Note  that(v=0,  ;j=0)  is  an  eigenvalue  pair  for  the  scalar  wave 
functions,  but  it  is  not  for  the  vector  wave  functions.  Since 

®£2 

contains  derivatives  of  the  scalar  functions  and  the  scalar 
eigenfunctions  that  correspond  to  (v=0,  y=0)  are  constant, 

Me2(v=0,  paO)  vanishes. 

It  is  aJ.so  possible  to  define  two  other  vector  wave  functions. 
Mg  ,  and  Ng  These  do  not  satisfy  equation  (4.2)  so  they  are  not 
acceptable  as  solutions  of  this  problem.  They  are  useful,  however, 
and  will  be  used  in  the  dyadic  Green's  function  derived  in  the  next 
chapter.  These  functions  are  defined  as  follows: 


M. 


|£1  gj^jj2  9  +  jji2  sin2  ^ 


1  -  cob2  0  o'  (e)  (4.8) 


|£1 


V'*’ '  ' 


1  -  k'2  cos2  <j)  ® 


v‘®’  Vs'*’ " 


(r 


Kr  k^  sin^  0  +  k'2  sin2  (}i 


1  -  k'2  cos^  ij) 


(4.9) 


0ej^2(0)  ^  I  1  -  cos2  6  0ej^2(e)  e 

0  0 


3£2'  '  e£2' 


77 


These  are  related  to  M.  and  N  in  tbc  following  manner: 


M-  =  i  7  X 
§11  K 


N-  =  —  7  X 
§£2  >c  ^12 


=  ±  7  X  N 

§12  ,c  g12 


"eu  ‘  ^  ^  -  «eu 


(It.lO) 


(4.12) 


(4.13) 


All  four  of  the  vector  wave  functions  defined  by  equations 
(4.6),  (4.7)»  (4.8),  and  (4.9)  possess  orthogonality  properties  on  the 
surface  of  a  sphere.  These  will  be  discussed  in  the  next  chapter, 
where  they  are  used  to  find  the  dyadic  Green’s  function.  To  facilitate 
the  investigation  of  orthogonality,  the  following  auxiliary  vector  wave 
functions  are  introduced  at  this  point. 


1  -  cos2  e 


§11  Ik^  sin2  e  +  k’2  sin2  4,  "gu' 


1  -  k'2  cos2  (J>  /  V  I 

- 0  (e)  ^  u)  e 

k2  sin2  0  +  k’2  sin2  !()  °  ■  ° 


(4.14) 


r  X  m- 


|l  -  k’^  cos2  (j) 


C,(4>)  J  (4.15) 


k2  sin2  6  +  k’2  sin2  ^  §1^ 


|l  -  k2  cos^  8  t  ,  .  ,  ,  - 

^7-''  =  0 

k^  sin2  0  +  k'2  sin2  (Ji  ° 


-  -aii? 
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O  o^-*-  O*'* 

The  auxiliary  vector  wave  functions  m  »  r  x  nig  ,  Tg  are  defined 

q12  qJ,2  g£2 

similarly  except  that  Neumann  functions  are  used  instead  of  the 
Dirichlet.  These  functions  are  related  to  the  original  vector  wave 
functions  in  the  following  manner. 


=  Zy  (xr)  m- 

g£2  '^£,2  g£2 


(U.17)- 


N 


g£l 


(ter) 

icr 


U 


(r  Zyj^,  (xr))'  _ 

+ -  (r  X  m.„, ) 


xr 


(Ii.l8) 


and  similarly  for  Mg  and  Ng  . 

q£1  q£2 

For  a  more  complete  discussion  of  vector  wave  functions ,  see 
Morse  and  Feshbach[l5] .  Vector  wave  functions  are  also  discussed  by 
Stratton [ l6 ] ,  Van  Bladel[l7],  and  Spence  suid  WellsflS], 


- 


If 


i  J 

;)  ( 


fr  ^ 


CHAPTER  V 

DETERMINATION  OF  THE  DYADIC  GREEN'S  FUNCTION 

In  this  chapter,  the  dyadic  Green’s  function  is  derived.  The 
derivation  is  given  in  some  detail  since  the  Green's  function  is 
composed  of  previously  tinknown  *  ector  wave  functions.  In  the 
process,  some  properties  of  these  vector  wave  fuifictions  are  determined. 

Introduction 

In  order  to  find  the  total  field  due  to  an  arbitrary  current  J 
in  the  vicinity  of  the  perfectly  conducting  plane  angular  sector,  the 
following  set  of  equations  must  be  solved. 

7xVxE-k2e  =  ..  (5.1) 

n  X  E  =  0  (5.2) 

on  the  plane  angular  sector,  E  must  satisfy  the  radiation  condition 
as  r  ->  «,  and  as  was  mentioned  in  Chapter  IV,  there  are  also  edge  and 
tip  conditions  which  must  be  satisfied. 

The  solution  to  this  set  of  equations  can  be  written 

E  (R)  =  J^r  (R,R’)  •  J  (R')  dv  (5.3) 

+  /  [E  (R’)  X  n  •  V  X  F  (R',R)]  ds 
s 
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vhere  R  is  the  field  point,  R'  is  the  source  point,  and  v  is  a  volume 
enclosing  the  source  current  J .  The  surface  S  is  the  plane  angular 
sector.  Ordinarily  the  surface  integral  is  zero  because  of  equation 
(5-2).  However,  if  the  plane  angular  sector  contains  slots,  n  x  E 
will  not  vanish  in  the  apertiire  of  the  slots.  Thus  the  dyadic  Green's 
fxmction  obtained  in  this  paper  can  be  used  not  only  to  determine  the 
fields  radiated  b^''  an  arbitrary  cuirent  distribution  in  the  vicinity 
of  a  plane  angular  sector,  but  also  to  determine  the  fields  ra’.iated 
by  slots  in  a  plane  angular  sector,  r  (R,R')  is  the  dyadic  Green's 
function.  It  must  satisfy  the  following  equations. 

V  X  V  X  r  (R,R')  ~  k^T  (R,R')  =  -  'F  6  ( |R-R’ j )  (5.U) 

where  e  is  the  unit  dyadic. 

n  X  r  (R,R’)  =  0  (5.5) 

on  the  plane  angular  sector,  r  (R,R')  must  satisfy  a  radiation 
condition  as  r  •♦•  »,  and  P  (R,R')  must  satisfy  conditions  at  the  edges 
and  tip. 

In  order  to  make  the  derivation  of  P  (R,R')  as  simple  as  possible, 
consider  the  following  operation  on  equations  (5.**)  and  (5.5).  Let 
an  arbitrary  unit  vector,  a,  be  dot  multiplied  to  the  right  of  the 
two  equations.  The  result  is 

V  X  V  X  G  ('r,R'  )  -  k2  g  (R,R'  )  =  _  a  6  (  |R-R'  | )  (5.6) 

n  X  G  (R,R')  =  0  (5.7) 

G  (R ,R ' )  =  Y  (R,R ' )  •  a 


where 


(5.8) 
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The  vector  Green's  function  G  (R,R*)  is  the  field  due  to  a  unit 
vector  point  source  a  located  at  R'.  Once  G  (R,R')  is  determined 
it  can  be  compared  with  equation  (5*8)  to  determine  T  (R,R'). 

It  was  determined  in  the  last  chapter  that  an  arbitrary 
divergenceless  vector  function  could  be  represented  as  a  sum  over 
the  complete  set  of  vector  wave  functions.  Mg  ,  and  .  Thuc 

0^2  gtl 

G  (R,R' )  can  be  written 

G  (R,R*)  =  I  [aq  (R’ )  Mqg  (R)  +  bq  (R ’ )  (R) ]  (5-9) 

q 

where  the  subscript  q  replaces  the  e^  subscript  used  previously,  and 
is  meant  to  include  all  of  these  functions.  The  subscript  1  or  2 
means,  as  before,  that  Dirichlet  or  Neumann  functions,  respectively, 
are  used  in  the  vector  wave  functions.  Equation  (5*9)  automatically 
satisfies  equation  (5.7).  The  coefficients  aq  (R')  and  bq  (R*)  will 
be  determined  in  order  to  satisfy  equation  (5.6). 

Derivation 

The  singularity  at  R'  separates  space  into  two  regions.  Let  the 
dividing  surface  S'  be  a  sphere  with  radius  r'  except  for  an  infinitesi¬ 
mal  cut  around  the  quarter  plane.  (See  Figure  7)  The  volume  within  S' 
is  called  region  I  and  that  outside  of  S'  is  called  region  II.  Also 
define  a  surface  S  with  the  same  shape  as  S'  but  at  a  radius  r.  The 
surface  S  is  in  region  II.  V/ithin  region  I,  the  vector  Green's 
fxinction  has  the  boundary  condition  on  r  that  the  function  must  be 
regular  at  r  =  0.  This  implies  that  the  Bessel  functions  used  in  the 
vector  wave  functions  must  be  of  the  first  kind.  Thus  in  region  I, 


Regions  used  JLn  the  determination 
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G-  (R,R*)  =  l  [&q  (R')  Mq2  *  \  ^ql 

q 

vhere  the  superscript  I  means  that  (sr)  =  jy  (Kr).  In  region  II, 
the  boundary  condition  on  r  is  that  the  function  must  represent  an 
outward  propagating  wave  for  large  r.  This  implies  that  the  Bessel 
functions  must  be  Henkel  functions  of  the  second  kind.  Thus  in 
region  II 

(R,R')  =  J  [aq  (R‘)  (R)  +  bq  (R ' )  (R)]  (5.11) 

<1 

{ 2 ) 

where  the  superscript  II  means  that  {<t)  =  h^  '  (Kr). 

In  order  to  determine  the  unknown  coefficients  aq  (R*)  and  bq  (R'), 
consider  the .following  equation,  known  as  Green’s  second  identity  for 
vectors. 

(B*Vx7xA-A*Vx7xB)dv  (5.12) 

=/  (Ax7xB-Bx7xA)*nds 
s 

A 

where  n  is  a  unit  normal  to  the  surface  S  which  is  pointing  out  of 
the  enclosed  volvune  V.  Let  the  surface  be  the  same  as  the  surface  S 
in  Figure  7.  Then  V  includes  all  of  region  I  and  part  of  region  II. 

Let 


A  =  G  (R,R') 

(5.13) 

B  .  mJj  (B) 

(5.11*) 
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The  volume  integral  is  readily  evaluated.  From  equations  (5.6) 
and  (5.13), 

Vx7xA=:VxVxG  (R,R' )  =  G  (R,R‘ )  -  a  6  (  |R-R'  \ )  (5.15) 


Prom  equations  (4.1)  and  (5.1^0, 


7  X  7  X  B  =  V  X  V  X  (R)  =  <2  (R) 

P2  p2 

Thus  the  left  hand  side  of  equation  (5*12)  is 

/  ‘  °  )  -  Mp2  *  a  5  (  |R-R’  1 ) 


(5.16) 


(5.17) 


-  G  (R,R’)  •  M^2  dv  =  -  a  •  (R’) 


The  evaluation  of  the  surface  integral  is  not  quite  as  simple. 
Consider  the  following  vector  identities. 

Ax7xB*n  =  n  •Ax7xBs=nxA'7xB  (5.18a) 

Bx7xA  •n  =  n  •Bx7xA  =  nxB  •7xA  (5.18b) 


Using  these  in  the  right  hand  side  of  equation  (5.12),  the  following 
is  obtained. 

/  [n  X  G  (R,R')  •  7  X  if 5  (R)  (5.19) 

's 

-  n  X  i^2  *7x0  (R,R')]  ds 

On  the  part  of  S  that  encloses  the  quarter  plane 

n  X  G  (R,R')  =  0  (5.20) 

and 

n  X  M^2  =  0  (5.21) 


so  that  the  integration  is  over  the  sphere  with  radius  r.  On  this 
sphere 
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n  X  G  (R,R’)  =  r  X  (R,R')  =  I  [a^  x  (R))  (5.22) 

+  hq  (R')  {r  X  Ig  (R))3 

n  X  Mp2  (R)  =  ^  X  j^2  (5.23) 

and  with  the  help  of  equations  (^.lO)  and  (U.ll) 

y  X  f^2  ^5*2**^ 

V  X  G  (R,R*)  =  7  X  G^^  (R,R')  (5.25) 

*  tc  n«q  +  \ 

q 

Using  these  results,  the  svirface  integral  is 


(5.26) 

+  bq  (R»)  (r  X  nJJ  (R)  •  (R)  -  r  x  (R)  •  (R)]}  ds 

=  K  ^aq  (R')  Aq  +  bq  (R')  Bq] 

q 

First  consider  Bq 


Bq  =  /  [r  X  (R)  •  (R)  -  ^  x  (R)  •  (R)]  ds 


(5.27) 
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Using  the  axixiliary  vector  wave  functions  introduced  in  the  previous 
chapter,  Bq  can  he  written 

(r  hj^^  (ter))*  (r  (icr))‘  +  r^  (ter)  hv  ,  (ter) 

B_  =  -  _ Pf _ qi 

ic2  1 

-  ,J 

(5-28) 

■C  ^  ^  P  d  (}.d  0  =  f  (r)  Cq 


where  the  weight  function  p  is 


_  heh({)  sin^  6  •<•  k*2  sin^  t}i 


^  1  -k2  cos2  6  1  -  k*^  cos^ 


(5.29) 


Using  the  definitions  of  the  auxiliary  wave  functions,  Cq  is 


Cq  =  /q  ll  [e4i(e)  0p2(e)  $*2(4.) 


(5.30) 


-  eqi(0)  $^3^(<^)  0^2^®' 

If  the  ql  functions  and  p2  functions  are  both  even  or  hoth  odd,  the 
integrand  is  odd  with  respect  to  4>  snd  Cq  =  0.  Thus  the  ql  functions 
must  he  even  and  the  p2  functions  odd  or  vice  versa  if  Cq  is  to  he 
nonzero.  Separating  the  0  and  ^  integrals  Cq  can  he  written 

Cq  =  $qi(<|))  /J  0’i(o)  0P2(e)  <10 

-  %(<fr)  <tp2(<l')  d4>  0qi(0)  0^2(0)  ^0 


(5.31) 
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Integrating  the  first  0  integral  hy  parts,  the  following  is  obtained. 


(5.32) 


-  /q  ®ql(®)  0p2(®)  <1® 

The  end  point  contributions  are  zero  because  0qi('Jr)  =  0,  and  either 
0ql(6)  or  6p2(9)  is  odd  so  that  the  contribution  at  0  =  0  is  zero. 
Combining  terms,  Cq  can  be  written 

Cq  =  -  /o  eql(0)  0p2^®)  ^8  [f  (♦ql^'<‘)  ^5-33) 


1  ^ 

The  term  in  the  square  brackets  is  just  '5>qi(4>)  ^p2(^)|  ^  which  is 
zero  because  of  the  periodicity  of  the  ♦((Jj)  functions.  Thus  Cq  =  0 
and  it  follows  that  Bq  =  0  for  all  q. 

Now  consider 

Aq  =  /g  [?  X  iig  (R)  •  (r)  -  f  X  (R)3  <^5  (5-34) 


Again  using  the  auxiliary  vector  wave  functions  and  manipulating  the 


dot  and  cross  products,  Aq  can  be  written 


Aq  = 


.(2) 


hy^g  (xr)(r  jy^g  (icr))'  -  jy^g  (icr)(r  h^|g  (xr) ) ' 


icr 


(5.35) 
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Using  the  definition  of  m,  the  surface  integral  can  he  written 

L  'v  ■  v'  °  4 


0^2(6)  0’2(e) 


(he)- 


0p2(e)  ♦;2(*)  0^2(6) 

{h^)2 


ds 


(5.36) 


where  the  metric  coefficients  have  been  used  for  compactness.  In 
order  to  evaluate  equation  (5*36)  a  slight  digression  is  necessary. 
Consider  the  following  vector  identity. 


V  •  (’4-p2  V  -J-qg)  =  V  ^  V  '  ^ 


q2 


(5.37) 


Integrate  throughout  a  volume  like  region  I  in  Figure  T,  use  Gauss' 
theorem  and  ♦q2  =  -  *q2  to  obtain 

/g'l'p2  V  )J)q2  •  nds  =  -  <2  'l>^2  %2  '  '^q2 


.38) 


Expanding  the  functions  into  component  form,  the  left  hand  side  can  be 
written 


=  Rp2(r)  Rq2(r)  4®p2^®^ 


where  the  contribution  to  the  integral  along  the  quarter  plane  is  zero 
and  the  surface  S  is  now  spherical.  The  surface  integral  is  the 
orthogonality  relation  for  the  Lame  functions  so  it  is  zero  unless 
q  =  p,  and  the  left  hand  side  of  equation  (5«38)  is  zero  unless 
q  =  p.  The  first  terra  on  the  right  hand  side  of  equation  (5-38)  can 


be  written 


4  '^q2  V  "4  ^4  V  V 
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(5.1<0) 

The  reason  the  voliame  integral  can  be  separated  in  this  way  is  that 

dv  =  hr  dr  (he  h^  d9  d(J))  =  hr  dr  ds  =  dr  ds  (5.J*l) 

The  metric  coefficient  hr  is  not  a  function  of  0  or  (|)  so  the  surface 
integral  can  be  integrated  independently  of  r.  But  the  surface  integrsQ. 
is  the  orthogonality  relationship  again  so  it  is  zero  uni  ess  q  =  p,  and 
thus  the  first  term  on  the  right  hand  side  of  equation  (5*38)  is  zero 
unless  q  =  p.  The  conclusion  to  be  drawn  at  this  point  is  that  the 
gradient  of  the  Lame  functions  possess  an  orthogonality  relationship 
over  a  volume,  i.e. , 

ly  ^  ^  ^ci2  ~  °  ^  ^  P  (5.1»2) 

Again  using  the  condition  on  the  metric  coefficients,  this  can  be 
written 

/  7  I^p2  •  V  Ipq2  dv  =  /^  (  /g  V  i|<p2  *  V  ds)  dr  (5.1*3) 

For  q  9^  p,  the  volume  integral  is  zero,  independently  of  the  radius 
of  the  spherical  volume.  This  implies  that  the  surface  integral  must 
be  identically  zero,  and  thus  the  gradient  of  the  Lame  functions  also 
possess  an  orthogonality  relationship  on  the  surface  of  a  sphere,  i.e., 

4  (V  >lp2  •  V  t|;q2  )  ds  =  0  q  #  p  (5.41*) 

Expanding  the  gradients  in  component  form,  this  can  be  written 
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4  ^  ^  '*'p2  •  ^  )  as  =  Rp2(r)  /s®p2(6)  '^'pgU)  0q2^®)  '*’q2('J’^  ^s 


-V<')B,,(r)4  v<*)  yu) 

L  (Oe)^ 

^  V<^>  ^p2t*>  V(8) 

(h*)2 


(5.1*5) 


Again  the  first  term  on  the  right  hand  side  is  the  orthogonality 
relationship  for  the  Lame  functions.  Thus  it  has  "been  shown  that  for 


q  Pj 


To;;(e)  04^(6)  ^  9^^(9)  r^(»)  e^^(e) 

(ha)^  (h^,)^ 


=  ^2  4  *  “q2^  ° 


(5.1*6) 


and  that  the  auxiliary  vector  wave  functions  mq2  are  orthogonal  on 
the  surface  of  a  sphere.  Thus  in  equation  (5.35)  Aq  can  he  written 


0  q  5^  p 


hj,p2  (''r)(r  ~  l^Vp2^''^))'  g 

An  ~  - -  -  - —  . .  ■  - -  —  r 


r  o^gie)  ^  ep;(e) 

®  he  h^  ^ 


(5.1*7) 


ds  q  =  p 
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Noting  that 


"  ^'’p2  ^Vp2 

Kr 


(5.U8) 


4p2  Jvp2  -  ^Vp2  ^ipg’ 


.(2) 


vtcr)2 


and  writing  out  tlii  metric  coefficients,  it  is  seen  that 


.TT  ,V 


A.  =  -^  /  / 

^  <2  U  ^0 


1  -  k2  cos2  e  ♦p2(*))^ 


1  -  k'^  cos^  (() 


(5.49) 


1  ~  cos^  <Ji 


1  -  k^  coa^  6 


d0  d(j> 


and  thus 


(!■) 

Ap2 


where  Ap2  is  the  surface  integral  in  equation  (5*52) 


(5.50) 


Ap2  =  -  ^  Ap 


(5.51) 


Repeating  the  above  procedure  starting  with  equation  (5.12)  and 
letting  B  =  Np-]^  (R)  instead  of  1^2  (^)»  I’p  (^')  is  determined  to  be 


hp  (R')  = 


i  K  wJi  (R')  -a 


Tl 


(5.52) 
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where  Ap^  has  the  same  form  as  Ap2,  the  only  diTference  being  that 


Dirichlet  functions  are  used  instead  of  Neumann.  Thus, 


(R,R')  =6^1 
q 


Mg  (K)  Mg  W)  -  I 


"»  ”51 


(5.53) 


Comparing  this  with  equation  {5*8)»  it  is  seen  that 


r  (R,R')  =  i  K  I 


mJs  (R)  M^gCR’) 


"S  > 


|r1  >  |r’ 


(5.5l») 


where,  with  a  =  1  or  2, 


IT  IT  |l  -  6 

=  /,  /„  . p 

1  -  k'2  cos^  (j> 


1  -  k'^  cos^  (p 


1  -  k^  cos^  0 


(0q„(0) 


(5.55) 


To  find  F  (R,K’)  for  |r1  <_  |k'|,  the  same  procedure  could  be  repeated, 
this  time  integrating  throughout  a  volume  external  to  the  surface  S, 
and  letting  the  surface  S’  be  at  R',  where  |R'1  ^  jRj,  or  the 


symmetry  property  of  Green’s  functions  could  be  used.  The  result  is 
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r  (R,R')  =  j  K  I 


-I 

^q2 


(R)  Mq2  (RO 


^q2 


(5.56) 


(R)  nJ^  {R‘) 

H-  - 1 - 

Aql 

A  point  concerning  the  normalization  constants  is  worthy  of 
note.  It  is  shown  in  Appendix  A  that  the  occur  5.n  discussing 
the  properties  of  the  two-dimensional  Sturm-Liouville  type  Lame 

•V- 

operator.  Using  the  notation  in  the  Appendix 


|r|  |r' 


^qa  ~  ^  yqct*  ^  ^qa  ^ 


But, 


^  yqa»  I-yqa  ^  , 

“  ^qa 

^  yqa»  P  yqa  > 


(5.57) 


(5.58) 


Thus,  if  the  normalization  of  the  Lame  functions  had  "been  chosen 
such  that 

^  yqo»  P  yqa  ^  ~  1  (5-59) 

i.e.,  if  the  Lame  functions  had  been  orthonormalized,  then 


^qa  ”  ^qa  ~  '^qa  ^'’qa  '**  (5.60) 

From  the  point  of  view  of  the  applied  mathematician,  this  would  be  an 
acceptable  normalization  of  the  Lame  functions.  This  is  especially 
true  for  the  Lame  polynomials.  Much  confusion  is  caused  by  the  lack 
of  a  standard  when  reading  the  work  of  various  authors.  The  same 
problem  exists  with  the  Legendre  polynomials.  They  are  well  defined. 


9l« 

however  they  ai*e  not  orthonormal.  This  is  the  reason  for  the  rather 
complicated  normalization  constant  found  in  eigenfunction  expansions 
in  spherical  coordinates. 

The  normalization  constants,  A  ,  are  given  in  Table  6.  In 

qa 

order  to  orthonormalize  the  Lame  functions  presented  in  Tables  1 
through  U,  it  is  necessary  to  multiply  each  0  and  ^  function  by 
(  These  numbers  are  also  given  in  Table  6. 


Table  6  -  Normalization  Constants 
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A;  Dirichlet  Problem 


V 

V- 

Ai 

0.296 

0.090 

2.940 

0.601 

1.425 

0.915 

35.103 

0.560 

1.000 

0.000 

8.378 

0.699 

1.130 

-  0.455 

14.201 

0.642 

2.480 

2.670 

269.941 

0.423 

2.000 

1.500 

10.050 

0.879 

2.290 

0.215 

20.797 

0.776 

2.000 

-  1.500 

10.050 

0.879 

2.o4o 

-  1.705 

194.581 

0.423 

3.495 

5.440 

3073.608 

0.267 

3.000 

3.873 

19.175 

0.889 

3.410 

1.535 

67.639 

0.687 

3.000 

0.000 

5.741 

1.202 

3.145 

-  0.825 

37.245 

0.769 

3.000 

-  3.873 

19.175 

0.889 

3.010 

-  3.940 

1012.967 

0.330 

4.499 

9.225 

26302.376 

0.175 

4.000 

7.190 

48.982 

0.799 

4.470 

3.790 

370.670 

0.507 

4.000 

2.190 

6.276 

1.336 

4.280 

0.335 

39.904 

0.867 

4.000 

-  2.190 

6.276 

1.336 

4.050 

-  2.575 

153.  U4 

0.605 

4.000 

-  7.190 

48.982 

0.799 

4.005 

-  7.205 

4648.698 

0.256 

5.500 

14.013 

206778.639 

0.115 

5.000 

11.489 

157.321 

0.661 

5.490 

7.110 

2517.012 

0.345 

5.000 

5.196 

10.424 

1.302 

5.499 

2.110 

136.722 

0.715 

5.000 

0.000 

4.265 

1.629 

5.150 

-  1.170 

62.106 

0.845 

5.000 

-  5.196 

10.4?4 

1.302 

5.015 

-  5.335 

970.112 

0.420 

5.000 

-11.489 

157.321 

0.661 

5.000 

-11.493 

73566.927 

0.142 

Table  6  -  (continued) 
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V 

Ax 

6.500 

19.805 

164299.188 

0.131 

6.000 

16.783 

533.201 

0.530 

6.499 

11.455 

11646.480 

0.254 

6.000 

9.114 

23.400 

1.157 

6.465 

4.850 

1172.142 

0  450 

6.000 

2.832 

4.679 

1.731 

6.282 

0.445 

1300.687 

0.433 

6.000 

-  2.832 

4.679 

1.731 

6.055 

-  3.380 

470.909 

0.549 

6.000 

-  9.114 

23.400 

1.157 

6.005 

-  9.160 

4881.640 

0.305 

6.000 

-16.783 

533.201 

0.530 

6.000 

-16.783 

427489.535 

0.099 

7.500 

26.597 

523205.875 

0.105 

7.000 

23.075 

841.478 

0.508 

7.500 

16.820 

62918.  UO 

0.178 

7.000 

14.000 

62.554 

0.973 

7.492 

8.695 

411996.648 

0.111 

7.000 

6.444 

7.341 

1.662 

7.393 

2.660 

1053.201 

0.493 

7.000 

0.000 

3.448 

2.007 

7.157 

-  1.500 

1947.103 

0.416 

7.000 

-  6.444 

7.341 

1.662 

7.020 

-  6.660 

1614.970 

0.432 

7.000 

-14.000 

62.554 

0.973 

7.000 

-14.010 

8308.581 

0.287 

7.000 

-23.075 

841.478 

0.508 

7.000 

-23.075 

992383.008 

0.087 

8.500 

34.389 

203797.938 

0.i4l 

8.000 

30.367 

5077.147 

0.345 

8.500 

23.191 

115954.960 

0.162 

8.000 

19.876 

156.869 

0.823 

8.497 

13.590 

866869.250 

0.098 

8.000 

10.949 

14.942 

1.482 

8.463 

5.885 

2928.812 

0.407 

8.000 

3.441 

3.745 

2.094 

8.275 

0.535 

464459.938 

0.113 

8.000 

-  3.441 

3.745 

2.094 

8.073 

-  4.170 

758.318 

0.557 

8.000 

-10.949 

14.942 

1.482 

8.005 

-11.015 

1130.930 

0.502 

8.000 

-19.876 

156.869 

0.823 

8.000 

-19.879 

6703.902 

0.322 

8.000 

-30.367 

5077.147 

0.345 

8.000 

-30.367 

902867.680 

0.09^' 

Table  6  -  (continued) 
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V 

n 

Al 

Oa/Ai: 

9.500 

9.000 

38.660 

677.440 

0.604 

9.500 

— 

— 

— 

9.000 

26.751 

645.618 

0.611 

9.50 

— 

— 

— 

9.000 

16.431 

37.872 

1.242 

9.000 

T.SkO 

5.650 

1.998 

9.000 

0.000 

2.915 

2.357 

9.000 

-  7.640 

5.650 

1.998 

9.000 

-16.413 

37.872 

1.242 

•9.00 

~  -16.4 

«  - 

9.000 

-26.751 

645.618 

0.611 

9.00 

~  -26.8 

-- 

9.000 

-38.660 

677.440 

0.604 

9.000 

-38.660 

607247.125 

0.110 

B :  Neumann  Problem 

1.000 

0.500 

16.750 

0.588 

0.814 

-  0.194 

12.233 

0.589 

1.000 

-  0.500 

16.750 

0.588 

2.000 

1.732 

120.627 

0.472 

i.595 

0.795 

9.843 

0.805 

2.000 

0.000 

20.094 

0.739 

1.955 

-  1.552 

12,890 

0.818 

2.000 

-  1.737 

120.627 

0.472 

3.000 

3.950 

1058.871 

0.326 

2.520 

2.625 

14.021 

0.892 

3.000 

0.950 

45.070 

0.718 

2.803 

-  0.350 

6.998 

0.111 

3.000 

-  0.950 

45.070 

0.718 

2.990 

-  3.890 

22.573 

0.853 

3.000 

-  3.950 

1058.871 

0.326 

4.000 

7.211 

9332.023 

0.215 

3.505 

5.430 

30.317 

0.850 

4.000 

2.646 

178.674 

0.578 

3.617 

1.260 

6.567 

1.263 

4.000 

O.OGJ 

40.897 

0.836 

Table  6  -  (continued) 
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V 

V- 

Ai 

(Ai/Ai)v 

3.940 

-  2.315 

6.817 

1.300 

4.000 

-  2.646 

178.674 

0.578 

3.998 

-  7.195 

56.527 

0.771 

4.000 

-  7.211 

9332.023 

0.215 

5.000 

11.494 

76l8i<.075 

0.141 

4.500 

9.218 

85.544 

0.733 

5.000 

5.363 

1047.489 

0.411 

4.530 

3.680 

8.176 

1.323 

5.000 

1.369 

73.363 

0.800 

4.795 

-  0.495 

5.281 

1.515 

5.000 

-  1.369 

73.363 

0.800 

4.904 

-  5.230 

10.556 

1.296 

5.000 

-  5.363 

1047.489 

0.411 

5.000 

-11.491 

176.504 

0.642 

5.000 

-U.494 

76184.075 

0.i4l 

6.000 

16.784 

617657.477 

0.090 

5.500 

14.012 

284.757 

0.595 

6.000 

9.165 

6926.720 

0.279 

5.507 

7.060 

15.340 

1.236 

6.000 

3.507 

236.801 

0.649 

5.627 

1.692 

4.679 

1.680 

6.0C0 

0.000 

62.880 

0.904 

5.927 

-  3.015 

4.486 

1.739 

6.000 

-  3.507 

236.801 

0.649 

5.997 

-  9.120 

26.646 

1.120 

6.000 

-  9.165 

6926.720 

0.279 

6.000 

-16.783 

637.161 

0.507 

6.000 

-16.784 

617657.477 

0.090 

7.000 

23.076 

70189.938 

0.168 

6.500 

19.804 

1031.706 

0.466 

7.000 

14.014 

46539.14? 

0.186 

6.501 

11.440 

69.934 

0.914 

7.000 

6.708 

1135.994 

0.471 

6.540 

4.675 

5.719 

1.714 

7.000 

1.770 

102.166 

0.860 

6.792 

-  0.630 

169.150 

0.748 

7.000 

-  1.770 

102.166 

0.860 

6.980 

-  6.550 

4.261 

i.901 

7.000 

-  6.708 

1135.994 

0.471 

7.000 

-14.002 

188.385 

0.738 

7.000 

-14.014 

46539.142 

0.186 

7.000 

-23.075 

2412.174 

0.390 

7.000 

-23.076 

701895.938 

0.094 

Table  6  -  (continued) 
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V 

Ai 

(Ai/Ai)^ 

8.000 

30.368 

943368.000 

0.093 

7.500 

26.597 

3927.403 

0.357 

8.000 

19.880 

294906.848 

0.125 

7.500 

16.816 

4.848 

1.904 

8.000 

11.036 

6352.965 

0.326 

7.510 

8.625 

9.765 

1.599 

8.000 

4.334 

292.966 

0.704 

7.632 

2.095 

94.074 

0.915 

8.000 

0.000 

85.737 

0.957 

7.923 

-  3.685 

6.299 

1.830 

8.000 

-  4.334 

292.966 

0.704 

7.995 

-10.965 

5.347 

1.915 

8.000 

-11.036 

6352.965 

0.326 

8.000 

-19.877 

819.559 

0.544 

8.000 

-19.880 

294906.848 

0.325 

8.000 

-30.368 

iI<267.483 

0.267 

8.000 

-30.368 

943368. COO 

0.093 

9o000 

38.660 

221153.000 

0.142 

8.500 

34.389 

11763.660 

0.288 

9.000 

26.752 

510399.188 

0.115 

8.500 

23.190 

462.985 

0.646 

9.000 

16.438 

36518.688 

0.223 

8.503 

13.570 

36.122 

1.223 

9.000 

8.004 

1233.913 

0.520 

8.545 

5.630 

36.915 

1.219 

9.000 

2.158 

131.879 

0.909 

8.787 

-  0.770 

2369.904 

0.436 

9.000 

-  2.158 

131.879 

0.909 

8.980 

-  7.710 

8.552 

1.799 

9.00c 

-  8.004 

1233.931 

0.520 

9.000 

-16.420 

9.861 

1.738 

9.000 

-16.438 

36518.688 

0.223 

9.000 

-26.751 

86.694 

1.009 

9.000 

-26.752 

510399.188 

0.115 

9.000 

-38.660 

78706.980 

0.184 

9,000 

- 38.660 

221153.000 

0.i42 

CHAPTER  VI 


DISCUSSION  OF  THE  FIELDS  AND  CURRENTS 
FOR  SOME  SPECIAL  CASES 


This  chapter  contains  several  numerical  examples  of  the  exact 
solution  given  in  Chapter  V.  First,  the  exact  fields  and  current  for 
an  infinitesimal  dipole  source  are  derived.  These  are  then  evaluated 
in  the  vicinity  of  the  tip  of  the  quarter  plane.  Several  source 
locations  and  orientations  are  considered.  The  dominant  "behavior 
near  the  tip  of  a  wide  angle  sector  and  a  narrow  angle  sector  are 
also  discussed.  The  reciprocal  case  of  a  source  close  to  the  tip  of 
the  quarter  plane  is  also  considered,  and  several  far  field  patterns 
are  given. 

Fields  and  Currents  for  a  Unit  Dipole  Source 

The  dyadic  Green's  functions  derived  in  Chapter  V  are  to  "be  used 
in  equation  ( 5 • 3)  in  order  to  determine  the  vector  E  field  at  R  due  to 
a  current  source  at  R'.  If  the  source  is  a  unit  dipole  source  at  Rq, 
i . e . ,  if 

juu  J  (R')  =  a  6  (R'-Ro)  (6.1) 

where  a  is  just  a  unit  vector  at  R^,  the  E  field  is  simply 
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r  —I  _  « 

1  _TT  —  ^q2  ^ 

3kI  (b)  {-2^ - ) 

q  I  Aq2 

-JI  ,  \l  ^  , 

+  N  1  (R)  { -  )  r  ^  r, 

^ql 


(6.2a) 


E  (R)  = 


i  <1  M  2  (R)  {—3 -  ) 


_T  —  ^ql  ( ^ 

+  (R)  (  jj—  )  r  < 


(6.2b) 


The  H  field  is  easily  determined  using  Maxwell’s  equations, 
and  the  properties  of  the  M  and  H  vectors  given  by  equations  (^.lO)  - 


(U.13). 


f— II  _  ^0?  ^^o)  * 

-  K  Yo  I  Nq2  (R)  (—— - 

q  ^  h2 


_TT  -  \l 

+  mJ^  (R)  - 

ql  Aqi 


)  r  >  r<, 


(6.3a) 


H  (R)  =  —  V  X  E  (R)  =' 

uiy 


M?p  (Ro)  •  a 


„  _I  ,  «q2 

Yc  I  Nq2  (R)  (-  --Y- 
q  h2 


nJi  (Ro)  -a 


-I  -  ‘'■ql 

.  m’,  (B)  (-^ 


)  r  < 


(6.3b) 


Yo  is  the  adraittance  of  the  medium 


Y 


o  ~ 


(6.1») 


The  current  on  the  quarter  plane  can  be  found  by  evaluating  the 
tangential  H  field  at  the  quarter  plane.  In  order  to  do  this  it  is 
necessary  to  consider  the  even  and  odd  parts  of  the  H  field. 

Equation  (6.3a)  is 


H  (R)  = 


K  Yo  I 
i 


'£2 


-II  ^o£2  ^  ,  —II  ,  ^epT^^o) 

+  N  (R)( - — - -  )  +  M^J^{R)(_J:1 - 

‘*^2  A,,. 


a 


(6.5) 


Recalling  the  notation  that  was  introduced  in  Chapter  IV,  the 
subscript  q  is  used  to  simplify  the  group  of  subscripts  Thus  if 
it  is  necessary  to  identify  even  and  odd  functions,  the  shorthand 
notation  q  can  not  be  x’sed  and  it  is  necessary  to  use  the  more  bulky 
notation  in  equation  (6.5). 

The  surface  curren';  on  the  quarter  plane  is 

J(R)  =  n  X  [H  (r,  0=Tr,  0  _<  ({i  <_  tt)  -  H  (r,  0=n,  tt  iji  ^  2ti)]  (6,6) 

where  n  is  a  unit  normal  vector  and  n  =  y.  _ 


Performing  the  calculations,  it  is  found  that  NeJ2(R)  and 

do  not  contribute  to  the  current.  The  actual  calculation  is  very 

complicated  because  of  the  different  unit  vector  directions  on 

opposite  sides  of  the  quarter  plane,  and  because  a  point  is  described 

by  (r,  0  =  Ti,  on  one  side  of  the  quarter  plauie  and  by 

(r,  0  «  V,  2ti  -  ({i)  on  the  other  side.  The  resiilt,  however,  is 

_jj  _  _ _ 

obvious  after  one  realizes  that  Hej^2(R)  and  Mo^^(R)  ai*e  composed  of 
Lame  polynomials.  This  implies  that  these  vector  wave  functions  are 
free  space  solutions  of  the  vector  wave  equation  (no  scattering  body) 
and  are  therefore  continuous  at  the  quarter  plane.  The  contribution 
of  the  remaining  two  terms  to  the  surface  current  is 


J(R}  =  -2k  Yo  I 
I 


^012 


-  v(v+l)  hv^'(Kr) 


Kr 


,  (*)  ;  .  (r  |l  t  sing  j  i  , 

sin  <j) 


a 


e£l 


hi^^Kr)  ^eoT(<i>) 


Kr 


ai' 
sin  (|) 


(6.7) 


For  the  sake  of  notational  economy,  the  eigenvalue  v  has  not  been 
subscripted.  This  will  occur  whenever  it  is  clear  that  the  eigenvalue 
to  be  used  corresponds  to  an  identified  eigenfunction. 

Equation  (6.7)  was  derived  from  equation  (6.3a)  and  is  thus  good 


for  r  ^  Tq .  For  r  £  r^ ,  the  surface  current  is 


J(R)  =  -2k  Yc  I 


iv(icr) 


Kr 


10l» 


^£2 


(r  jv  (Kr))'  ^  sin2  ^  ^ 

+  -  ^  ^*  (^)  r  ) 

sin  4> 


Kr 


*  -^7 -  — -  r 


eu 


sin  <}) 


^  1 


(6.8) 


The  variable  <j)  in  equation  (6.7)  and  (6.8)  is  restricted  such  that 

0  <Ji  £  TT. 

Up  to  this  point,  no  approximations  have  been  made.  In  principle, 
eqixations  (6.2),  (6.3) »  (6.7)  and  (6.8)  for  E,  H,  ai).d  J  can  be  used 
to  determine  the  exact  fields  and  currents  anyvhere  due  to  a  dipole 
source  with  any  location  and  any  orientation.  Admittedly  the  equations 
are  complicated  but  it  is  practicable  to  use  them  in  calculations 
where  the  source  point  or  observation  point  is  no  more  than  one  or 
two  wavelengths  from  the  tip.  Moreover,  it  may  be  possible  to  deve3.op 
asymptotic  expressions  starting  with  these  expeinsions. 


Fields  and  Currents  Near  the  Tip  of  the  Quarter  Plane 

Consider  what  happens  as  Kr  becomes  very  smell.  The  E  field  for 
ro  >  r  is  given  by  equation  (6.2b).  The  R  dependent  terms  of  the  E 
field  are  Mq2(R)  and  Nq2^(R).  First,  consider  Mq2(R)-  From  equation 
(it. 17) 


"q2 


(6.9) 
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For  small  Kr,  the  Bessel  function  can  he  written 


Jv  (fr): 


|T  (Kr)^ 

2^+1  r(v+i.5) 


(6.10) 


where  r(x)  is  the  Gamma  function.  The  lowest  Neumann  eigenvalue  is 


V  -  0.8].lt.  For  this  eigenvalue 

Mjj2(R)  ^  0.1t28  (v  =  0.8lU) 


(6.11) 


Note  that  the  subscript  £  has  been  suppressed.  This  will  occur  in  all 
of  the  following  equations  when  it  is  clearly  understood  which 
eigenvalue  is  being  used.  The  next  highest  eigenvalue  is  v  =  1.  The 

A 

two  corresponding  Lame  functions  are  even,  and  the  vector  wave  functions 


-I  Kr  1  ®e2  (v  =  1,  P  =  0.5) 

Mg2(R)~  "Y  _ 

^  me2  (v  =  1,  p  =-0.5) 


(6.12) 


Next,  Nqi(B)  is  examined  when  xr  is  small.  From  equation  (l+.l8) 

-I  Jv  (kp)  t  .  (r  jv  (Kr))'  ^ 

=  “T;: -  ^ql  -  (r  X  mai) 


(6.13) 


For  small  Kr 


(r  jv  (Kr))' 


nr(vn)(Kr)^'-^ 

2''  (2v+l)r(v+0.5) 


(6.1U) 


The  lowest  Dirichlet  eigenvalue  is  v  =  0.296.  The  corresponding  Lame 
iijiction  is  even,  and  the  vector  wave  function  is 


0-'^T6  (Kr)"°''^°^  [£ei(v  =  0.296) 


+  1.296  (r  X  mei  (v  =  0.296))  ] 


(6.15) 
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The  next  eigenvalue  is  v  =  1.  The  corresponding  Lame  function  is  odd, 
and  the  vector  wave  function  is 

^  J  2  (r  •  inci  (v  =  1))1  (6.16) 

The  dominant  term  for  icr  very  small  is  obviously  Ne3^(R).  The 
normalization  constant  associated  with  this  term  is  (from  Table  6) 


=  2. 9*^0 


(6.17) 


Thus  the  E  field  very  close  to  the  tip  of  the  quarter  plane  can  be 
written 

E(R)  =  j  K  (0.261))(Nei  (Ra)  •  a)(Kr)'°‘'^°^ 


[  +  1.296  (r  X  mg]^)  ]  ter  <  <  1 


(6.18) 


where  it  is  understood  that  the  eigenfunctions  to  be  used  in  the  vector 
wave  functions  correspond  to  the  eigenvalue  v  =  0.296. 

Next  consider  the  H  field.  Using  the  small  argument  approximation 
for  the  Bessel  functions,  and  evaluating  the  first  few  terms,  it  is 
seen  that  the  dominant  behavior  of  the  H  field  near  the  tip  is 

H(R)  =  -K  Yo  (0.02lt)(NF^(Rc)  •  a)(icr)'°'^®^ 

C2 


[  £^2  +  1.8lU  (r  X  ra^g)  1  xr  <  <  1 


(6.19) 


The  eigenfunctions  in  this  equation  correspond  to  the  eigenvalue 
\)  =  0.8ll».  The  normalization  constant  is 


Ao2  =  12.233 


(6.20) 
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The  dominant  term  in  the  surface  current  is  the  same  as  for  H 
and  is  given  hy 


J(R)  =  -2  K  Yo  (0.031)(mJ2^Ro)  *  1) 


p 

A  1  +  sin'^  ij> 


Kr  <  <  1 


(6.21) 


Equations  (6.l8),  (6.19),  and  (6.21)  give  only  the  singular 
components  of  the  fields  and  current  near  the  tip  of  the  queirter 
plane.  These  equations  govern  the  dominant  field  and  current  behavior 
whenever  they  do  not  vanish.  For  certain  source  locations  and 

^  „  —jj 

orientations  a  •  and  a  •  will  he  zero  and  it  will  he  necessary 
to  investigate  other  terms  in  the  field  and  current  expansions.  This 
will  ho  done  later. 

Next,  the  singular  fields  are  examined  at  the  edges  and  on  the 
surface  of  the  quarter  plane.  From  equation  (6.l8),  the  E  field,  for 
a  constant  value  of  r  is 


E(R)  =  A  0g^(0)  $gi((|))  r 


a  II  +  sin*- 


sin2  0  4  sin2 


'*'el 


B 


1  +  sin^  9 


sin2  0  +  sin^  (fi 


9ei(0)  I'elU)  e 


Kr  <  <  1 


(6.22) 


where  A  and  B  are  constants.  The  edges  of  the  quarter  plane  correspond 
to  (0  =  TT,  (|)  =  O)  and  (0  =  tt,  4»  =  u).  Approaching  the  edge  on  the 
({»  =  0  sector  (see  Figure  8),  it  is  seen  that  the  r  component  goes  to 


L 


zero,  the  0  component  is  zero,  and  the  component  goes  to  infinity. 


Similar  behavior  is  found  for  the  edge  (9  =  ti,  6  =  n).  Approaching 

the  edges  from  the  quarter  plane  itself  (0  =  tt)  ,  it  is  seen  that  the 

^  ^ 

r  component  is  zero,  the  4i  component  is  zero,  and  the  0  component 
goes  to  infinity  (see  Figure  9).  This  edge  behavior  is  consistent 
with  the  results  from  half-plane  theory.  On  the  quarter  plane 
itself,  the  parallel  (r,  ({>)  components  are  zero,  and  the  perpendicular 

A 

(0)  component  is  finite  except  at  the  edge. 

From  equation  (6.19)  the  K  field  is 

H(R)  =  C  ©02(6)  ^ 

D  1  +  sin^ 

+  - -  ©02^®) 

sin^  0  +  sin^  <J>' 


D 


1  +  sin^  e 


sin^  6  +  sin^ 


002(e)  '^o2('^)  ® 


Kr  <  <  1 


(6.23) 


where  C  and  D  are  constants.  Approaching  the  (0  =  ti,  ij)  =  O)  edge 
from  the  (|)  =  0  sector,  it  is  seen  that  the  r  component  is  zero,  the 

A  ^ 

4)  component  goes  to  infinity,  and  the  0  component  is  zero.  Similar 
behavior  is  found  for  the  other  edge.  Approaching  the  edges  from  the 
quarter  plane  itself  (0  =  v),  it  is  seen  that  the  r  component  goes  to 
zero,  the  component  goes  to  infinity,  and  the  0  component  is  zero. 
This  edge  behavior  for  the  H  field  is  also  consistent  w;th  half-plane 
theory.  On  the  quarter  plane  itself,  the  parallel  components  (r,  ^) 


ilO 

are  finite  except  at  the  edges,  where  the  r  component  is  zero,  and 
the  (f  component  is  infinite.  The  perpendicular  component  (6)  is  zero 
It  was  mentioned  in  Chapters  IV  and  V  that  the  electromagnetic 
fields  must  satisfy  edge  and  tip  conditions.  These  conditions  are 
necessary  for  uniqueness,  and  they  do  not  occur  in  the  case  of 
scattering  bodies  with  well  defined  surface  normals.  A  good  summary 
of  this  subject  can  be  found  in  Jones  [^,19].  Jones  describes  the 
conditions  imposed  by  several  authors  and  concludes  that  although 
their  viewpoints  are  different,  their  results  are  in  agreement.  For 
the  quarter  plane  problem,  perhaps  the  simplest  way  of  stating  the 
condition  is  to  require  that  no  energy  be  radiated  from  the  tip  or 
edges.  Thus  the  condition  can  be  reduced  to  satisfying 


Re  f  E  X  1*  •  ds  =  0 
s 

on  a  surface  enclosing  the  tip  (or  edge)  as  this  surface  shrinks  to 
the  tin  (or  edge).  This  can  be  done  by  inspecting  the  order  of  the 
singularities  of  the  proper  field  components  as  the  edges  and  tip 
are  approached.  The  tip  and  edge  conditions  are  satisfied  by  the 
solutions  given  here. 

From  equation  (6.21),  the  surface  current  on  the  quarter  plane 


IS 


J(R)  =  E  <tQ2((^)  <^  +  F 


1  +  sin‘= 


sin  (|> 


<{j^2((}>)  r  Kr  <  <  1 


(6. 


Figures  10,  11,  and  12  are  sketches  of  the  r  and  (j)  components  of  the 
surface  current,  and  the  total  surface  current  flow  near  the  tip  of 


2h) 


the  quarter  plane. 


Fig.  10.  Variation  of  J^.  for  >:r<<l 
(From  equation  6.2k) 


Fip:.  11.  Variation  of  J.  fo 
(From  equation  6.2U) 
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All  of  the  fields  end  currents  just  described  are  the  dominant 
terms  whenever  they  are  not  zero.  In  order  to  determine  when  these 
terms  vanish,  it  is  necessary  to  examine  the  source  terms.  The  source 
term  for  the  E  field  is 

.(2). 


jl  «  0.381t  ho,296 

N  (Ro)  •  a  =  - 


+  (r»  ho^296  (»ro))’ 

Kr^  sin^  Go  +  sin2  ({iq 


(  |l  +  sin^  (j)o  0ei(®o)  ♦eiUo)  *  a 
+  1 1  +  sin2  Go  0ei(®o)  *ei(^o)  §0  ’  a  )  (6.2?) 


First,  consider  a  source  in  the  direction,  i.e.,  a  =  Xg  .  Then 

(2) 

_II  _  />  (0.304)  hjQ  296  'J’ei(^o) 

Nei(Ro)  •  a  =  - - - 

<Tn 


(6.26) 


In  general,  this  term  is  not  zero  for  any  Go  or  410  except  of  course 

A 

Gq  =  IT.  Next  consider  a  source  in  the  ({>0  direction.  Then 


-II,-,  ^  (ro  ho.296 

Nei(Ro)  •  a  = - ^ - 

Krc 


1  +  sin2  0^1(60)  <I'el('J'o^ 

sin2  Go  +  sin^  ({i^ 

(6.2?) 


The  eigenfunction  'tei(*{>o)  is  oi*  form 

CO 

<*’el(<f'o)  =  I  B2ni  cos  2m  <j)jj 
ra=0 


so  that 


00 


111* 


^el^^o)  =  -  I  2m  B2ni  sin  2m  (j>o 

m=l 


Thus  this  terni  is  zero  for  <t)o  =  Tt/2,  3ti/2,  which  corresponds  to  the 

A  A 

z  =  0  plane.  The  vector  a  =  is  perpendicular  to  this  plane.  These 
results  imply  that  equation  (6.l8)  for  the  E  field  close  to  the  tip  of 
the  quarter  plane  is  not  valid  for  this  particular  source  location  and 
orientation.  Equation  (6.27)  is  also  zero  for  (J)^  =  0,  tt.  A  source  in 

A 

the  (}io  direction  at  (J>o  “  0,  ir  is  a  source  in  the  plane  of  the  quarter 
plane  (y  =  0  plane)  but  not  on  the  quarter  plane  itself,  and  oriented 
perpendicular  to  it.  It  is  not  surprising  that  the  current  vanishes 
for  this  case,  since  this  type  of  source  produces  no  scattered  field. 
It  is  also  zero  for  9©  =  if.  Next  consider  a  source  in  the  Qq 
direction.  Then 


N;i(Ro) 


-  (3^0  ho?J96  (f'l'o))*  1  +  sin^  Oc  0el(eo)  ^el^'f’o- 

a  =  — - -  -  ■  - - - 


KTn 


3in2  So  +  sin2  (j)^ 


(6. 


This  term  is  zero  for  9©  =  0.  Again  this  corresponds  to  a  so\u:ce 
located  in  the  plane  of  the  quarter  plane  and  oriented  perpendicul'or 
to  it;  on  the  other  hand,  equation  (6.28)  is  not  zero  for  a  source  on 
the  quarter  plane  itself. 

Since  equation  (6.I8)  does  not  describe  the  dominant  E  field  when 
the  source  is  in  the  z  =  0  plane  and  oriented  perpendicular  to  it .  it 
is  necessary  to  look  for  the  next  most  dominant  term.  This  term 
corresponds  to  NgQ^(R)  with  v  -  1.130.  The  E  field  is 


13.5 


E(R)=jic(0.011)  (Nei(Ro)  *  a)(icr)°*^^  [  £el  +  2.13  (r  x  rngl)] 

xr  <  <  1  (6.29) 

The  normalization  constant  is 

Aei  =  llt.201  (6.30) 

Note  that  this  field  is  not  singuleir  at  the  tip.  The  edge  "behavior, 
however,  is  the  same  as  before. 


The  source  term  for  the  dominant  H  field  and  current  is 


-  1  +  sin2  ().o  ©02(90)  90  •  a  )  (6.31) 


This  term  is,  of  course,  zero  for  a  source  in  the  r^  direction.  For 

A 

a  source  in  the  4*0  direction, 

M^2(Ro)  •  a  =  (<ro)  — ^ - f I-?-..  P.2  -  ^  .  (6.32) 


sin2  00  +  sin2  (J)^ 


This  term  is  not  zero,  in  general,  for  any  source  location  except 
4>o  =  0,  which  has  already  been  discussed,  and  of  course  for  60=1?. 

A 

For  a  source  in  the  6^  direction, 

.  ;  =  .  *«  ®o2(8.) 

I  sin2  Oq  +  sin2  (Jiq 


(6.33) 
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This  term  is  zero  for  ‘J'o  =  ti/2j  Thus  equations  (6.19)  and  (6.21) 

do  not  give  the  dominant  "behavior  of  the  H  field  and  current  when  the 
source  is  in  the  z  =  0  plane  and  in  the  6©  direction.  Since  the  source 
term  is  also  zero  for  a  source  in  the  r^  direction,  it  is  zero  for  any 
source  in  the  z  =  0  plane  and  oriented  psirallel  to  it.  This  is 
reasonable  when  one  considers  the  symmetrical  nature  of  the  source  and 
the  antisymmetrical  nature  of  the  current  shown  in  Figure  12.  The 
dominant  current  behavior  for  this  type  of  source  corresponds  to  the 
eigenvalue  v  =  0.296.  The  cvirrent  is 

2kYo  (0.172)  (Ne?(Ro)  •  a)(Kr)°'^^^  <?eiU)  r 
j(R)  = - f - 

sin  i|) 


<r  <  <  1  (6.3*t) 


The  corresponding  H  field  is 


H(R)  =  -  K  Yq  (0.261»)(NeJ(Ro)  *  a)(Kr)°'^^^ 


KT  <  <  1  (6.35) 


These  terms  are  not  singular  at  the  origin.  The  edge  behavior  of  the 
H  field  is  the  same  as  before,  however,  and  is  consistent  with  half¬ 
plane  theory.  The  surface  current  is  in  the  r  direction  only,  and 
goes  to  zero  at  the  tip.  Its  behavior  for  Kr  <  <  1  is  sketched  in 
Figure  13,  and  the  current  flow  is  sketched  in  Figure  lit. 

The  dominant  fields  and  currents  near  the  tip  of  the  quarter 
plane  for  the  various  source  locations  are  summarized  on  pages  118, 
119,  and  120. 


Summary:  Example  1 


Fig.  15,  Source  in  the  z=0  plane  and 
perpendicular  to  that  plane 


0  X3  ^  *** 

E(R)=jK(0.011)  (Nei(Rt>)  *  a)(Kr)  '  [lel  +  2.13  (r  x  mel)] 

V  =  1.13  ,  Kr  <  <  1 


(6.29 


H(R)  =  -K  Yo  (0,02U)(M^2{Ro)  '  a)(icrr°’^®^ 


(  1^2  +  1.8lU  (r  X  5  ^  -  0.6l4  ,  xr  <  <  1  (6.19 


J(R)  -2  K  Yo  (0.031) (M^gd^o)  *  a)(Kr)”®’^®^ 


1  +  sin^  •}> 


*  +  1-228 


(6.21 


V  =  0.8lli 


Kr  <  <  1 


SuBMTx:  Exciq>le  2 


E(R) 


•  H(R) 


y 


Fig.  16.  Source  In  the  z»0  plane  and 
parallel  to  that  plane 


«  J  K  (0.2610(KejL 

[  4-  1.296  (r  X  Bel)  1  ''  '  0*296  ,  icr  <  <  1 


B  -  K  Yq  (0.26U)(Ne^{Ro)  •  a)(»cr)°’^^^ 


Summary:  Example  3 


E(R)  « 


H(R) 


J(R)  »  -2 


y 


Fig.  17.  Source  location  and  orientation 
arbitrary  to  the  extent  that  it  is  not 
in  the  z=0  plane 


j  K  (0.264)(Ne][  (Ro)  •  a)(K:r) 


[  +  1.296  (r  X  ]  V  s  0.296  ,  xr  <  <  1 


-K  Ye  (0.024)(M^^(R,)  . 


[  +  1.8llt  (r  X  m^g)  ]  V  =  0.814  ,  icr  <  <  1 


Yo  (0.03l)(Mi[iJ(Ro)  •  a){Krr°‘-"®^ 


02 


1  -t-  sln^  ^  I  /  V  * 


♦c2<«)  ^  +  i-228  ^ 


V  »  0.8l4  ,  Kr  <  <  1 
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Fields  and  Currents  Near  the  Tip  of  Angular  Sectors 
vith  =  0.1  and  =  0.9  for  a  Unit  Dipole  Source 

The  dominant  behavior  of  the  E  field  in  the  vicinity  of  the  tip 

is  determined  by  the  lowest  even  Birichlet  eigenvalue,  and  the 

dominant  behavior  of  the  H  field  and  current  is  governed  by  the  lowest 

odd  Neumann  eigenvalue.  In  Chapter  III,  these  eigenvalues  were 
2  2 

determined  for  k  =  0.1  and  k  =  0.9. 

When  the  parameter  k^  -  0.1,  the  corner  angle  of  the  plane  angulsn: 
sector  is  0.795if,  thus  as  k^  becomes  smaller  the  E  field  behavior 
should  be  similar  to  that  for  the  half  plane.  For  =  0.1, 

Vei  =  O.UOTj  and  it  is  evident  that  (ier)'*®^“^,  corresponding  to 
equation  (6.I8),  is  approaching  (icr)”®*^  as  the  plane  angular  sector 
approaches  a  half  plane.  In  order  to  compare  this  singular  part  of 
the  E  field  with  half-plane  theory,  consider  the  2*0  ((J)  »  it/2,  3it/2) 

A 

plane.  In  this  plane  the  unit  vector  r  corresponds  to  the  radia!l 

A  A 

vector  D  in  cylindrical  coordinates.  The  unit  vector  ((>  in  the  sphero- 

A 

conal  system  is  the  same  as  -z  in  the  cylindrical  system,  and  the 

A  A 

sphero-conal  unit  vector  0  approximates  the  cylindrical  unit  vector  (f. 
(See  Figures  I8  and  19.)  Only  the  r=p  and  0  components  of  equation  (6.18) 
exist  in  this  plane.  Thus  in  the  cylindrical  coordinate  system,  the 
singular  fields  are  the  p  and  <{>  components.  They  vary  as  (xp)"^^^, 
which  84?rees  with  half-plane  theory.  For  a  general  discussion  of  field 
behavior  at  an  edge,  sea  Jones [19]. 

The  dominant  current  behaves  like  equation  (6.21).  The  eigenvalue 
ccrresponding  to  this  equation  for  k^  =  o.l  is  =  0.6l3.  This 
eigenvalue  is  also  approaching  0.5  as  k^  approaches  zero.  Thus  the  r 
dependence  of  the  dominant  current  is  also  approaching  (t<r)“^^^. 


Fig.  18.  A  plane  angular  sector  with  a  large 
corner  angle  and  the  unit  vectors  in  the 
z=0  plane 


Fig.  19.  A  half  plane. in  a  cylindrical 
coordinate  system  and  the  unit  vectors 
in  the  z=0  plane 
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A  sketch  of  this  current  is  given  in  Figure  20.  The  current  is 
predominantly  pavr-allel  to  the  edge  and  varying  approximately  as 
which  is  in  agreement  with  half -plane  theory. 

The  above  discussion  is  rather  qualitative  and  hy  no  means 
describes  the  exact  behavior  of  the  fields  and  current  as  the  plane 
angular  sector  becomes  a  half  plane.  The  purpose  of  the  discussion 
is  only  to  indicate  that  as  -»■  0,  the  various  field  and  current 
components  do  appear  to  be  approaching  the  correct  half -plane  theory, 
and  that  the  half-plane  problem  could  be  studied  if  desired. 

The  lowest  even  Dirichlet  eigenvalue  and  the  lowest  odd  Hevimann 
eigenvalue  wei'e  also  determined  for  k^  =  0.9.  This  corresponds  to  a 
plane  angular  sector  with  a  comer  angle  of  0.205ir.  This  case  will 
not  be  considered  in  detail  but  will  be  used  only  to  indicate  the 
variation  with  r  of  the  fields  and  currents  near  the  tip  of  a  sharp 
sector.  The  eigenvalue  Ve^  =  0.171,  so  that  the  dominant  E  field, 
which  varies  as  is  approaching  (xr)"^.  The  eigenvalue 

v<32  =  0.970,  so  that  the  dominant  H  and  J,  vhich  vary  as 
are  approaching  (<r)^.  These  are  only  limiting  values,  because  the 
vector  wave  functions  corresponding  to  these  eigenvalues  go  to  zero 
for  k*^  =  1.  The  reason  for  this  is  that  the  scattering  body  is 
becoming  smaller  as  k^  approaches  1,  and  for  k^  =  1  there  is  no 
scattering  body.  The  singular  terms  do  not  exist  for  this  case. 

The  behavior  of  the  dominant  fields  and  currents  is  tabulated  in 
Table  7  for  k^  =  0.1,  0.5,  and  0.9. 


The  S  Field  Near  the  Tlij  of  the  Quarter  Plane 


for  a  Plane  wave  Incident 

In  order  to  compare  the  present  solution  with  eai*lier  work,  the 
illumination  of  the  quarter  plane  by  a  plane  wave  will  be  treated. 
Consider  what  happens  as  the  source  is  moved  far  from  the  tip,  i.e., 
as  Kr^  -i-  For  this  case,  the  lai*ge  argument  approximations  can  be 
used  for  the  Hankel  functions.  In  example  1  on  page  118,  the 
source  term  for  the  E  field  is 


I  V  1 

— TT —  (r©  hj^  3^2  0ei(9o)  ^el^'i’o) 


|1  +  sin2  00 

Correct  to  order  l/xro ,  the  r©  dependence  is 


(^0  ^1^13 e“«5''^o 


and  the  E  field  due  to  this  soui-ce  can  be  written 

m  =  j(0.107)eJ“-5«5»  Ei  °.n(9.) 

|l  +  sin2  00 

[  ^01  2.13  (r  X  ni0i )  ]  *1*0  ~  1^/2,  3'Jr/2, 


V  =  1.13,  KU  <  <  1,  Kro>  >  1 


where 


is  the  incident  field  at  the  tip.  E^  also  can  be  considered  as  the 
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field  of  a  pliane  wave  at  the  tip.  In  particular,  consider  a  plane 

A  A 

wave  propagating  in  the  -y  direction  with  the  E  field  in  the  z 
direction.  The  total  field  in  the  vicinity  of  the  tip  is 

E(R)  =  j  (0.096)  Ei  (icr)°‘^^  (le^^  +  2.13  (r  x  %!)) 

V  =  1.13,  Kr  <  <  1,  icr^  >  >  1  (6.1t0) 

If  the  plane  wave  is  propagating  in  the  -x  direction  with  the  sajce 
polarization,  the  E  field  is 

E(R)  =  j  (0.078)  Ei  e’5°-5^^’^(Kr)°‘^^  (Ig^^  +  2.13  (r  x  5^^^)) 

V  =  1.13,  xr  <  <  1,  Kro>  >  1  (6.1tl) 

Next  consider  example  2.  For  a  plane  wave  propagating  in  the  -y 
direction  with  x  polarization,  the  total  E  field  in  the  vicinity  of 
the  tip  is 

E(R)  =  j (1-065)  E^  (Ig^  +  1.296  (r  x  inei)) 

V  =  0.296,  xr  <  <  1,  xr^  >  >  1  (6.U2) 

Equations  (6.38)-(6.U2)  ideally  yield  no  new  information  concerning 
the  behavior  of  the  fields  very  close  to  the  tip.  For  the  dominant 
field  terms,  the  range  of  the  source  for  any  of  the  three  cases 
considered  on  pages  II8,  119 , and  120  affects  o^^ly  amplitude  and 
not  the  form  of  the  fields.  However,  several  authors  have  studied  the 
behavior  of  the  scedar  fields  near  the  tip  for  scalar  plane  v/ave 
illumination.  Braunbeck[5]  considered  the  form  of  the  eigenfunction 


solution  of  L&place's  equation.  Using  physical  reasoning;,  he 
douennined  that  the  solution  varies  as  r“  (a  <  0.5).  Radlov(3] 
used  an  integral  transform  technique  to  solve  the  scsJ.'ar  problem.  He 
interprets  his  result  as  being  the  scattered  E  field  due  to  an  incident 
plane  wave  with  the  E  field  parallel  to  the  quarter  plane.  He  estimates 
that  the  E- field  varies  roughly  as  r®*^5.  Radlow  also  mentions 

0  3 

another  calculation  in  which  Noble  estimates  the  variation  to  be  r  . 
The  eigenvalue  in  the  present  paper  which  corresponds  to  these 
exponents  is  Vei  =  0.29^^.  Radlow' s  treatment  of  the  scalar  problem 
appeal's  to  be  correct  and  the  results  ha  obtains  are  a  valuable 
contribution  to  this  problem;  however  his  interpretation  of  his 
solution  as  one  component  of  the  electric  field  is  clearly  in  eiror. 

The  vector  solution  given  here  shows  that  the  scattered  field  depends 
on  the  polarization  of  the  incident  plane  wave  and  t)iat  it  varies 
as  r''"*. 

Far  Zone  Fields  and  Currents  Due  to  a  Unit  Dipole 
Source  Close  to  the  Tip  of  the  Quarter  Plaa~.c 

Next  the  reciprocal  case  is  studied;  i.e.,  the  far  zone  fields 
are  determined  for  the  source  located  cD.ose  to  the  tip  of  the  quarter 
plane.  Perhaps  the  simplest  way  to  determine  these  fields  is  to  start 
with  equation  (6.2a).  Let  r^  be  very  small  and  use  the  small  argiiment 
approximation  for  the  Bessel  functions  with  sirgument  icr^ .  Then  let  r 
be  very  large  and  use  the  large  argument  approximation  for  the  Bessel 
functions  with  argument  kt.  The  procedure  is  essentially  the  same  as 
before,  except  that  r  and  r^  are  interchanged. 
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Equation  (6.2a)  is 

_  _TT  _  (Re)  '  a 

E(R)  =  j  K  J  Mq2  (R)  (  - : -  ) 

q  ^q2 


-II  -  '*01  ^ 

tNji(R)(-^^ - )  r>r, 

Aqi  j 


(6.U3) 


The  source  terms  are  examined  for  r^  very  small,  and  it  is  seen  that 
the  dominant  term  is 

~I  \  .  (ro  Jv(Kro))‘ 

(Ro)  *  a  =  -  *  a  +  - - = - 


(roX  mei)  •  a 


{6M) 


with  V  =  0.296.  Using  the  small  argument  approximation  for  the 
Bessel  fxmctions,  this  becomes 


N?  (Ro)  •  a  1.236  (  tg,  ‘at  1.296 


(r^x  mg]^)  •  a  ) 


(6.1*5) 


The  R  dependent  term  is 


(R)  = 


hi^^(Kr)  —  ,  (rhv^^(Krj)' 


^ei 


(r  X  mei) 


(6.U6) 


Using  the  large  argument  approximations  for  the  Hankel  functions  and 
neglecting  terms  of  order  higher  than  l/Kr,  this  becomes 

^  jO. 11*8*1  ^ 

No,  (R)'^  e""  -  (rxmei) 


(6.1*7) 
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Nov  consider  hhe  source  to  be  located  in  the  z  =  0  plane  and  in 

A 

the  -6o  direction.  This  corresponds  to  example  2  on  page  119*  In 
particulars,  let  the  source  be  located  on  the  line  6©  =  if/2, 

(})c  -  Tr/2,  (-00  =  x).  The  E  field  is 

E(R)'^  j  2.157  e’^°*^^®^(icro)**°'^°**  -  (r  x  mei)  (6.1*8) 

UliT 


In  the  2=0  plane  this  is 

E(R) -- J2.301*  e‘^°‘^^®’'(icro)~°‘'^°**  ®  (6.1*9) 


This  pattern  is  normalized  to  unity  at  0  =  n  and  plotted  in  Figure  21. 
C.  T.  Tai(20]  has  computed  the  pattern  in  this  plane  for  the  same 
source  but  with  a  half  plane  scattering  body.  The  two  patterns  are 
essentially  the  same.  As  a  matter  of  fact,  the  pattern  has  the  same 
general  shape  for  any  sector  angle  provided  Kro  <  <  1* 

In  the  X  =  0  plane,  the  E  field  is 


^  j0.ll*8Tf,  .-O.JOl* 

E(R)'^  J  1.666  e'^  (xro) 


l*iTr 


-  ^  *  (0.51*0) 


'J>ei(4>)  6 
!  1  +  sin^  ij) 


(6.50) 


Both  components  of  this  pattern  are  plotted  separately  in  Figures  22 
and  23.  These  patterns  are  normalized  relative  to  Figure  21.  Note 
that  the  magnitude  of  E(}i  in  Figure  22  is  quite  small  ^  in  fact  the 
peak  amplitude  is  more  than  20db  down  from  the  peak  amplitude  of  the 
field  in  the  z  =  0  plane. 
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The  y  =  0  plane  is  composed  of  four  sectors,  $  =  n ,  6 
ij)  =  0,2n,  and  0  =  ti.  For  ({>  =  0,  the  E  field  is 

-jicr 


=  0, 


E(R)-  -  j(2.157)  ^ 


itirr 


1.  ,t..sing..9_  Q.  (g)  g 
sin  0  ^ 


For  0=0,  the  S  field  is 


E(R)^-  J(2.157)  ~ 


-jKr 


Uiir 


1  +  sin^  <{>  /.X  2 


sin  4)  I 


(6.51) 


(6.52) 


For  4>  =  IT,  the  E  field  is  given  by  equation  (6.5l)»  The  fieids  for 
<|)  =  0,  0  =  0,  and  <))  =  it  are  all  parallel  to  the  y  =  0  plane.  For 
0  =  IT,  the  parallel  field  is  of  course  zero.  A  field  perpendicular 
to  the  quarter  plane  does  exist,  however.  It  is 


E(R) 


j(l.U06)(Kro ) 


-0.70J» 


JO-IUBt: 

e 


l4Tir 


4'el(4')  0 
Isin  4)1 


(6.53) 


This  field  exists  on  both  sides  of  the  quarter  plane,  the  E  vectors 
po5nting  either  into  or  out  of  the  sector  on  both  sides.  The  patterns 
given  by  equations  (6.5l),  (6.52),  and  (6.53)  are  normalized  relative 
to  Figure  21  and  plotted  in  Figure  2U.  Note  that  the  edges  of  the 
quarter  plane  appear  to  have  a  guiding  effect  on  the  field.  This 
might  have  some  useful  applications. 


Fig.  21.  Far  field  L'l  the  z  =  0  plane  due  to  an  infinitesimal  dipole  source 

cloiie  to  the  tip  of  the  quarter  plane 


22.  <J>  component  of  the  far  field  in  the  x  =  0  plane  due  to  an  infinitesimal  dipole  source  close 

to  the  tip  of  the  quarter  plane  (source  in  the  z  =  0  plane  and  parallel  to  it) 


13"? 


•CH 

r- 


23.  e  con^onent  of  the  far  field  in  the  x  =  0  plane  due  to  an  infinitesimal  dipole  source 
close  to  the  tip  of  the  quarter  plane  (source  in  the  2=0  plane  and  parallel  to  it) 
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The  dominant  far  zone  current  is 


J(R)'-  -2j  Yo  (1.1*06)  e 


JO.IWti,  ,-0.70U 


i^itr 


^'ei('}') 
sin  <|i 


(6.5I;) 


At  a  constant  r,  the  behavior  of  J(R)  is  the  same  as  that  shown  in 
Figure  13.  Equation  (6.3*<)  is  valid  for  the  current  near  the  tip 
due  to  this  source,  so  Figures  13  and  lU  show  the  behavior  of  the 
current  there.  Notice  that  the  current  in  the  far  zone  is  decaying 
as  1/r  and  the  current  near  the  tip  is  varying  as  rO‘296. 
currents  have  the  same  behavior  with  respect  to  (fi. 

Next,  consider  the  source  location  to  be  on  the  quarter  plane, 
still  in  the  z  =  0  plane  and  still  in  the  -8o(y)  direction.  This 
source  corresponds  to  a  vertical  dipole  on  the  quarter  plane  and 
close  to  the  tip.  The  only  thing  that  is  changed  is  that  So  =  it 
instead  of  ir/2.  Equations  (6.1»8)-(6.5^)  are  unchanged  except  in 
magnitude.  Multiply  each  of  these  equations  by  (2.21)  to  obtain  the 
far  fields  and  current  for  this  source,  figures  21»2U  give  the 
patterns.  Actually,  these  figures  give  the  patterns  for  any  source 
located  in  the  z  =  0  plane  and  peirallel  to  it,  and  very  close  to  the 
tip  of  the  quarter  plane.  Of  course,  the  amplitude  of  the  field  is 
zero  for  the  obvious  case  of  a  source  on  the  quarter  plane  and 
parallel  to  it.  The  amplitude  is  also  zero  for  the  source  on  the  x 
axis  and  oriented  in  the  y  direction.  This  is  the  configuration  that 
produces  no  scattered  field.  In  Figures  21  to  24  it  is  interesting 
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to  note  that  the  magnitude  of  E  increases  as  the  electric  current 
dipole  approaches  the  tip.  The  domiueuit  term  of  the  electric  field 
for  these  cases  is  entirely  due  to  currents  induced  on  the  quarter 
plane.  Thus  E  can  he  interpreted  as  a  field  scattered  from  the 
quarter  plane  due  to  a  current  dipole  in  close  proximity  to  the  tip. 

Now  consider  the  source  location  shown  in  example  1  on  page  ll8. 
Tiet  the  source  he  in  the  z  =  0  plane  and  perpendiculeir  to  it.  In 
particular  consider  a  source  in  the  direction  and  on  the  (po  =  ii/2, 

A  A 

9o  =  Tr/2  line  =  z).  The  dominant  term  for  the  E  field  corresponds 
to  the  eigenvalue  Vgi  =  1.13.  The  far  zone  E  field  is 


E(i)-  J(0.199)  (r  x  mei)  (6.55) 


Uirr 


In  the  z  =  0  plane  this  is 


E(R)-  J(0.291) 

I  5 


(6.56) 


This  pattern  is  nor.malized  to  Tinity  and  plotted  in  Figure  25.  The 
pattern  in  Figure  25  also  agrees  with  Tai's  half  plane  pattern, 
which  is  to  he  expected,  because  the  genereil  pattern  shape  is  good 
for  all  sector  angles. 

In  the  x  =  0  plane,  the  far  zone  field  is 


E(R)  '-j(0.2l»3) 


l+7rr 


0.332  e 

1  +  sin^  <t> 


(6.57) 


Both  components  of  this  pattern  are  normejLized  relative  to  Figure  25 
and  plotted  separately  in  Figures  26  and  27-  The  quarter  plane 
appears  to  have  very  little  effect  on  Fi/jure  26.  The  pattern  in 
Figure  27  is  of  a  cross-polarized  field  and  is  dovn  about  lOdb  from 
the  normally  polarized  field. 

In  the  y  =  0  plane,  the  far  field  for  the  ij)  =  0  sector  is 

E(R)^j (0.199) 

4itr 

.iJLiiHii.  ©^^(e)  e  (6. 

sin  6 


For  the  0-0  sector,  it  is 


E(R)'^  J(0.199)  eJ°-5657r(^j,j0.13 

kvr 


1  sin^  i{) 
I  sin  (|)| 


(6. 


For  =  ir,  the  field  is  given  by  equation  (6.58).  These  fields 
behave  in  the  same  way  as  in  the  previous  example,  i.e.,  they  are  all 
parallel  to  the  y  =  0  plane.  For  0  =  tt,  there  is  a  perpendicular 
field  given  by 


E(R)^''  -j(0.195) 


^J0.565.(^^^)0.13 


A 

'i’el(<t>)  9 
jsin  (ti| 


(6. 


Again,  this  field  exists  on  both  sides  of  the  quarter  plane,  the 
E  vectors  pointing  either  into  or  out  of  the  sector  on  both  sides. 
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The  patterns  given  by  equations  (6.58),  (6.59) »  and  (6.60)  are 
normeilized  relative  to  Figure  25  and  plotted  in  Figure  28.  Again 
note  the  guiding  effect  of  the  edges  of  the  quarter  plane. 

The  current  is 


J(R)/^  2  j  Yq  (0.195)  e'^°‘^^^’^(icrc)°‘^^ 


Ifirr 


‘t*el(4»)  r 
sin  (j) 


(6.61) 


At  a  constant  r,  the  behavior  of  J(R)  is  similar  to  that  shovm  in 
Figure  10.  Near  the  tip,  equation  (6. 2k)  is  valid  for  this  source, 
so  Figures  10,  11,  and  12  show  the  behavior  of  the  current  there. 

The  current  in  the  far  zone  is  decaying  as  l/r;  the  current  near  the 
tip  is  singxilar,  varj'ing  as  (tcr)”^*^®^. 

One  other  source  location  is  considered  corresponding  to  this 
case.  Let  this  source  be  located  on  the  x  axis  and  oriented  in  the 
z  direction.  Equations (6. 55)-(6.6l)  have  the  same  form,  only  their 
magnitude  being  changed.  Multiply  these  equations  by  (0,8l9)  to 
obtain  tne  far  fields  and  currents  for  this  source.  Figures  25-^58 
give  the  patterns  for  this  source  and,  in  fact,  for  any  source  in  the 
z  “  0  plane  and  perpendicular  to  it,  and  very  close  tc  the  tip  of  the 
quarter  plane.  Again  it  is  mentioned  that  the  amplitude  0“  the  field 
is  zero  for  the  obvious  case  of  the  source  on  the  quarter  plane.  In 
Figures  25-28  it  is  interesting  to  note  that  the  magnitude  of  E 
decreases  as  the  electric  current  dipole  approaches  the  tip.  The 
dominant  term  of  the  electric  field  for  these  cases  again  is  entirely 


component  of  the  far  field  in  the  x  =  0  plane  due  to  an  infinitesimal  dipole  source 
the  tip  of  the  quarter  plane  (source  in  the  z  =  0  plane  and  perpendicular  to  it) 


.  ?7.  component  of  the  far  field  in  the  x  =  0  plane  due  to  an  infinitesimal  dipols  source 
close  to  the  tip  of  the  quarter  plane  (source  in  the  z  =  0  plane  and  perpendicular  to  it) 


fields  in  the  y  ■=  0  plane  due  to  an  infinitesimal 


to  currents  induced  on  the  quarter  plane,  and  E  can  be  interpreted 
a  scattered  field  for  these  cases  too. 


CHAPTER  VII 


CONCLUDING  REMARKS 

In  this  reporta  rigorous,  general  solution  for  the  diffraction  of 
an  electro3iagnetic  wave  by  a  perfectly-conducting  plane  angular  sector 
is  derived  in  the  form  of  a  dyadic  Green’s  function.  Thus,  for  example, 
the  excitation  of  the  plane  angular  sector  by  neighboring  antennas  or  by 
slots  in  the  sector  can  be  determined.  This  problem  has  not  been  solved 
previously  and  although  formal  solutions  to  the  corresponding  scalar 
problem  exist,  no  numerical  results  have  been  presented. 

The  dyadic  Green's  function  is  expressed  in  terms  of  a  complete 
set  of  vector  wave  functions  of  the  sphero-conal  coordinate  system.  The 
vector  wave  functions  contain  Lame'  functions,  which  have  been  studied 
numericariy  in  the  case  of  the  quai'ter  plane;  the  results  are  presented 
in  tables  of  the  Lame'  functions  and  their  associated  eigenvalues. 

Except  for  the  Lame'  polynomials ,  the  Lame'  functions  and  bheir  eigen¬ 
values  have  not  been  tabulated  previously. 

Numerical  results  for  the  behavior  of  the  fields  and  currents  near 
the  tip  of  the  quarter  plane  are  obtained  and  compared  with  estimates 
and  conjectures  made  by  other  authors - 

The  patterns  of  electric  current  dipoles  very  close  to  the  tip  are 
presented.  In  the  perpendicular  plane  of  symmetry,  the  pattern  is  the 
same  as  that  of  the  half  plane  similarly  excited.  The  patterns  in  the 
vicinity  of  the  edges  of  the  quarter  plane  indicate  that  a  substantial 
amount  of  energy  is  guided  away  from  the  source  along  these  edges.  In 
general,  a  dipole  very  close  to  the  quarter  plane  tip  excites  a  strong 
surface  current  flow  along  the  edges  and  near  the  tip;  consequently, 
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the  pattern  of  the  electric  current  dipole  is  severely  altered  by  the 
presence  of  the  quarter  plane. 

The  examples  that  have  bean  given  are -actually  quite  simple.  They 
are  restricted  such  tl'iat  either  the  field  point  or  the  source  point  is 
very  close  to  the  tip  of  the  quarter  plane;  in  fact,  so  close  that  only 
the  first  or  second  term  in  the  expansion  is  used.  From  the  tables  in 
Chapter  III  it  is  seen  that  approximately  200  terms  of  the  vector  wave 
function  expansion  have  been  ccanputed.  Using  these  results  it  should 
be  possible-  to  calculate  the  fields  of  sources  up  to  one  wavelength 
from  the  tip  with  good  accuracy.  At  this  distance  it  may  be  possible 
to  deteimine  the  diffraction  coefficient  approximately  and  thereby 
extend  the  solution  by  means-  of  the  geometrical  theory  of  diffraction 
so  that  the  fields  of  sources  remote  from  the  tip  can  be  calculated. 

All  of  this,  however,  must  await  future  work. 


4 


^  i- 


l  : 


I  i 


APPENDIX  A 

THE  SELF-ADJOINT  AND  POSITH/E  DEFINITE 
PROPERTY  OF  THE  TWO-DIMENSIONAL 
STURM-LIOUVILLE  OPERATOR 


The  differential  equation  of  interest  can  he  written 


Ly-Xpy=0 


where  the  Sturm-Liouville  operator  is 


L  = 


1  -  cos^  (j) 


8  /  ,  ,2  2  «  3  \ 

:  —  {  1  -  k  cos^  6  —  ) 


de 


-■  . . —  —  (  |l  -  k*^  cos^  (|)  ~  ) 

1  -  k2  cos2  0  •}> 


and  the  weight  function  is 


k^  sin^  6  +  k'^  sin^  tf> 


-  k^  cos2  9  1  -  k’^  cos^  if 


Let  y  he  an  eigenfunction 


y  =  0(e)  4>((}i) 


corresponding  to  an  eigenvalue 


X  =  v(v+l) 


(Al) 


(A2) 


(A3) 


(AM 


(A5) 
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Equation  (Al)  is  the  same  as  equation  (3.5)-  The  weight  function  is 
the  same  p  defined  by  equation  (5.29). 

The  self-adjoint  property  is  defined  by 

<w,  Ly>-<y,  Lw>ssO  (A6) 

where  <  >  denotes  a  scalar  product,  and  y  and  w  are  eigenftinctions . 


<  w,  L  y  > 


[T  -  cos2  0 


—  (  [iT-  k*^  cos^  <}>  -^  )  d9  d^) 


Performing  the  different at ion  and  rearranging  the  integrals,  this  can 


be  written 


<w,Ly>  =  -  f  - i -  f  |l  -  oos2  t  w 

ll  -  k'2  cos2  *  “  L 


k^  cos  0  sin  0  ay  it 

— : — . .  w  rr  do  d6  -  f  - L - 

|l  -  k^  cos2  0  |l  “  k*^  cos*^  0 


It - - - -  gcy  k»2  gog  ^  gji^n  ^  gy 

/  1  -  k'^  cos2  (j)  W  S-  +  . .  - - -  w  T—  i 


^  +  . .  . \  dQ  (a8) 

|l  -  k'2  cos2  (j>  J 


Integrating  the  first  term  of  the  first  0  integral  and  the  first  term 
of  the  second  (}i  integral  by  parts,  this  becomes 
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<  V,  L  y  > 


3y 

w  — 

ae 


e=Tr 

6«0 


-  /  1 1  - 
0 


cos2  e 


(— )(^) 
^ae^'ae' 


de 


d(J> 


1 


cos 
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For  both  the  Dirichlet  and  the  Nevunann  problems,  there  are  four 
possible  combinations  of  w  sind  y.  Either  they  are  both  even,  both 
odd,  V  even  and  y  odd,  or  w  odd  and  y  even.  For  all  of  these 
combinations  for  both  problems,  the  endpoint  contributions  are  zero, 
or  the  entire  integral  is  zero.  In  any  case,  <  w,  L  y  >  can  be 
written 


<  w,  L  y  > 


lo 


[l  -  cos2  0 
jl  -  k*2  cos2  (j> 


(^)(^) 
ae  39 


1  -  k'2  cos^  <fi 
+ - 

1  -  k^  cos^  0 


a(j)' 


d0  d({) 


(AlO) 


It  is  easily  seen  from  the  symmetrical  nature  of  the  equation  with 
respect  to  w  and  y,  that  it  is  also  equal  to  <  y,  L  w  >,  and  therefore 

<w,  Ly>  =  <y,  Lw>  (All) 

and  L  is  a  self-adjoint  operator. 


In  order  to  determine  the  positive-definite  property  of  L, 


consider  equation  (Al),  Scalar  multiplication  by  y  yields 


<  y»  b  y  > 

<  y»  p  y  > 


(A12) 


where  y  is  the  eigenfunction  corresponding  to  the  eigenvalue  X. 

/  /  y^  (k^  sin^  e  +  k'^  sin^  if) 

<  y,  p  y  >  =  de  d^  (A13) 

1  -  k2  cos^  0  1  -  k*2  cos^  (}i 


f""  1  -  k^  cos^  6  9y  p 

<  y,  L  y  >  =  /  /  — {||)^ 

1  -  k*2  eos^  <{( 


1  -  k'^  cos^  (|)  ,3y,2 

- .  {-^r  d0  d.^ 

1  -  k2  cos2  0 


(Alll) 


As  a  point  of  interest  note  that 


<  y?  b  y  >  =  A 


(AI5) 


where  A  is  the  normalization  constant  employed  in  the  dyadic  Green’s 
function. 

Inspection  of  equations  (A13)  and  (Al4)  shows  that  the  integrands 
are  always  positive  and  consequently 


<  y»  b  y  > 

<  y.  p  y  > 


=  A  >  0 


(AI6) 
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Thus  the  operator  L  is  positive  definJ.te  and  \  is  always  greater  than 
zero.  Recall  that  A  =  v(v+l)  so  that  equation  (Al6)  implies  that 

V  i  0  (A17) 

or 

V  <  -1 

Note  that  equation  (Al)  is  symmetric  around  v  =-1/2,  chat  is,  the 
differential  equation  is  the  same  if  v  is  replaced  hy  -(vHj.).  In 
other  words,  for  v  =-0.U  or  -0.6,  or  for  v  =  0  or  -1,  or  for  v  =  a 
or  -(a  +  l),  the  differential  equation  is  the  same,  and  consequently 
the  solutions  are  the  same.  Thus,  the  eigenvalues,  v,  need  to  be 
considered  only  for  v  ^  -1/2,  or  for  v  ^  -1/2.  For  convenience, 

V  ^  -1/2  is  chosen.  As  a  result  of  this  choice  and  equation  (A17), 
only  positive  values  of  v  occvir  in  this  analysis. 
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APPENDIX  B 

DETERMINATION  OF  THE  EIGENVALUES 

It  is  stated  in  Chapter  III  that  the  eigenveilues  are  determined 
hy  simultaneously  solving  two  continued  fraction  equations.  The 
purpose  of  this  appendix  is  to  explain  how  this  is  accomplished. 

Consider  the  even  Dirichlet  problem.  In  order  to  have  solutions 
of  equations  (3. 21),  Jt  is  necessary  that  the  eigenvalues  satisfy 
equation  (3.26).  This  equation  is  an  infinite  continued  fraction 
containing  two  unknowns,  v  and  h.  For  a  derivation  of  this  equation, 
the  reader  is  referred  to  Ince's  paper[lO].  The  difficulty  here  is 
that  the  recurrence  relations  have  no  starting  point.  The  coefficient 
subscripts  approach  both  plus  and  minus  infinity.  Ince  has  managed  to 
determine  a  differential  equation  which  must  have  the  same  typ^  of 
solutions  as  the  ones  used  here,  but  with  solutions  that  can  be 
written  as  series  starting  with  Aq  and  going  toward  A„.  In  this 
way,  the  recurrence  relations  can  be  written  as  a  matrix  equation  euid 
a  tractable  determinant  can  be  identified.  Diagonalization  of  this 
determinant  leads  to  the  continued  fraction  equation  (3.26)  which  is 
written  here  as  equation  (Bl).  An  eigenvalue  equation  will  be  derived 
in  detail  later  to  illustrate  this  procedure. 
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1  +  (£v-l)  (2v+3)  k2/9 

h  = -  + - 

^  1  +  k^  -  ?jh/9 


16  (2v-3)  (2V+5)  k^/225  36  (2v-$)  (2vH7)  k^/1225 

+  1  +  k2  -  lth/25  +  1  +  k^  -  Uh/U9 


kr^  (2v-2i’+l)  (2v+2r4l)  k^/(Ur2-l)^ 


+  1  +  k^  - 


Uh 

■(2r+l)2 


and  y  is  given  "by 


y  =  h  -  v{v+l)  k^ 


The  conventional  notation  for  a  continued  fraction  expansion  is  used 
here,  i.e.,  the  term  following  the  lowered  plus  is  added  to  the 
denominator  of  the  preceding  term.  The  dummy  index  r  is  simply  an 
integer.  The  continued  fraction  terminates  for  v  equal  to  half 
integers.  For  v  =  1/2,  h  is  constant.  For  v  =  3/2,  the  equation  is  a 
second-order  polynomial  with  two  roots;  for  v  =  5/2,  it  is  a  third- 
order  polynomial  with  three  roots;  etc.  Actually,  there  are  an  infinite 
number  of  roots  for  each  value  of  v.  For  example,  for  v  =  1/2,  the 
equation  can  also  he  satisfied  hy  equating  the  entire  denominator  of 
the  second  term  to  zero;  for  v  =  3/2,  the  denominator  of  the  third 
term  is  zero;  etc. 

For  V  =  1/2;  h  =  0.375,  and  y  =  0.  This  gives  one  point  on  a 
plot  of  y  versus  v.  For  v  =  3/2,  the  lowest  root  is  h  =  1.009  and 
y  =  0.866.  If  values  of  y  can  he  determined  for  v  between  1/2  and  3/2, 


! 
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it  will  be  possible  to  plot  a  curve  of  y  versus  v  for  the  lowest  root 
of  equation  (Bl).  This  is  done  in  the  following  manner.  For  v  =  0.6, 
a  value  of  h  between  h  =  0.375  and  h  =  1.009  is  assumed.  The  first 
100  terms  of  the  continued  fraction  are  evaluated  for  this  (v ,  h)  pair. 
If  the  difference  between  the  left  hand  side  and  right  hand  side  of  the 
equation  is  greater  then  10“^,  the  value  of  h  is  increased  or  decreeused 
by  10“^  and  the  continued  fraction  is  evaluated  again.  This  is  repeated 
until  either  the  difference  changes  signs  or  the  magnitude  of  the 
difference  is  less  than  10”^.  When  this  occurs,  h  is  increased  or 
decreased  by  10“^.  The  evaluation  is  continued  until  the  difference 
changes  signs  or  the  magnitude  of  the  difference  is  less  than  10"^. 

When  this  occurs,  h  is  incremented  by  10“^.  This  time,  when  the 
difference  changes  signs  or  its  magnitude  is  less  than  10”^,  the  last 
value  of  h  is  assumed  to  be  correct.  The  value  of  y  is  computed  and 
another  point  is  added  to  the  plot  of  y  versus  v ♦  The  value  of  v  is 
then  increased  to  0.7  and  the  process  is  repeated  again.  By  incre¬ 
menting  V  in  increments  of  0.1  between  v=0  and  v=9  and  finding  y  for 
each  incrementation,  a  plot  of  y^  versus  v  is  obtained.  The  super¬ 
script  1  on  y  is  used  to  indicate  the  lowest  root,  and  the  subscript  0 
indicates  that  it  is  a  root  of  the  6  equation.  Table  8  gives  the 
values  of  y^  for  v  =  0  to  9j  =  1/2,  and  Figure  29  shows  the  curve. 

Curve  yf  is  obtained  by  starting  with  the  second  root  at  v  =  3/2. 
This  curve  is  computed  in  the  same  way  as  curve  y^.  Note  that  the 
second  root  of  equation  (Bl)  for  v  =  1/2,  which  can  be  obtained  by 
equating  the  denominator  of  the  second  term  of  equation  (Bl)  to  zero, 

can  8Q.SO  be  obtained  by  continuing  y^  back  from  v  =  3/2  to  v  =  1/2. 

0 
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This  is  mentioned  only  to  indicate  that  all  of  the  y's  which  are  roots 
for  each  value  of  v  can  he  obtained  by  using  this  method.  It  is  seen 
later  that  it  is  not  necessary  to  determine  all  of  the  roots  for  each 
value  of  V. 

2  . 

To  be  certain  that  Ug  is  the  curve  of  the  second  root  and  not  the 

third  or  fourth,  it  is  necessary  to  investigate  the  area  between  curves 

W*  and  The  same  procedure  is  used  to  find  h  for  a  fixed  value 

of  V.  It  is  found  that  there  are  no  values  of  h  lying  between  curve 

y^  and  y^  that  will  satisfy  equation  (Bl).  It  is  also  found  that  there 
6  0 

1  12 

are  no  roots  lying  beneath  curve  y.  and  thus  y;r  and  y  are  actually  the 

6  9  6 

curves  of  the  first  and  second  roots  of  equation  (Bl). 

The  rest  of  the  roots  are  found  in  the  same  way  and  are  tabulated 
in  Table  8  and  plotted  in  Figure  29* 

Some  computational  difficulties  arise  in  computing  the  higher 
order  roots.  The  continued  fraction  varies  rapidly  for  assumed  values 
of  h  close  to  the  actual,  value,  and  the  root  is  sometimes  difficult  to 
find.  For  example,  for  a  difference  in  an  assumed  value  of  h  of  as 
little  as  10”^,  it  is  possible  for  the  continued  fraction  to  change 
from  being  smaller  than  h  to  being  greater  than  h  and  back  again  to  an 
even  smaller  value  than  the  first.  V/hen  this  happens,  the  computer 
fails  to  find  the  root.  Whether  or  not  the  root  is  determined  in 
cases  like  this  depends  to  a  large  extent  on  the  initial  guess.  This 
accounts  for  the  blanks  in  Table  8  and  also  the  lack  of  four  place 
precision  in  some  of  the  roots. 

The  choice  of  100  terms  of  the  continued  fraction  was  somewhat 
arbitrary.  It  was  found  that  the  difference  in  accuracy  using  50  terms 
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instead  of  100  terms  was  usually  insignificant,  so  that  50  terms  is 
probably  sufficient;  however,  100  terms  were  used  tc  ‘nsure  confidence. 
The  extra  computer  time  for  100  terms  is  insignificant. 

The  next  step  in  de■^erminining  the  eigenvalues  is  to  solve  the 
continued  fraction  equations  generated  by  solutions  of  equation  (3.22). 
There  are  two  of  these  given  in  Chapter  III  by  equations  (3*32)  and 
(3.33). 

The  derivation  of  equation  (3.32)  will  now  be  given  as  an  example. 
All  of  the  other  eigenvalue  equations  used  in  this  paper,  except  the 
special  one  mentioned  at  the  beginning  of  this  appendix,  can  be  derived 
using  the  procedure  outlined  here.  The  recurrence  relation  is  given  by 
equation  (3.30).  With  some  simplification  in  notation,  this  set  of 
equations  can  be  written  in  matrix  form  as  follows: 
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In  order  to  have  a  non-trivial  solution  it  is  necessary  that  the 
determinant  be  equal  to  ssro.  Thus  the  eigenvalue  equation  is 
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For  the  purpose  of  the  derivation  assume  that  the  matrix  in  equation 
(b4)  can  be  truncated  as  a  7  x7  matrix.  Starting  at  the  lower  right 
hand  comer,  use  standard  techniques  to  di-agonalist  the  matrix.  First, 

at" 

multiply  column  seven  by  and  subtract  from  column  six.  This 

°12 

creates  a  zero  in  colimin  six,  row  seven.  Kow  multiply  row  seven  by 
c  10 

;r^  and  subtract  from  row  six.  This  createri  a  zero  in  colximn  sev^n, 
^12 

row  six.  The  matrix  now  has  the  form 
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Note  that  the  upper  left  6x6  matrix  hsis  the  same  fona  as  the 
previous  7x7  matrix.  Repeat  the  process  on  this  matrix,  and  on  the 
succeeding  ones  Tintil  the  entire  matrix  is  diagonalized.  The  final 
form  of  the  matrix  is  diagonal  with  each  tens  a  continued  fraction. 

The  term  in  the  upper  left  hand  comer  is 

aU  C2  *6  %  ag  C£, 

-  ■bii  -  '^6  -  bg 

^12  ^10 

~  ^12  _  .  .  .  •  (B6) 
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a.2  CQ 


^10  ^8 
-  ■^10 


The  second  term  is  just  the  denominator  of  the  second  term  of 
equation  (b6).  It  is 


^2 


%  °2  H  %  ^8  C6  no  ^3 


bi^  -  b(S  -  b8 


bio 


^12  °i0 
^12 


(B7) 


The  rest  of  the  terms  follow  this  same  pattern,  i.e. ,  the  third  term 
is  just  the  denominator  of  the  second  term  of  equation  (B7),  etc. 
Although  the  matrix  was  truncated  for  the  purpose  of  the  derivation, 
the  pattern  of  the  continued  fraction  for  the  infinite  matrix  is 
obvious  from  the  above  equations.  The  eigenvalue  equation  is  obtained 
by  requiring  that  the  first  terra  be  equal  to  zero.  Writing  out  the 
a^,  bjjj,  and  Cp,  terms,  rearranging  and  simplifying,  the  eigenvalue 
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equation  can  be  written  as  equation  (3.32)  which  is  written  below 
as  equation  (B8). 

-  (v2-i)(v)(v+2)  k’V8  (v2.9)(v-2)(v+1*)  k'V(l6)2 

n  =  - - - 

2  -  k'2  -  nA  -  2  -  k'2  -  n/i6 

(v2-23)(v-lt)(vt.6)  k'V(48)2 

2  -  k'2  -  r»/36  -  .  .  .  . 

(v2-(2r-l)2)( v-2r+2 ) ( v+2r )  k ’ *♦  / ( 8r ( r-l ) ) ^ 

2  -  k'2  -  riAr2  -  .  .  .  .  (b8) 


and  y  is  given  by 


y  -  1/2  (-n  +  v(v+l)  k'^) 


(B9) 


The  continued  fraction  terminates  for  integer  values  of  v.  Using 

the  integers  as  starting  points,  curves  of  y  versus  v  are  computed 

using  the  same  procedure  as  before.  These  curves  are  labeled  y^, 

y?»  Vxi  ebc.  The  roots  for  k'2  =  1/2  are  tabulated  in  Table  9  and 
<P  <P 

the  curves  are  dra>m  in  Figure  29.  Equation  (3.33)  is 


n  =  v(v+l)  k'2/2  +  2 


k.2  _  (v^-^)(v-l)(v+3)  k'V(6)2 

2  -  k’2  -  n/9 


(v^-16)(v-3)(v+5)  k'V(30)^ 
2  -  k'2  -  n/25 


( v2-ltr2 )  ( v-2r+l )  ( v+2r+l )  k '  ^  /  ( 2  ( Ur2-1 ) )  2 


2  -  k'^  -  n/(2r+l)^ 


(BIO) 


and  ji  is  given  by  equation  (B9).  This  continued  fraction  also  termi¬ 
nates  for  integer  values  of  v.  The  roots  are  determined  as  before  and 

are  tabulated  in  Table  9  and  the  curves  are  drawn  in  Figure  29.  These 

2  4  6  ^ 

curves  are  labeled  m.,  p.,  p.,  etc. 

9  9  9 

The  intersections  of  the  6  and  (J)  curves  in  Figure  29  give  the 
(v,  p)  pairs  that  satisfy  both  equations.  These  ai'e  the  eigenvalues 
of  the  even  Dirichlet  problem  and  are  tabulated  in  Table  1.  Even 
though  the  eigenvalues  are  determined  graphically,  they  should  be 
quite  accurate.  Accuracy  is  achieved  b>  plotting  each  intersection  on 
a  scale  large  enough  to  use  the  four  I'lace  accuracy  of  p.  It  is  felt 
that  the  eigenvalues  ietermined  by  using  this  graphical  procedure  are 
correct  to  three  places  for  most  cases;  however,  accuracy  is  only 
claimed  to  +  5  x  10“3.  Better  accuracy  can  be  achieved,  if  desired, 
by  using  finer  increments  of  v. 

I’rom  the  pattern  formed  by  the  curves,  it  is  evident  why  eill  of 

the  eigenvalue  pairs  for  each  continued  fraction  are  not  needed.  The 

family  of  0  cuirves  is  moving  upward  and  the  family  of  9  curves  is 

moving  downward  as  the  order  of  the  root  is  increased.  The  area 

beneath  the  p^  curve  and  the  area  above  the  p}  curve  contain  roots 
0  9 

of  each  continued  fraction  but  no  roots  that  can  simultaneously  satisfy 
the  0  and  9  equations.  Thus  it  is  not  necessary  to  solve  the  continued 
fraction  equations  in  these  regions. 

It  is  seen  from  the  curves  that  the  number  of  eigenvalues  is 
increasing  as  v  increases.  If  the  integer  n  is  defined  as  it  is  in 
Chapter  III,  n  -  i/2  £  v  <_  n  +  1/2,  there  are  n  +  1  eigenvalues 
for  each  value  of  v.  This  agrees  with  the  number  of  eigenvalues  for 
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the  even  Lame  polynomials.  It  is  interesting  to  notice  the  slope  of 
the  curves.  The  slope  of  all  the  curves  is  decreasing  as  o  =  0  is 
approached.  Actually,  the  slope  will  he  zero  at  v  =  -1/2,  and  the 
curves  will  he  symmetrical  around  a  vertical  line  at  v  =  -1/2.  This 
is  due  to  the  symmetry  mentioned  in  Appendix  A. 

Depending  on  the  form  of  the  eigenfunctions,  there  are  two 
continued  fraction  equations  for  the  6  eigenvalues  of  the  odd  Dirichlet 
problem.  Equation  (3.*t3)  is 


n  =  -  v(v+l)  k^/2  +  2  - 


(v^-U)(v-l)(v«-3)  k  V(6)2 
2  -  k  2  -  n/9 


(v^-l6)(v-3)(vt-$)  k^/(30)^ 
2  -  k2  -  n/25 


( v^-Ur^ )  ( v-2r+l )  ( v+2ri-l )  k^  /  ( 2  ( Ur2-1 ) )  ^ 

2  -  k2  -  q/(2r+l)2 


(Bll) 


and  u  is  given  hy 

y  =  1/2  (n  -  v(v+l)  k2)  (B12) 

The  roots  are  labeled  y^,  y|,  y^,  etc.  Their  negative  values  are 

6  0  0 

tabulated  in  Table  10  and  the  ciirves  are  shown  in  Figure  30. 

Equation  (3.^*5)  is 


(v2-9)(v-2)(v+l»)kV(8)2  (v2-25)(v-li)(v+6)kV(‘t8)2 

n  =  4(2-k  ) - - 

2  -  k2  -  n/l6  2  -  k2  -  ,1/36 

( v2- ( 2r+l ) 2 ) ( v-2r ) ( v+2r+2 )  k^/ ( Ur ( 2r+2 ) ) ^ 


2  -  k2  -  n/(2r+2)2 


(B13) 
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and  n  is  given  by  equation  (B12).  These  roots  are  labeled 
Pg,  etc.  Their  negatives  are  tabulated  in  Table  10  and  the  curves 
are  shovn  in  Figure  30.  The  continued  fraction  equations  for  the  (|) 
eigenvalxxes  are  the  same  as  equations  (Bll)  and  (B13)  with  replaced 
by  k'^.  For  k^  =  k'^  =  1/2,  the  equations  are  the  same,  and  the  only 
difference  is  in  the  sign  of  p. 

p  =  1/2  (-n  +  v(v+l)  k‘^)  (BlU) 

Thus  the  family  of  <|(  curves  is  the  mirror  image  of  the  family  of  0 
curves.  The  eigenvalues  satisfying  both  sets  of  equations  correspond 
to  integer  values  of  v,  and  the  eigenfunctions  ai'e  Lame  polynomials. 
There  are  n  eigenvalues  for  each  v  =  n,  so  that  the  Dirichlet  problem 
has  a  total  of  2n  +  1  eigenvalues  for  each  value  of  v.  This  is 
indicated  by  the  grouping  of  the  eigenvalues  in  Table  1. 

Accuracy  is  better  for  the  odd  Dirichlet  problem  than  for  the 
even  Dirichlet  problem,  since  it  is  known  that  the  intersections  occur 
when  V  is  an  integer,  and  the  eigenvedue  p  has  been  computed  for 
integer  values  of  v.  Thus  no  graphical  interpolation  is  necessary. 

This  remark  also  applies  to  the  eigenvalues  of  the  even  Neumann 
problem. 

The  even  Neumann  problem  is  also  similar  to  the  odd  Dirichlet 
problem  in  that  the  continued  fractions  for  the  0  and  ({>  solutions  are 
the  same.  The  continued  fractions  for  the  <}>  equation  are  equations 
(b8)  and  (BIO) which  have  already  been  solved.  The  curves  for  the  0 
equation  are  the  mirror  image  of  the  <i>  curves.  The  roots  of  the  ({i 
equations  and  the  negative  roots  of  the  0  equations  are  tabulated  in 
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Table  9*  The  curves  are  shown  in  Figure  31.  The  eigenvalues 
correspond  to  v  =  n,  an  integer;  and  there  are  n  +  1  y*s  I'or  each  v. 
They  are  tabulated  in  Table  3. 

The  continued  fraction  for  the  6  equation  of  the  odd  Neimann 
problem  is  given  by  equation  (Bl).  This  equation  has  already  been 
solved  and  the  roots  tabulated.  The  curves  are  shown  in  Figure  32. 

The  continued  fractions  for  the  <j)  equations  are  equations  (Bll)  and 
(BI3)  with  replaced  by  k*^.  The  roots  are  tabulated  in  Table  10, 
and  the  curves  are  shown  in  Figure  32.  The  eigenvalues  are  tabulated 
in  Table  3.  There  are  n  eigenvalues  for  eacli  v,  where 
n  -  1/2  1.  V  £  n  +  1/2 ,  makii a  total  of  2n  +  1  eigenvalues  for  each 
V  of  the  Neumann  problem. 

The  lowest  roots  of  the  even  Dirichlet  problem  Euid  the  odd 
Neumann  problem  for  k^  =  0.1  and  k^  =  0.9  were  also  determined.  The 
roots  are  tabulated  in  Tables  11,  12,  and  13,  and  the  curves  are 
drawn  in  Figures  33  and  31*.  The  eigenveJ.ues  are  given  in  Table  5. 

Values  of  v  and  h  for  the  cases  when  the  continued  fractions  in 
this  appendix  terminate  have  been  computed  by  Ince[ll]  and  Arscott[lU]. 
They  are  in  agreement  with  the  numbers  presented  here.  Ince  has 
also  discussed  some  properties  of  the  continued  fractions,  such  as 
convergence  and  asymptotic  characteristics.  The  convergence  property 
is  the  main  reason  for  the  use  of  the  variable  n  in  some  of  the 
continued  fractions.  The  reader  is  cautioned  that  Ince  has  used 
differential  equations  different  from  the  ones  used  in  this  paper.  A 
simple  change  of  variables  converts  his  equations  into  the  equations 
used  here;  however,  a  direct  comparison  of  the  eigenvalues  is  not  easy. 
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Eigenvalues  for  the  odd  Dirichlet  problem 
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Table  b  -  Lttenvalues  o:  the  0  Er,uellnn  :or  tbe  Even  Cirichlee  Proxies  end  the  Odd  rie-aern  Frorlea 


Ei(7envaluf9  ol  Uie  K'^ufttion  lor  the  Ev'*n  Dirichlet  <ind  n^usaxm  Problems.  (Eigenvalues  of  the 
0  Equation  for  the  t^en  Keusann  Proles  are  Civen  by  -p) 


Table  10 


Eigenvalues  of  the  5  Equation  for  the  Odd  Dlrichlet  and  Ne'Eiann  Problena.  (Eigenvalues  of  the  J.70 
0  Equation  for  the  Odd  Dlrichlet  frOblea  are  Given  by  -u) 


V  Ui  ;i2  U3  Ue  Us  Uo  Uv 

0.0  -  0.7178 

•3.1  -  0.6789 

0.2  -  0.6329 

0.3  -  0.9797 

0.4  -  0.5'.93 

0.5  -  0.4515 

0.6  -  0.3763 

0.7  -  0.2936 

0.8  0.2034 

0.9  -  0.1055 

l.O  0.0000 

1  •-  0.1133 

1.2  0.2346 

1.3  0.3633 

1.4  0.5011 

1.5  0.6465 

1.6  0.8003 

1.7  0.9624 

1.8  1.1330 


1.9 

1.3122 

-  1.6168 

2.0 

1.5000 

-  1.5000 

2.1 

1.6965 

-  1.3775 

-  1.245" 

2.2 

1.9020 

2.3 

2.1163 

-  1.1150 

2.4 

2.3396 

-  0.9747 

2.5 

2.5720 

-  0.8284 

2.6 

2.3136 

-  0.6757 

2.7 

3.0644 

-  0.5167 

2.8 

3.321*6 

-  0.3511 

2.9 

3.5940 

-  0.1790 

-  4.0173 

3.0 

3.8-'29 

0.0000 

-  3.8729 

3.1 

4.1613 

0.1858 

-  3.7245 

3.2 

4.4592 

0.3787 

•  3.5719 

3-3 

4.7667 

0.5787 

-  3.4150 

3-6 

5.0838 

0.7860 

-  3.2538 

3.5 

5.4105 

1.0008 

-  3.0882 

3.6 

5.7470 

1.2231 

-  2.9180 

3-7 

6.0931 

1.4530 

-  2.7433 

3.8 

6.4491 

1.69<« 

-  2.5639 

3.9 

6.8148 

1.9366 

-  2.3716 

-  7.3712 

4.0 

7.1930 

2.1903 

-  2.1904 

-  7.1904 

4.1 

7.5757 

2.4523 

-  1.9562 

-  7.0062 

4.2 

7.9709 

2.7226 

-  1,7968 

-  6.8187 

••.3 

8.3759 

3.0013 

-  1.5921 

-  6.6278 

4.4 

8.7909 

3.2835 

-  1.3820 

-  6.4337 

'•.i 

9.2158 

3.5843 

-  1.1664 

-  6.2361 

4.6 

9.6506 

3.3388 

-  0.9451 

-  6.0350 

•*.7 

10.0952 

4.2022 

-  0.7179 

-  5.8306 

4.8 

10.5499 

4.5244 

.  0.4848 

-  5.6627 

6.9 

11.0145 

4.8557 

-  0.2455 

-  5.4112 

-ll./lOl 

5.0 

11.4890 

5.1961 

0.0000 

-  5.1'.‘61 

-11.4891 

5.1 

11.9735 

5.5456 

C.25I9 

-  4.</n5 

-11.2651 

5.2 

12.4679 

5.9043 

0.5109 

-  4.7551 

-11.0379 

5.3 

12.9723 

6.2724 

0.7755 

-  4.5238 

-10.8a’/8 

5.‘< 

13.4867 

6.6498 

1.0475 

-  4.2-788 

-10.5747 

5.5 

14.0111 

7.0367 

1.3265 

-  4.0647 

-10.3386 

5.6 

14.5454 

7.4330 

1.6127 

-  3.8267 

-10.0996 

5.7 

15.0897 

7.8389 

1.9061 

-  3.5844 

-  9.8576 

5.8 

15.6440 

8.2544 

2.2070 

■  3.3379 

-  9.6127 

5.9 

16.2083 

8.6795 

2.5154 

-  3. "^70 

-  9.3650 

-17.0-150 

6.0 

16.7826 

9.1142 

2.8315 

-  2.8316 

-  9.1143 

-10.7828 

6.1 

17.3669 

9.5587 

3.1555 

-  2.5715 

-  8.8607 

-16.5175 

6.2 

17.9611 

10.0128 

3.4874 

-  2.3067 

-  8.6042 

-16.2492 

6.3 

18.5654 

10.4768 

3.8274 

-  2.0370 

-  8.3448 

-15.9780 

6.4 

19.1796 

10.9505 

4.1757 

-  1.7623 

-  8.0824 

-15.7039 

6.5 

19.8039 

11.4341 

4.5323 

-  1.4324 

-  7.8170 

-15.4268 

6.6 

20.4381 

11.9275 

4.8973 

-  1.1972 

-  7.5487 

-15.1470 

6.7 

21.0824 

12.4307 

5.2709 

-  0.9064 

-  7.2773 

-14.8643 

6.8 

21.7366 

12.9439 

5.6532 

-  0.6101 

-  7.0028 

-14. '.>789 

6-9 

22.4008 

13.4669 

6.044j 

-  0.3081 

-  6.7251 

-14.2908 

-23.3790 

7.0 

23.0751 

14.0000 

6.4442 

0.0000 

-  6.4443 

-14.0000 

-23.0753 

7.1 

23.7593 

14.5427 

6.8531 

0.3140 

-  6.1602 

-13.7065 

-22.7685 

7.2 

24.4535 

15.0954 

7.2711 

0.6344 

-  5.8727 

-13.4104 

-22.4587 

7.3 

25.1578 

15.6532 

7.6981 

0.9612 

-  5.5817 

-13.1117 

-22.1459 

7.4 

25.3720 

16.2309 

8.1344 

1.2945 

-  5.2872 

-12.8104 

-21.8302 

7.5 

26.5963 

16.8135 

8.5800 

1.6345 

-  4.9892 

-12.5066 

-21.5115 

7.6 

27.3305 

17.4061 

9.0340 

1.9813 

-  4.6874 

-12,2002 

-21.1900 

7.7 

28.0747 

18.0087 

9.4992 

2.3352 

-  4.3817 

-11.8913 

-20.8657 

7.8 

28.8289 

18.6212 

9.mo 

2.6963 

-  4.0722 

-11.5798 

-20.5387 

7.9 

29.5932 

19.2437 

10.4563 

3.0647 

-  3.7585 

-11.2658 

-20.2090 

-30.7131 

8.0 

30.3674 

19.8762 

10.9492 

3.4405 

-  3.4406 

-10.9493 

-19.8764 

-30.3678 

8.1 

31.1517 

20.5187 

11.4516 

3.8240 

-  3.1184 

-10.6302 

-19.5414 

-30.0194 

8.2 

31.9459 

21.1712 

11.9638 

4.2152 

-  2.7918 

-10  3086 

-19.2036 

-29.6680 

8.3 

32.7501 

21.8337 

12.4856 

4.6143 

-  2.4605 

-  9.9844 

-18.8634 

-29.3135 

S.4 

33.5644 

22.5061 

13.0172 

5.0215 

-  2.1245 

-  9.6576 

-18.5206 

-28.9561 

8.5 

34.3866 

23.1886 

13.5586 

5.4367 

-  1.7836 

-  9.3282 

-16.1753 

-28.5958 

8.6 

35.2228 

23.8811 

14.1097 

5.8603 

-  1.-1376 

-  8.9961 

-17.8276 

-28.2325 

8.7 

36.0671 

34.5836 

14.6707 

6.2923 

-  1.0865 

-  8.6613 

-17.4774 

-27.8664 

8.8 

36.9214 

25.2961 

15.2415 

6.7327 

-  0.7299 

-  8.3236 

-17.1249 

-27.4975 

8.9 

37.7856 

26.0186 

15.8221 

7.1818 

-  0.3678 

-  7.9832 

-27.1258 

9.0 

33.6598 

26.7510 

16.4127 

7.6396 

O.COOO 

-  7.6398 

-16.4127 

-26.7513 

Problem  for  -  0.1  and  =  0.9 


k2  = 

=  0-1 

k2 

V 

V 

V- 

0.0 

0.2372 

0 

0.1 

0.2318 

0.1 

0.2 

0.2253 

0.2 

0.3 

0.2179 

0.3 

-0.0480 

0.4 

0.2094 

0.4 

0.5 

0.2000 

0.5 

0.6 

0.1895 

0.6 

0.7 

0.1780 

0.7 

-0.4112 

0.8 

0.1653 

0.8 

-0.5417 

0.9 

0.1516 

0.9 

-0.6895 

1.0 

0.1369 

1.0 

-0.8552 

Table  12  ••  Eigenvalues  of  the  cp  Equation  for 
the  Even  Dirichlet  and  Neumann 
Problems  for  k’^  =  0.1  and 
k»2  =  0.9 


k»2  = 

=  0.1 

k»2 

=  0.9 

V 

1^ 

V 

1^ 

0.0 

0.0000 

0.0 

0.0000 

0.1 

0.0054 

0.1 

0.0369 

0.2 

0.0119 

0.2 

0.0825 

0.3 

0.0193 

0.3 

0,1377 

0.4 

0.0277 

0.4 

0.2038 

0.5 

0,0372 

0.5 

0.2821 

0.6 

0.04V7 

0,6 

0.3738 

0.7 

0.0592 

0.7 

0.4803 

0.8 

0.0717 

0.8 

0.6029 

0.9 

0.0853 

0.9 

0.7425 

1.0 

0.1000 

1.0 

0.9000 

APPENDIX  C 


DETERiMINATION  OF  THE  EIGENFUNCTIONS 

The  recurrence  relations  for  the  6  solutions  of  the  even 
Dirichlet  and  odd  Neumann  problems  are  the  same.  The  relationship 
is  given  by  an  infinite  set  of  three  term  equations,  the  set  extending 
toward  infinity  in  both  directions.  One  way  of  presenting  this  set 
of  equations  is 

•  •  o  •  • 


(Cla) 

a_lA_2+b_2.^-l‘*'  c^^Aq 

=  0 

(Clb) 

“^0  '^-1*  '’0  *0 

II 

0 

(Clc) 

a,  A 

+  b. 

A  +  c,  A  =  0 

(Cld) 

1  0 

1 

1  12 

a 

A  +b  A  +c  A=0 

(Cle) 

? 

1  2  2  2  3 

Since  there  is  no  starting  point,  the  evaluation  of  the  coefficient 
is  not  simple.  However,  a  continued  fraction  expression  similar  to 
the  ones  used  for  the  eigenvalues  can  be  determined.  For  example, 
equation  (Cld)  can  be  written 

^  -  ill  _ 

Al  ai  A1/A2 
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"but,  from  equation  (Cle),  k-^/k2  can  be  written. 
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Al  b2  c2/a2 

A2  a2  A2/A3 


(C3) 


This  can  be  repeated  as  many  times  as  desired;  however,  it  is  alreadi-- 
evident  that  the  infinite  continued  fraction  will  have  the  form 

^  ^  c3/a3  cl+/aU  ^ 

Al  "  ai  b2/a2  +  b3/a3  ^  blt/a^  +  b5/a5  + .  ^ 

It  is  also  possible  to  start  with  equation  (Clc)  and  go  in  the  other 
direction  to  determine  that 


~  -  ^12+  a_j^/c_i  a>2/c_2 

Aq  cq  t-i/c_i +b_2/c_2  ^-3/^-3 

a_3/c_3 

+  h^k/o-h  + .  (C5) 

By  equating  (C4)  and  the  inverse  of  (C5)  a  relationship  between 
the  eigenvalues  v  anc.  can  be  obtained.  This  is  in  essence  what  was 
done  in  order  to  determine  the  eigenvalues.  The  eigenvalue  equations 
are  actually  consistency  equations,  since  they  must  be  satisfied  in 
order  to  have  a  consistent  relationship  between  the  coefficients  of 
the  eigenfunctions. 

The  coefficients  a^,  bg,,  aind  Cm  are  functions  of  v,  y,  and  c*. 
Since  the  eigenvalues  have  been  determined  already,  am,  bm,  and  c® 
are  known  quantities  and  it  is  only  necessary  to  calculate  the  value 
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of  the  continued  fraction  in  equation  {Ck)  tc  determine  Aq/Ai.  For 
convenience  Aq  can  be  set  equal  to  unity,  so  that  Ai  is  knovn;  next 
A1/A2  can  be  determined,  then  A2  is  knovm.  This  can  be  repeated  for 
as  many  coefficients  as  are  desired.  In  order  to  compute  A.^,  one 
could  use  equation  (Clc)  to  VTite 

^  ^  (C6) 

Ao  ac  Ao/Ai 

eind  use  the  information  just  obtained  from  the  computation  of  the 
coefficients  with  positive  subscripts,  or  one  could  use  equation 
(Clb)  and  generate  a  system  of  equations  like  equation  (C5)*  In 
either  case,  the  common  factor  is  Aq.  It  was  found  that  by  starting 
with  Aq  and  working  toward  coefficients  with  large  magnitude  subscripts, 
a  cumulative  error  was  generated.  This  occurs  because  the  coefficients 
with  small  subscripts  are  generally  larger  than  the  coefficients  with 
large  subscripts;  thus  a  small  percentage  error  in  Ai  might  be  very 
l6irge  relative  to  A5.  For  this  reason  it  is  better  to  start  with  An, 
n  large,  and  work  back  to  Aq.  This  is  the  way  the  eigenfunction 
coefficients  were  computed  in  this  paper. 

The  series  converge  rather  quickly,  so  it  was  decided  that 
approximately  20  teims  should  be  more  than  sufficient.  Initially 
All  is  assumed  tc  be  unity.  Then  Aiq/Ah  is  determined  and  thus  Aiq 
'.s  known.  This  is  repeated  for  Aq/AiQj  As/Aq,  etc.  until  Aq  is 
determined.  Next  A_ii  is  assumed  to  be  unity.  A_iq/A_ii  is  determined 
and  thus  A_iq  is  luiown.  This  process  is  repeated  until  Aq  is 
determined  again.  The  Aq  found  by  starting  at  A_ii  is  set  equal  to 


the  Aq  cedculated  the  first  time  and  all  of  the  negative  subscripted 
coefficients  are  scaled  accordingly.  Then  using  all  of  the 
coefficients  from  A_2^j^  to  the  eigenfunction  is  normalized  so 
“hat  it  has  unit  magnitude  or  \mit  slope  at  zero  argument  according 
to  whether  it  is  even  or  odd  respectively.  The  normalized  coefficients 
which  are  greater  than  5  x  10“^  are  tabulated  in  the  tables  in 
Chapter  III. 

The  continued  fractions  used  in  the  calculations  described  in 
the  previous  paragraphs  are 


where 


^-1  ^  ^  Cjh/ Sjj  ®m+l 

Ajn  aji  ^+l/ ®m+l  ^m+2/ ®m+2 


k  1* 

,  (l+m-l)2  ,k2  ^  v(v+l)  k2 

bm  =  (  “i; - ( Y  -  1)  + - - - +  p  ) 


k2  ,  (ltm+l)(4m+3)  , 

Cm  =  T"  '  - i - “  v(v+l)) 

4  4 


(C7) 

(C8a) 

(C8b) 

(C8c) 


Equations  (C7)  and  (C8)  are  just  the  set  of  equations  (3.25)  written 
in  a  form  suitable  for  calculating  the  coefficients.  Although  these 
equations  are  good  for  all  m,  they  are  used  only  for  m  ^  1,  i.e., 
they  are  used  to  calculate  only  the  coefficients  with  positive 
subscripts.  For  the  coefficients  with  negative  subscripts,  the  set 
of  equations  (3.25)  are  written  in  a  different  form. 
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Affl+l  ^  %/%  ®ra-l/°m-l 

Am  Cm  ^-l/Cffl-l  +  ^m-2/cin-2 

affl-2/cni-2 

+  'bm-s/'^m-S  +•••*.  (C9) 

where  e^,  hm*  and  Cm  are  the  same  as  before  and  the  equation  is  used 
for  m  £  -  1. 

Equations  (C7),  (C8),  and  (C9)  are  used  to  calculate  the 
coefficients  of  all  of  the  eigenfunctions  of  the  even  Dirichlet  and 
odd  Neumann  problems,  except  those  for  which  v  is  a  half  integer. 

For  these  cases,  either  a®  or  Cm  is  zero  for  some  value  of  m,  and 
either  equation  (C7)  or  (C9)  becomes  singular.  For  v  =  2n  +  1/2, 

n  an  integer,  it  can  be  shown  that  =  0,  which  implies  that  eJLl 

of  the  coefficients  with  m  less  than  -n-1  are  zero.  For 
V  =  2n  -  1/2,  it  can  be  shown  that  Ajj+i  =  0,  which  implies  that  all 
of  the  coefficients  with  m  greater  than  n  +  1  are  zero.  For  each  of 
these  cases,  one  of  the  infinite  continued  fractions  becomes  a  finite 
continued  fraction.  For  example,  for  v  =  7.5>  Am  =  0  for  m  >  i*. 

Tlius  it  is  assumed  that  Alj  =  1  and  equation  (C7)  is  used  to  find  A3, 

A2,  Ai,  and  Aq.  The  rest  of  the  calculation  is  the  same  as  before. 

The  recurrence  relations  for  the  (J)  eigenfunctions  of  the  even 
Dirichlet  problem  and  the  odd  Neumann  problem  are  not  the  same,  but 
the  method  of  determining  the  coefficients  is.  For  both  of  these 
problems,  the  eigenfunctions  are  represented  by  infinite  series,  but 
tmlike  the  previous  case  the  summation  begins  at  m  =  0.  There  are 


V  ( \)+T  ^  Ic  ^  ^  It  *  ^ 

+  - -  p)  +  B2m+3  —  ((2m+3)(2m+2)  -  v(v+l))  =  0 

(Cllc) 

The  easiest  way  to  solve  this  set  of  equations  is  to  set  Bi  equal 
to  unity  and  solve  equation  (Clla)  for  B3.  Then  equation  (Cllh)  can 
h^  solved  for  B5,  etc.  As  was  the  case  with  the  0  eigenfunctions, 
starting  with  a  large  coefficient  and  working  toward  the  smaller 
ones  generates  a  cumulative  error.  By  using  all  of  the  recurrence 
equations  except  the  first  one,  ein  infinite  continued  fraction  can  be 
generated.  For  example,  if  equations  (Cll)  are  written  in  the  form 

bi  Bi  +  Cl  B3  =  0  ^  (C12a) 

a3  Bi  +  b3  B3  +  C3  B5  =  0  (C12b) 

a5  B3  +  b5  B5  +  C5  B7  =  0 


(C12c) 


the  equation  for  B1/B3  can  be  written 
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Bj  b3  ^  c3/r,3  c5/a5  c7/a7 

B3  a3  ^5/^5  *■  ^0/^,9  + . 

and  similarly  for  B3/B5,  B^/®- ,  etc^ 

Note  that  equation  (Cg*.,  can  be  written 

%  _  _  Cl 

B3  bi 


(C13) 


(Cll») 


Equation  (ClU)  can  be  equated  with  equation  (C13)  to  obtain  the 
eigenvalue,  or  consistency,  equation.  The  eigenvalue  equation, 
equation  (3.33),  was  actually  derived  by  using  standard  procedures 
on  the  coefficient  determinant.  Manipulating  the  determinant  and 
equating  (C13)  and  (Cllt)  are  essentially  two  different  ways  of 
handling  the  same  mathematical  problem;  the  final  result  using  either 
approach  is  the  same.  The  determinantal  method  is  more  straight¬ 
forward  and  systemmatic,  whereas  the  method  described  in  this 
appendix  lends  some  insight  into  why  the  eigenvalue  equation  must  be 
satisfied. 

From  equations  (Cll),  the  continued  fraction  for  the  coefficients 
is 

^2m-l  ^  _  ^2m-<-l  ^  ^2m+1^^2m+l 

B2m+1  a2in+i  b2ni+3/ 


^2m+3/^2m+3  ^2a+5^  ^-2m+^ 

■*’  ^2m+5/^2m+5  ^2m+T'^^2ra+7  . 

tn  _>  1  (CI5) 


where 


k'2 

a2m+l  =  ~T~  ((2m-l)(2!n)  -  v(v+l))  (Cl6a)  ■ 

*} 

^2m+l  =  ((2di+1)2  {  - - 1)  +  ^  -  u)  (Cl6h) 

2  2 

k'2 

C2m+1  =  "7“  ( (2m+3)(2m-.‘-2)  -  v(v+l))  (Cl6c) 

To  solve  for  the  coefficients,  is  set  equal  to  unity,  B39/Bi;i  is 
calculated,  and  thus  B39  is  known.  This  is  repeated  for  B3Y/B39, 
B35/B3Y,  etc.,  until  B^  is  knovm.  Then  using  till  of  the  coefficients, 
Bi  through  Blji,  the  eigenfunction  is  normalized.  Tabulated  values 
appear  in  Chapter  III. 

Convergence  properties  of  the  eigenfunction  expansion  caii  be 
determined  by  considering  the  continued  fraction  expression  for  the 
ratio  of  two  successive  coefficients,  such  as  An/Ajj_i,  as  n 
approaches  infinity.  The  subscript  n  can  be  m  or  2m,  depending  on 
the  type  of  eigenfunction.  Ince  has  investigated  this  and  determined 
that 

lim  An  o 

-  <  Yr  . 

n  -»■  ®  An-i 

He  concludes  that  the  series  do  converge,  and  that  the  convergence 
is  more  rapid  than  a  power  series  expansion.  For  a  more  complete 
discussion  of  this  subject,  see  Ince[lO,ll]. 

When  V  is  an  integer,  the  above  procedure  for  determining  the 
eigenfunction  coefficients  breaks  down,  because  either  ap^+i  or 
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c^m+i  is  zero  for  some  m.  For  integral  values  of  v,  the  eigen¬ 
functions  are  Lame  polynomials  and  the  series  are  no  longer  infinite. 
All  of  the  0  and  eigenfunctions  of  the  odd  Dirichlet  and  even 
Neumann  problem  fall  into  this  category.  Although  these  eigenfunctions 
are  written  in  eight  different  forms  (see  equations  (3.35)“(3.38)  and 
(3.5l)-(3.5it))  with  eight  different  sets  of  recurrence  relations,  the 
method  of  determining  the  coefficients  is  the  same.  A  typical 
recurrence  relation  is  given  hy  equation  (3.^2),  which  appears  below 

^2m-2  ~  ((2m)(2m-l)  -  v(v+l)) 

U 

.Aa.  ((2«)2  -1) . 

2  2 

+  A2ni+2  ~  ((2m)(2m+l)  -  v(v+l))  =  0 

*T 

Ao  =  0  m  ^  1  (C17) 

This  set  of  equations  is  very  similar  to  the  set  of  equations  (Cll). 

Tlie  first  equation  contains  two  coefficients,  A2  and  A4;  the  second 
equation  contains  Ag,  A^,  and  Ag;  and  the  rest  follow  a  three  term 
recursive  relationship.  The  main  difference  between  this  problem  and 
the  previous  problem  is  that  there  ere  only  a  finite  number  of 
coefficients.  When  v  =  9,  the  largest  value  considered  here,  the 
maximum  number  of  coefficients  is  four;  A2,  Al*,  Ag,  and  Ag.  These 
coefficients  are  foiind  simply  by  solving  each  equation  one  at  a  time 
starting  with  the  first.  With  a  simplified  notation,  the  set  of 
equations  (C17)  can  be  written  (for  v  =  9) 
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b2  A2  +  C2  Alj  =  0 

(Cl8a) 

ah  A2  +  bij  Ah  +  ch  Ag  =  0 

(Cl8b) 

• 

Ah  +  bg  Ag  +  c6  Ae  =  0 

(Cl8c) 

with  the  ajn's,  bm's,  and  Cm's  easily  identifiable  from  equation 

(CI7).  With 

A2  *  1 

02 

(C19a) 

alt  H 

(C19b) 

.0  -  “6  '>6 

*8--  ^A6 

(C19c) 

The  coefficients  are  then  normalized  and  tabulated.  The  fact  that 
the  eigenfunction  has  a  finite  nvmber  of  terms  is  indicated  in  the 
tables  by  the  statement  Am  =  0,  m  >  8.  There  are  seven  other  sets 
of  recurrence  relations  for  the  odd  Dirichlet  and  even  Neumann 
problem.  They  all  contain  a  finite  number  of  terms,  they  all  have 
integer  value*;  of  v  for  eigenvalues,  and  they  are  all  evaluated 
using  this  technique. 

The  same  problem  concerning  accuracy  exists  here  as  before, 
i.e.,  a  cumulative  error  is  generated.  The  problem  could  be  cast  in 
the  form  of  a  continued  fraction,  but  it  was  felt  that  this  would 
unduly  complicate  it.  The  ninnber  of  coefficients  involved  is  small 
enough  that  the  accumulated  error  is  usually  insignificant.  However, 


it  is  recommended  for  the  sake  of  accuracy  that  the  continued  fraction 
method  be  used  if  eigenfunctions  are  computed  for  v  >  9* 

Arscctt  has  tabulated  the  Lame  polynomials  for  the  values  of  v 

from  one  to  thirty,  and  for  =  0.1,  0.2,  . ,  0.9.  They  are  not 

in  the  same  form  as  the  Lame  polynomials  calculated  here,  but  they 
can  be  converted  with  the  help  of  some  trigonometric  identities. 
Several  of  his  polynomials  were  converted  and  compared  with  the 
functions  computed  in  this  paper.  The  comparison  indicated  good 
accuracy  for  most  cases,  however  it  did  show  some  evidence  of  the 
cumulative  error  problem.  For  v  =  9j  the  eigenfunction  contains 
foiu:  or  five  terms.  For  some  of  these  eigenf'jnctions,  the  magnitude 
of  the  coefficients  is  decreasing  as  the  subscript  is  increasing. 

For  these  eigenfunctions,  the  last  coefficient  usually  lacks 
precision;  however,  this  is  relatively  unimportant  because  it  is 
always  small.  This  lack  of  precision  is  a  consequence  of  subtracting 
two  large  numbers  which  are  nearly  equal. 


APPENDIX  D 


VARIATIONAL  METHOD  FOR  DETERMINING 
EIGENVALUE?  AND  EIGENFUNCTIONS 


It  vas  shown  in  Appendix  A  that 

<  y»  L  y  >  _ 

<  y»  p  y  > 


(Dl) 


where  L  is  the  two-dimensional  St-urm-Liouville  operator  defined  hy 
equation  (A2),  p  is  the  weight  function  defined  by  equation  (A3), 

A  is  an  eigenvalue,  and  y  is  the  corresponding  eigenfunction.  It 
was  pointed  out  in  Appendix  A  that  the  pertinent  eigenvalues  are  real 
and  positive  or  zero.  The  method  employed  here  is  essentially  the 
Rayleigh-Ritz  method  with  some  modifications  due  to  the  two-dimensional 
nature  of  the  problem  ani  due  to  coraputa-tional  difficulties.  Thus 


<  w,  L  w  > 

-  ^1 

<  w,  p  w  > 


(D2) 


where  w  is  any  function  in  the  domain  of  L,  and  A^^  is  the  lowest 
eigenvalue.  Equation  (D2)  is  the  basis  for  a  method  of  determining 
the  eigenvalues  and  eigenfunctions. 

Consider  the  following  eigenfunctions  of  the  odd  Neumann  problem: 


0(e)  =  -  Aq  sin  0/2  +  ki  sin  30/2  -  sin  50/2 

+  A2  sin  70/2  -  k_2  sin  90/2  +  . 

l8i» 


(D3) 
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<{>(()))  =  sin  <j!  +  B3  sin  3ij>  +  B5  sin  54>  + . 

w  (9,  41)  =  0(0)  4>((i>)  (D5) 

As  a  first  approximation,  let 

Bi  =  Ao  =  1  (D6a) 

A_^  =  A2  -  A_2  =  B3  =  B5  =  0  (D6B) 


Thus, 

w  (e,  41)  =  (-sin  0/2  +  Ai  sin  36/2)  sin  4>  (B?) 

Insert  this  into  equation  (D2)  and  evaluate  the  two-dimensionsil 
integrals.  The  result  is 

(A*.  T  +  Ai  U  +  V 

- - = -  >  Xi  (D3) 

(Ai)2  X  +  Ai  Y  +  Z 

where  T,  U,  V,  X,  Y,  and  Z  follow  from  the  integral  evaluation. 
Minimizing  the  left  hand  side  leads  to  a  quadratic  equation  in  A^. 

A  basic  assumption  in  this  method  is  that  Aq  =  1  is  the  largest 
coefficient.  Thus  if  either  or  both  roots  for  Aj  are  larger  than  one, 
they  are  discarded.  If  both  are  less  than  one,  the  smallest  is  used. 
Inserting  this  value  of  Aj^  into  equation  (D8),  the  first  approximation 
for  Ai  is  found. 

Next ,  assume  that 

w  (e,  4))  =  (-  sin  0/2  +  Ai  sin  30/2  -  A_i  sin  50/2)  sin  4)  (D9) 
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The  procedure  is  repeated  for  the  new  unknown,  A_^,  and  a  second 
approximation  for  is  found.  This  procedure  is  repeated  until  all 
of  the  coefficients  A^,  A_2^,  A2,  A_2,  and  are  known. 

Next  the  process  is  repeated  again  for  small  changes  in  the 
coefficients,  i.e.,  w  (0,  (|i)  is  assumed  to  have  the  form 

V  (9,  (|))  =  (-sin  0/2  +  (Ai  +  AAi)  sin  30/2 

-  sin  50/2  +  A2  sin  70/2  -  A_2  sin  90/2) 

(sin  41  +  63  sin  34>  +  B5  sin  54>)  (DlO) 

where  aA^  is  the  unknown.  This  first  order  correction  is  determined 
using  the  same  procedure  as  before.  It  is  repeated  for  A  A.  AA2, 

AA„2>  AB3,  and  AB5,  The  process  could  be  repeated  again  for  a  second 
order  correction  if  desired.  This  was  not  done  for  the  eigenvalues 
and  coefficients  of  the  eigenfxmctions  given  in  Table  lU. 

The  computation  of  other  eigenveilues  is  similar.  The  first 
eigenvalue  was  found  by  assuming  that  Aq  and  B-^  are  the  largest 
coefficients.  Other  eigenvalues  are  found  by  assuming  that  other 
coefficients  are  maximum;  for  example,  the  second  eigenvalue  in 
Table  l^i  results  from  assuming  A^  and  Bj  to  be  maximvim. 

The  main  difficulty  with  this  method  is  that  it  is  based  on 
assumptions  which,  in  turn,  are  mostly  based  on  intuition.  Thus 
there  is  no  assurance  that  all  of  the  eigenvalues  will  be  foiuid^  or 
even  if  all  of  the  numbers  are  actually  eigenvalues.  The  method 
works  quite  well  for  the  lowest  eigenvalues,  however,  as  can  be  seen 
in  Table  ll4. 
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The  coefficients  in  this  table  were  normalized  so  that  the 
maximum  values  for  both  sets  of  coefficients  agree. 
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Table  l4  -  Compaxison  of  Eigenvalues  and 
Eigenfunctions  Determined  by 
the  Variational  Method  and  the 
Continued  Fraction  Method 


Variational 

Fraction 

V 

0.815 

0.8l4 

Ao 

-1.473 

-1.473 

Ai 

0.215 

0.214 

A_i 

0.030 

0,030 

Aa 

0.007 

0.007 

A-a 

0.001 

0.002 

Bi 

0.964 

0.964 

B3 

0.011 

0.010 

Bs 

0.001 

0.001 

V 

1.597 

1.595 

Ao 

0.183 

0.174 

Ai 

0.689 

0.689 

A-i 

-0.028 

-0.026 

Aa 

•0.005 

-0.005 

A-a 

-0.001 

-0.001 

Bi 

1.160 

1.160 

B3 

-0.049 

-0.048 

Bs 

-0.002 

-0.003 

V 

1.997 

1.955 

Ao 

-1.360 

-1.360 

Ai 

0.581 

0.568 

0.186 

0.179 

Aa 

-0.018 

-0.016 

A.2 

0.004 

0.004 

Ba 

0.496 

0.496 

B4 

0.001 

0.002 

V 

2.806 

2,803 

Ao 

0.414 

0.346 

Ai 

0.718 

0.718 

A--, 

-0.209 

-0  180 

Aa 

-0.092 

-0.091 

A-a 

0.004 

0.003 

Ba 

0.583 

0.583 

B4 

-0.038 

-0.038 

V 

2.984 

2.990 

Ao 

-1.720 

-1.720 

Ai 

1.032 

1.025 

A-i 

0.427 

0.429 

Aa 

-0.105 

-0.104 

A-a 

-0.010 

-0.011 

Bi 

0 

0.032 

B3 

0.322 

0.322 

188 


189 


REFERENCES 


1.  Kraus,  L.,  "Diffraction  by  a  Plane  Angular  Sector",  Thesis,  New 
York  University,  1955- 

2.  Radlow,  J.,  "Diffraction  by  a  Quarter  Plane",  Archive  for  Rational 
Mechanics  and  Analysis,  Vol.  8,  (1961),  pp.  139-158. 

3.  Radlow,  J.,  "Note  on  the  Diffraction  at  a  Corner",  Archive  for 
Rational  Mechanics  and  Analysis,  Vol.  19,  (1965),  PP.  62-70. 

4.  Jones,  D.  S.,  "Diffraction  by  an  Edge  and  by  a  Corner",  Quarterly 
Journal  of  Mechanics  and  Applied  Mathematics  ,  Vol.  V.,  Ft.  3, 

(1952),  pp.  363-377. 

5.  Braunbek,  W.,  "On  the  Diffraction  Field  Near  a  Plane-Screen  Corner", 
IEEE  PGAP,  Vol.  AP-4,  (1956),  pp.  219-223. 

6.  Kraus,  L.  and  Levine,  L.,  "Diffraction  by  an  Elliptic  Cone",  Com¬ 
munications  on  Pure  and  Applied  Mathematics,  Vol.  XIV,  (I961), 

pp.  49-68. 

7.  Morse,  P.  and  Feshbach,  H. ,  Methods  of  Theoretical  Physics,  McGraw- 
Hill,  New  York,  1953,  p.  154,  659. 

8.  Byerly,  W. ,  An  Elementary  Treatise  on  Fourier's  Series,  Dover,  New 
York,  1893,  p.  263. 

9.  Moon,  P.  and  Spencer,  D.,  Field  Theory  Handbook,  Springer-Verlag, 
Berlin,  I96I. 

10.  Ince,  E.,  "Ft’rther  Investigations  into  the  Periodic  Lame'  Functions", 
Proceedings  of  the  Royal  Society  of  Edinburgh,  Vol.  60,  Sect.  A., 

(1940),  pp.  87-99. 

11.  Ince,  E.,  "The  Periodic  Lame  Functions",  Proceedings  of  the  Royal 
Society  of  Edinburgh,  Vol.  60,  Sec.  A.,  (1940),  pp.  47-63. 

12.  Erdelyi,  et  al.,  Higher  Transcendental  Functions  Vol.  Ill,  McGraw- 
Hill,  1955,  Ch.  15. 

13.  Prasad,  G.,  A  Treatise  on  Spherical  Harmonics  and  the  Functions  of 
Bessel  and  Lame',  Benares  Mathematics  Society,  Benares  Citv,  India. 

1930. 

14.  Arscott,  F.,  Tables  of  Lame'  Polynomials,  MacMillan,  New  York,  I962. 

15.  Morse,  P.  and  Feshbach,  H.,  Methods  of  Theoretical  Physics,  McGraw- 
Hill,  New  York,  1953,  Ch.  I3. 


190 

16.  Stratton,  J. ,  Electromagnetic  Theory,  McGraw-Hill,  1941,  Ch.  VII. 

17.  Van  Bladel,  J.,  Electromagnetic  Fields,  McGraw-Hill,  1964. 

18.  Spence,  R.  and  Wells,  C.,  "Vector  Wave  Functions",  Communications 
on  Pure  and  Applied  Mathematics,  Vol,  IV., (1951)}  PP«  95-104. 

19.  Jones,  D.  S.,  The  Theory  of  Electromagnetism  Mai.Mil‘'an,  New  York, 
1964,  Ch.  9. 

20.  Tai,  C.  T.,  "Radiation  from  Current  Elements  and  Apertures  in  the 
Presence  of  a  Perfectly  Conducting  Half-Plaiie  Sheet",  Technical 
Report  No.  45,  July,  1954,  Stanford  Research  Institute,  Stanford, 
California. 


Security  Clasiificetion 


DOCUMENT  CONTROL  DATA  -  R&D 

{S$euritf  elo*$lfieation  ofUtlt,  My  ofabsuaet  anj  indixi-if  annotation  must  be  entered  when  tht  ovcToll  report  is  classified) 

loiNATiNO  ACTIVITY  ElectroScleiice  Laboratory  hwort  sicuhity  classification 

partment  of  Electrical  Ensineerlno.  The  Ohio  UNCLASSIFIED 


I.  RCaORT  TITLI 

ELECTROMAGNETIC  DIFFRACTION  BY  A  PERFECTLY-CONDUCTING  PLANE  ANGULAR 
SECTION 


4.  DMCRIPTIVIt  NCril  (Type  of  report  and  inclusive  doUe) 

Scientific _ Interim _ 


8.  authors;  (Last  nans,  first  name,  initial) 


Ramon  S.  Satterwhite 
Robert  G.  Kouyoumjian 


«.  REPORT  DATE  . 

23  March  1970. 


so.  CONTRACT  OR  4RANT  NO. 

Contract  AF  19 (628) -5929 

Project,  Task,  Work  Unit  Nos. 
5635-02-01 


7a  TOTAL  NO.  OF  PAGES 


IKlt] 


DoD  Element 


61102F 


DoD  Subelement  681305 


to.  AVAILABILITY/LIMITATION  NOTICES 


fa  ORt:..NATOR'S  REPORT  HUMBERTS; 

ElectroScience  Laboratory  2183-2 
Scientific  Report  No.  2 


AFCRL-69-0401 


1  -  This  document  has  been  approved  for  public 

release  and  sale;  its  distribution  is  luilimited. 


II.  SUPPLEMENTARY  NOTES  «*•  SPONSORING  MILITARY  ACTIVITY 

Air  Force  Cambridge  Research 

TECH,  OTHER  Laboratories  (CRD),  L.G.  Hans com  Fie Ic 

Bedford,  Massachusetts  01730 


IS.  ABSTRACT  In  thi«  reports  the  «x«ct  solution  (or  the  •Itctromegnetic  field  diffrected  by  a  perfectly  conducting  plane 

angular  eector  ii  detnrmired*  The  problem  te  a  three>dimen»lonal  vector  problem  and  the  aolutlon  le  preiented 
in  the  form  of  a  dyadic  Creen'c  functioni  vhich  le  the  moet  general  fos^m  cl  eolution  poetible.  Thue  the  vector 
ftelde  at  well  ae  the  current  on  tho  eector  may  be  determined  for  any  given  eource  excitation* 

The  corner  angle  of  the  plane  angular  eector  le  arbitrary,  varying  between  zero  and  v.  Special  caeee  of 
the  plane  angular  eector  are  the  half  plane  and  the  quarter  plane.  To  fin*'  the  fielde  for  larger  anglee,  euch  ae 
the  corner  of  an  aperture,  Babinet'e  principle  can  ^  ueed. 

The  dyadic  Creen'e  function  ie  eompoeed  of  vector  wave  functione,  which  in  turn  are  compoeed  of  ecalar 
wave  functione*  The  problem  ie  eolved  in  a  ephero'conal  coordinate  eyetem.  In  thie  eyetem>  the  plane  angular 
eector  iaone  of  the  cookdlNale  aurfacea,  ao  the  eeparation  of  variablee  technique  le  ueed  to  find  the  ecalar  wav'  I 

functione*  Tliey  enneiel  of  epherical  Beeaere  functione  and  Lame^  functione.  The  Lame'' functione  are  aolutune 
of  two  coupled  differential  equatione*  Theee  equatione  are  eolvei'  for  the  epecial  caee  of  a  quarter  .'^lane  icat* 
terer*  Ihe  firet  192  eigenvaluee  and  eigenfunettone  arc  computed  end  tabulated. 

Tbt  fielde  a.nd  currente  cloce  to  the  tip  of  the  quarter  plane  are  preeented.  Theee  fielda  and  currente 
have  been  the  aubject  of  much  conjecture  by  eeveral  authore*  It  le  ahown  that  the  dominant  field  behaves  ae 

where  the  low<jet  value  of  the  eigenvalue  v  ie  0.29b.  The  far  fielde  for  infiniteeimal  dipole  eourcee 
very  cloee  to  the  ti]<  are  aleo  deteimined  and  eeveral  patterne  are  prenented. 

Aa  with  any  exact  eolution  of  a  complex  problem,  the  reaulte  are  not  eimple.  It  le  felt,  however,  that 
the  exact  eolution  obtained  here  laye  the  foundation  for  eubeequent  vork.  For  inetance,  it  ahould  be  ooe^ible 
to  uee  thie  work  to  determine  an  aeymptotic  approximation  and  thue  derive  a  ’'diffraction  coeffiricnt"  for  the 
tip.  Thie  would  be  very  ueefut  in  Ktller'e  '’Geometrical  Theory  of  Diffraction."  Without  any  additional 
work,  there  are  aufficient  numerical  leeulte  preeented  m  thie  report  to  determine  the  fielde  within  approxi* 
rnately  one  wavelength  of  the  tip  for  a:<y  eource,  or  the  fielde  everywhere  for  a  eource  within  one  wavelength 
of  the  tip. 


FORM 
•'•'IJANS* 


UNCLASSl^-IED 

Security  Classification 


UNCLASSIFIED 


Security  Ciasstficatton 


Diffraction 

Plane  angular  sector 

Quarter  plane 

Fields  and  Currents 

Dyadic  Green’s  function 

Vector  wave 

Functions 

Eigenvalues 

Lame ’ -functions 

Sphero-conal  coordinate  system 


INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Enter  the  tame  and  eddreaa 
of  the  contractor,  aubcontractor,  grantee,  Department  of 
Diefense  activity  or  other  organization  (corporate  author) 
iaauing  the  report. 

2a.  REPORT  SECURITY  CLASSIFICATION:  Enter  the  over¬ 
all  aecurity  claaaification  of  the  report.  Indicate  whether 
"Reatricted  Data"  ia  included.  Marking  ia  to  be  in  accord¬ 
ance  with  appropriate  aecurity  regulationa. 

26.  GROUP:  Automatic  downgrading  ia  apecified  in  DoD 
Directive  5200.10  and  Arnmd  Forcea  Induatrial  Manuel. 

Enter  the  group  number.  Alao,  when  applicable,  ahow  that 
optional  markings  have  been  used  for  Croup  3  and  Group  4 
as  authorized. 

3.  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letters.  Titles  in  all  cases  should  be  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  show  title  classification  in  all  capitals  in  parenthesis 
immediately  following  the  title. 


covered. 

5.  AUTHOR(S):  Enter  the  namefs)  of  authorfs)  as  shown  on 
cr  in  the  report.  Enter  last  name,  first  name,  middle  initial. 
If  military,  show  rank  and  branch  of  service.  The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  date  of  the  report'as  day, 
month,  year,  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  i.e.,  enter  the 
number  of  pages  containing  information. 

76.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  report. 

8a.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

86,  8c,  4  8d.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  such  ss  project  number, 
subproject  number,  system  numbers,  tssk  number,  etc, 

9a.  ORIGINATOR'S  REPORT  NUMBER(S):  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  he  identified 
and  controlledhy  the  originating  activity,  Thiu  number  must 
he  unique  to  this  report, 

96,  OTHER  REPORT  NUMBER(S):.  If  the  report  has  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  the  sponsor),  also  enter  this  number(s), 


10.  AVAILABILITY/LIMITATION  NOTICES:  Enter  anv  limi¬ 
tations  on  further  dissemination  of  the  report,  other  than' those 
imposed  by  security  classification,  using  standard  statements 
Auch  as: 

(1)  "Oualified  requesters  may  obtain  copies  of  this 
report  from  DDC.” 

(2,  “Foreign  atmouncemenl  and  dissemination  of  this 
report  by  DDC  is  not  authorized." 

(3)  ‘‘U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.  Other  qualified  DDC 
users  shall  request  through 

_  _  »* 

(4)  "U.  S.  inilitary  agencies  may  obtain  copies  of  this  * 
report  directly  from  DDC.  Other  qualiffcd  users 
shall  request  through 

■  »> 

report  is  controlled.  Quail- 
fled  DDC  users  shall  requeot  through 

If  the  report  has  been  furnished  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi- 
cate  this  fact  and  enter  the  price,  if  known.  ^ 

11.  SUPPLEMENTARY  NOTES:  Use  for  eddltiona' explana- 
tory  notess 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  fpoy- 
fng /orj  the  resesrch  and  development.  Include  address. 

13.  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 

summary  of  the  document  indleatlvs  of  ik,. 
thou  ‘ .  ' 


y  _ 

end  with  an  indication  of  the  mllirary  'se'duriryVUsrflcanon 
Ihjormatlon  In  the  paragraph,  represented  as  (TS),  (S), 
(W>  or  (U), 

There  is  no  limitation  on  the  length  of  the  obstrset.  How¬ 
ever,  the  suggested  length  Is  from  150  to  22S  words. 

WORDSt  Key  words  are  technically  meaningful  terms 
or  short  phresee  that  enaracterlxe  a  report  and  may  lie  used  as 
Index  entries  for  cstaloglng  the  report.  Key  words  must  be 
selected  so  that  no  security  classification  Is  required.  Identi¬ 
fiers,  such  as  equipment  model  designation,  trade  name,  mili¬ 
tary  project  code  name,  geographic  location,  may  be  used  as 
key  words  but  will  be  followed  by  en  indication  of  technical 
context.  The  assignment  of  links,  rules,  and  weights  is 
optionsl.  ^ 


UNCLASSIFIED 

Security  Classification 


